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EINE CHRIA, DARIN ICH VON MEINEN ACADEMISCHEN LeBEN 
UND Wandel Nachricht GEBE. 

’Bin auch auf Unverstddten gewesen, und hab’ auch studirt. Ne, studirt hab’ ich nicht, 
aber auf Unverstddten bin ich gewesen, und weifi von allem Bescheid. Ich ward von 
ohngefdhr mit einigen Studenten bekannt, und die haben mir die game Unverstddt 
gewiesen, und mich allenthalben mit hingenommen, auch ins Collegium. Da sitzen die 
Herren Studenten alle neben ’nander auf Bdnken wie in der Kirch’, und am Fenster 
steht eine Hittsche, darauf sitzt ’n Profeflor oder so etwas, und fuhrt iiber dies und 
das allerley Reden, und das heissen sie denn dociren. Das auf der Hittschen safl, 
als ich d’rinn war, das war ’n Magister, und hatt’ eine grofie krause Pariique auf’m 
Kopf und die Studenten sagten, dafl seine Gelehrsamkeit noch viel grofler und krauser, 
und er unter der Hand ein so capitaler Freygeist sey, als irgend einer in Frankreich 
und England. Mochte wohl was d’ran seyn, denn ’s gieng ihm vom Maule weg als 
wenn’s aus’m Mostschlauch gekommen war, und demonstriren konnt ’ er, wie der Wind. 
Wenn er etwas vornahm, so fieng er nur so ehen ’n bisgen an, und, eh’ man sich 
umsah, da wars demonstrirt. So demonstrirt er z. Ex. daft ’n Student ’n Student 
und kein Rinoceros sey. Denn, sagte er, ’n Student ist entweder ’n Student oder ’n 
Rinoceros; nun ist aber ’n Student kein Rinoceros, denn sonst muflte ’n Rinoceros auch 
’n Student seyn; ’n Rinoceros ist aber kein Student, also ist ’n Student ’n Student. Man 
sollte denken, das verstiind sich von selbst, aber unser einer weifi das nicht besser. Er 
sagte, das Ding „dafi ’n Student kein Rinoceros, sondern ’n Student wdre“ sey eine 
Hauptstiitze der ganzen Philosophic, und die Magisters konnten den Riicken nicht fest 
genug gegenstemmen, dafi sie nicht umkippe. 

Weil man auf Einem Fufi nicht gehn kann, so hat die Philosophic auch den andern, 
und darin war die Rede von mehr als Einem Etwas, und das Eine Etwas, sagte der 
Magister, sey fur jedermann; zum Andern Etwas gehor’ aber eine feinere Nas’, und 
das sey nur fur ihn und seine Collegen. Als wenn eine Spinn’ einen Faden spinnt, 
da sey der Faden fiir jedermann und jedermann fur den Faden, aber im Hintertheil 
der Spinne sey sein bescheiden Theil, ndmlich das Andre Etwas das der zureichende 
Grund von dem Ersten Etwas ist, und einen solchen zureichenden Grund miifi’ ein 
jedes Etwas haben, dock hrauche der nicht immer im Hintertheil zu seyn. Ich hatt’ 
auch mit diesem Axioma, wie der Magister ’s nannte, ubel zu Fall kommen konnen. 
Daran hdngt alles in der Welt, sagt er, und, wenn einer ’s umstofit, so geht alles iiber 
und drunter. 

Denn kam er auf die Gelehrsamkeit, und die Gelehrten zu sprechen, und zog bey der 
Gelegenheit gegen die Ungelahrten lofi. Alle Hagel, wie fegt’ er sie! Dem ungelahrten 
Pobel setzen sich die Vorurtheile von Alp, Leichtdornern, Religion etc. wie Fliegen auf 
die Nase und stechen ihn; aber ihm, dem Magister, diirfe keine kommen, und kam’ 
ihm eine, Schnaps schliig’ er sie mit der Klappe der Philosophic sich auf der Nasen 
todt. Ob, und was Gott sey, lehr’ allein die Philosophic, und ohne sie konne man 


keinen Gedanken von Gott haben u. s. w. Dies nun sagt’ der Magister wohl aber nur 
so. Mir kann kein Mensch mit Grund der Wahrheit nachsagen, dafl ich ’n Philosoph 
sey, aber ich gehe niemahls durch’n Wald, dafi mir nicht einfiele, wer dock die Bdume 
wohl wachsen mache, und denn ahndet mich so von feme und leise etwas von einem 
Umbekannten, und ich wollte wetten dafi ich denn an Gott denke, so ehrerbietig und 
freudig schauert mich dabey. 

Weiter sprach er von Berg und Thai, von Bonn’ und Mond, als wenn er sie hdtte 
machen helfen. Mir fiel dabei der Isop ein, der an der Wand wachfit (i.Reg 4,33); 
aber die Wahrheit zu sagen, ’s kam mir doch nicht vor, als wenn der Magister so weise 
war, als Salomo. Mich diinkt, wer was rechts weifi, mufi, mufi ~ sah ich nur n’mahl 
einen, ich wollt ’n wohl kennen, malen wollt ich ’n auch wohl, mit dem hellen heitern 
ruhigen Auge, mit dem stillen grofien Bewufitseyn etc. Breit mufi sich ein solcher nicht 
machen konnen, am allerwenigsten andre verachten und fegen. 0 ! Eigendiinkel und 
Stolz ist eine feindseelige Leidenschaft; Gras und Blumen konnen in der Nachbarschaft 
nicht gedeyen. 

Mathias Claudius, Asmus omnia sua secum portans, oder: 

Sammtliche Werke des Wandsbecker Bothen, Erster Theil. 

Wandsbeck, beim Verfasser, i‘~l‘~l4, pp. Qsgg. 
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Introduction 


Overview 


This dissertation treats questions about the definition of “simplices” inside Riemannian 
manifolds, the comparison between those simplices and Euclidean ones, as well as 
Galerkin methods for variational problems on manifolds. 

During the las t three years, the “Riemannian centre of mass” technique described by 


KARCHEra ( 1Q77 I has been successfully employed to defi ne the notion of a si mplex in 


a Riem annian manifold M of non-constant curvature by Rustamov ( 20ioll . [SANDEFil 
(1201211 and others. This approach constructs, for given vertices pi € M, a uniquely 
defined “barycentric map” x : A —>■ M from the standard simplex A into the manifold, 
and calls x(A) the “Karcher simplex” with vertices pi. 

However, the question whether x is bijective and hence actually induces barycentric 
coordinates on x(A) remained open for most cases. We show that under shape regu¬ 
larity conditions similar to the Euclidean setting, the distortion induced by x is of the 
same order as for normal coordinates: dx is almost an isometry (of course, this can 
only work if A is endowed with an appropriately-chosen Euclidean metric), and Vdx 
al most vanishes. The estimate on dx could have already been deduced from the work 
of JoST and KarcherI ( iq 82I1 . but it is the combination with the Vdx estimate which 
paves the ground for applications of Galerkin finite element techniques. 

For example, the construction can be employed to triangulate M and solve prob¬ 
lems like the Poisson problem or the Hodge decomposition on the piecewise flat sim- 
plicial manifo ld instead of M. This leads to analogues of the classical estimates by 
Dziuk (Iiq88 ) and subsequ ent authors in the field of su rface PDe’s (we only mention 
Hildebrandt et al. 2006I and Holst and Ste^|2012| at this point), but as no em¬ 
bedding is needed in our approach, the range of the surface finite element method 
is extended to abstract Riemannian manifolds without modification of the computa¬ 
tional scheme. Second, one can approximate submanifolds S inside spaces other than 
]R™ (for example, minimal submanifolds in hyperbolic space), for which the classical 
“normal height map” or “orthogonal projection” construction from the above-mentioned 
literature directly carries over, and the error term generated by the curvature of M is 
dominated by the well-known error from the principal curvatures of S. 

Apart from classical conforming Galerkin meth ods, there are other discretisation 
ideas, e. g. the “discrete exterior calculus” (dec, see HlRANlI 200 in which variational 
problems such as the Poisson problem or the Hodge decomposition can be solved 
without any reference to some smooth problem. Gonvergence proofs are less developed 
in this area, mainly because albeit there are interpolation operators from discrete 
fc-forms to these interpolations do not commute with the (differing) notions 

of exterior derivative on both sides. We re-interpret dec as non-conforming Galerkin 
schemes. 
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Results 


Let (M, g) be a smooth compaet Riemannian manifold. Concerning the simplex defi¬ 
nition and parametrisation problem, we obtained the following (for German readers, 
we also refer to the official abstract on page[TTTj): 


(1) For given points ... ,pn € M inside a common convex ball, we consider the 

“barycentric mapping” x : A ^ M from the standard simplex into M defined by 
the Riemannian eentre of mass technique. Its image s := x{A) is called the (possibly 
degenerate) n-dimensional “Karcher simplex” with vertices pi. If A is equipped with a 
flat metric g® defined by edge lengths d(pi,pj) < h, where d is the geodesic distance in 
{M,g), and if vol(A,g®) > ah"' for some a > 0 independent of h (“shape regularity”), 
we give a estimate for the differenee — x*g between the flat and the pulled-back 
metric of order h^, as well as a first-order estimate for the difference V® — ® 

between the Euclidean and the pulled-back connection Jb.iyl 

(2) We give estimates for the interpolation of functions s IR and s —>■ A, where N 
is a second Riemannian manifold (j 7 j 4 | |7.r5j) . 


(3) Starting from the already existing theo ry o f Voronoi tesselations in R iemannian 
manifolds by ILeibon and LetscherI (|2nnrll and IBoissonnat et 'al. ( 2011 ), we define 
the Karcher-Delaunay triangulation for a given dense and “generie” vertex set (EH). 


(4) Concerning the Poisson problem on the spaee of weakly differentiable real-valued 
functions H^(M, IR), weakly differentiable real-valued differential forms H^O^(M), and 
weakly differentiable mappings into a second manifold H^(M, A), we prove error es¬ 
timates for their respective Galerkin approximations ( jlo.l3[ |10.17[ |13.14j ) . The same 
method giv es estimates for t he Ho dge decomposition in H^f 2 *^(M) if appropriate trial 
spaces as in Arnold et al. (200^ are chosen (iio.i5j). 


(5) We give proximity and metric comparison estimates for the “normal height map” or 
“orthogonal projection map” between a smooth submanifold and its Karcher-simplicial 
approximation, which is the classical tool for finite element analysis on surfaces in R^, 
but this time for submanifolds inside another curved manifold (iii.3[ Iii.i 8 |) . 


(6) We show that the differential of a Karcher simplex’ area functional with respect 
to variations of its vertices is well-approximated by the area differential of the flat 
simplex (A,5®) with (;® as above (\li.iiii . 


Concerning the convergence analysis of discrete exterior calculus schemes for a simpli- 
cial complex: 

(7) We define a (piecewise constant) interpolation ik : —>■ from discrete 

differential forms to a subspace of which turns the discrete exterior derivative 

into a “differential” d : with Stokes’ and Green’s formula for sim- 

plicial domains. This reduces convergence issues for dec from simplicial (co-)chains 
to approximation estimates between the non-conforming trial space (P^^n^,d) and 
(H^n^,(i). We estimate the approximation quality of P^^ forms in ( jq.iqi [9T2bj ) 

and compare the solutions of variational problems in and (lio.2HlER)l . 
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Structure and Method 


All the thesis is divided into three parts, one of which introduces notation, the main 
constructions another, its applications to standard problems in numerical analysis of 
geometric problems and surfaces pde’s (changing the usual setting from embedded 
surfaces to abstract (sub-)manifolds) the third. Hav ing in mind that “the introduction 
of numbers as coordinates [...] is an act of violence" ( WeylIiozloL p. 90), we try to stay 
inside the absolute Riemannian calculus as far as possible. Our main tool are Jacobi 
fields, which naturally occur when taking derivatives of the exponential map and its 
inverse. Whereas the standard situat i on for estimates on a Jacobi field J{t) are given 
values J(0) and j(0), see e. g. JostI ( 2niiL chap. 5), we will deal with Jacobi fields 
with prescribed start and end value, which is convered by (fairly rough, but satisfying) 
growth estimates 16.61 and |1^.4[ 
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Symbol List 


We listed here those symbols which occur in several sections without being introduced 
every time. Symbols with bracketed explanation are also used with a different meaning, 
which then will be defined in the section. Where it is useful, we added a reference to 
the definition. 


M 

m 

9 

P 

R 

r 

d 

Xp 

X 

inj, cvr 
Co, Cl 
h 

I? 

Co.i 

C^.i 

manifold, Mg is shortcut for (M, g) 
dimension of M 

Riemannian metric on M 

parallel transport (beside in section [3]) 

curvature tensor of M (fiT^ 

Christoffel symbols flTTSi). Christoffel operator (ji.l4j) 

geodesic distance function in M 

gradient of i dfn. ■) (il.SSi) 

barycentric mapping 

injectivity and convexity radius ([57^ 

global bound for ||i?|| and ||Vi?|| resp. 

mesh size 

fullness parameter 

:= Co + hCi 

■■= Co,ii?-2 

i? 

n 

r 

< 

simplicial complex (}4TT|) 

dimension of ^ 

elements, facets, simplices 

(realisation operator for simplicial complexes, 

< up to a constant that only depends on n 

A 

e-i 

In 

1 

BriU) 

standard simplex, Laplace-Beltrami operator 

Euclidean basis vector 
= (!,...,!) GR" 
unit matrix 

set of points with distance < r from U 

d 

6 

d 

V 

D 

L 

differential, exterior derivative 
(exterior coderivative, Kronecker symbol) 
partial / coordinate derivative, boundary of sets 
covariant derivative 

covariant derivative along curves (except section [jj) 

(weak Laplacian, [277]) , curve length functional 

1 • \f^ 

canonical Euclidean norm of R" 

I- l 

II- I 

l-l 

II II 

pointwise norm on bundles induced by g, volume of sets 
pointwise operator norm (iTTTI) 
integrated (or supremum) pointwise g-norm ([2^ 
integrated (or supremum) pointwise operator norm 
operator norm in function spaces (110.31) 
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Qk 

Lr 

\/\/k,r 

Hfc 

H§ etc. 
Hi,o HO,i 
Hi,i 

Hi+i 

P 

X 

i if , i iff 

L^X etc. 

fc-times continuously differentiable functions 
functions whose r’th power is Lebesgue-integrable 
functions that have k covariant differentials in L’’ (pT^ 

:= (except section [T^ 

functions in etc. with vanishing trace on the boundary 

forms a with weak da or 6 a of class resp. 

forms a with weak da and Sa of class 
forms with weak da and 6 a of class 

polynomial forms (jg.Op, functions (jlO.,'lj), vector fields (jiS.yj) 
vector fields of class C°° 
differential fc-forms of class C°° 

diff. forms with vanishing tangential/normal trace on the boundary 
vector helds of class etc. 

S 

TM\s 

TS^ 

V 

TT 

n 

t 

$ 

submanifold 

vector bundle over S with fibres TpM 
normal bundle of S' in M 
normal on S in M 
projection 

projection onto normal part 
projection onto tangential part 

normal height map p i—>■ exp^ Z for normal vector field Z 
geodesic homotopy p i—>■ exp^ tZ 
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A. Preliminaries 


Let us briefly recapitulate basic notions and concepts of the concerned mathematical 
fields: Riemannian manifolds and variational problems on these, simplex geometry and 
simplicial complexes. For the quick reader with experience in numerics on surfaces, a 
short look on the simplex metric in barycentric coordinates and our definition 

of simplicial complexes (j^75j) might be of interest. 


1. Riemannian Geometry 


For this section, we will keep close to the notations of JoST ( 2011 ) and Lee ( iqqtII . — 
Let (M, g) OT Mg for short be an m-dimensional Riemannian manifold. We write {X, Y) 
or g{X^ Y) instead of g{X, Y) for X,Y G TpM, mainly to prevent the use of too many 
round brackets. Whereas charts map open sets in M into IR™, we will mostly use 
coordinates {U,x), i.e. maps x from open sets U C R"* into M that are locally 
homeomorphisms. 

Throughout this thesis, we apply Einstein convention for computations in local co¬ 
ordinates or any other upper-lower index pair. Only when it explicity helps to clarify 
our statements, we note the evaluation of a vector field X or the metric g at a specific 
point p G M as X\p OT g\p respectively. 


Tangent Bundle and Norms. Coordinates {U, x) around p give rise to a basis 
or shortly di of TpM, and a dual basis dx® on T*M. The tangent-cotangent isomorphism 
is denoted by b and its inverse by jj. The natural extension of g to T*M has coefficients 
g'^^ with g^^gjk = (Kronecker symbol). On higher tensor bundles, g also naturally 
induces scalar products by g{v i^) v, w (E) w) ■= g{v,w)g{v,w) and similar for covector 
and mixed tensors. The space of smooth vector fields is denotes as X, the spaces of 
smooth alternating /c-forms as With •, we denote the Euclidean scalar product in 

R®®. 

We will denote the norm on all these bundles simply by | • | or | • |g, because we do 
not see ambiguity here. However, it differs from the operator norm of a tensor denoted 
as I • I. Both are equivalent, || • ||g < | • 1 ^ < c|| • ||g with a constant c that o nly dep ends 
on the dimension m and the rank of the tensor ( Golub and VAN LoanI Iiq 8 :: 

2.2-9). particular, operator and induced norm agree on 1-forms. 


eqn. 


Curvature 

In local coordinates {U,x), the metric g is a smooth field of positive definite m x 1.2 
m-matrices over U. A connection V on Mg is given in local coordinates by some 
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A. Preliminaries 


Christoffel symbols = T% via 

Va. 5 , = Vx^ = {X^d^Y^ + X^Y^T';^)dk (1.2a) 

for vector fields X, Y around p with coordinates X = X^'di and Y = Y^di respectively. 
It naturally induces a connection on higher tensor bundles, e. g. on the bundle of linear 
maps A : TpAI —>• TpM, by (Vy^)(Vl^) = Vv(AW) — ACVyW). There is a unique 
connection that is symmetric and compatible with g, the Levi—Civita connection 
oi Mg, whose Christoffel symbols can be computed by 

^ij = g^^d^gu + digji - digij). (1.2b) 


1.3 The Riemann curvature tensor R of XIg is defined by 


R{X, Y)z = VxXyZ — VyXxZ — V[X,F]-Z’. 

In local coordinates, it has coefficients 

RU = + r”fcrL - rr^rl,-, R{d,, d,)dk = R^^di (i-sb) 

and obeys the following (anti-)symmetries: 


{R{X, Y)Z, W) = -{R{Y, X)Z, W), 

{R{X, Y)Z, IT) = -{R{X, y)lT, Z), (1.3c) 

{R{X,Y)Z,W) = {R{Z,W)X,Y). 


Let us agree that V and D bind weaker than linear operators, so DtAW as above 
always means Dt{AW), not {DtA)W = AW. 

1.4 Along smooth curves c : [a; 6] —> M,t 1—> c(t), any connection uniquely induces a 

covariant differentiation Dt along c by 

V(t) := DtV{t) = {V'^ + cWW%)dk. (1.4a) 

A geodesic is a curve with vanishing covariant derivative, i. e. Dtc = 0 or, slightly 
inprecise, VcC = 0 . In coordinates. 




= -Pc’V't 


(1.4b) 


(note that we use the symbol c only for the covariant derivative of c, and we denote 
the coordinate derivative by a comma-separated subscript). If the parametrisation does 
not matter, we denote a curve c with endpoints p,q £ M as c : p q. The geodesic 
distance d(p, q) is the length of the shortest geodesic p'^ q. A Riemannian manifold is 
complete if any two points can be joined by a geodesic. For some neighbourhood B of 
p, we say that B is convex i f each two points q,r £ B have a unique shortest geodesic 
q'^ r in M which lies in B ( KarcherI lofiSh . 

Along a geodesic c : [a; 6] — >■ M, there is a parallel translation 
Tc(i)M for every € [a; 6], defined by = W{t) for the vector field W along 7 
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1. Riemannian Geometry 


with W{s) = V and W = 0. Parallel translation is an isometry, as = {W, W) = 

0 . The derivative of P with respect to a variation of c is computed in |7.8| 

As geodesics are unique inside a convex ball B, we will also write for G B. 

The unintuitive order of the evaluation points is inspired by the fact that some vector 
in TpM enters on the right, and a vector in TqM comes out on the left.—We remark 
that in general ^ instead pP’’^p^’‘iP’S’P is the holonomy of the loop 

p q r p. 

Assumption. Throughout the whole thesis, we will assume that Mg is a compact 1.5 
smooth m-dimensional manifold (without boundary, if not specified) with curvature 
bounds ||P|| < Co and ||VP|| < Ci everywhere. To keep definitions together, we give 
a “forward declaration”: When a radius (or a mesh size) r and a fullness parameter -d 
are defined, we will also use Cq := C'o'd”^ and Cop := Cq + rCi, analogously Cq i {C[ 
will not be used). 

Remark. Up to a factor of the bound ||P|| < Co is the same as requiring that the sectional 
curvature is bounded, becaus e if all sectional curvatures are bounded by TA, then ||P|| < |A 
(IBuser and KARCHERlfioSil . 6.1.1), which is the usual assumption in the works of Karcher, 

JosT et al. Of course, on the other hand K < Co- 


Second Derivatives 


Let Njap and Mgij be two smooth Riemannian manifolds with coordinates and 
u® respectively and / : TV —>■ M be a smooth mapping. Its first derivative is, at each 
p € TV, a linear map dpf : TpN Tf(pfM. Of course, the Levi Civita connections of 
M and TV induce a unique way to define the Hessian Vdf. For this purpose, df has 
to be considered as a section in E := T*N ® f*TM, a bundle over TV with fibres 
Ep =T*N X Tf(jqM. We want to give a coordinate expression for this. 

Definition. Let M and TV be two Riemannian manifolds, f : N ^ M smooth. The 
Hessian of / is V^df, a section of T*N O T*N O f*TM. 


Fact. The connection on the cotangent bundle T*N is defined by 

d{uj{X)) = a;(V™A) + (V'^*“w)(A:) for w G H^(A), A G X(A), 
cf. JostI ||20iiL eqn. 4.1.20). This gives 


1.6 


0 = d{du‘^{dfi)){d^) = du^{Xa,dp) + (V9^du“)(9p) 

= du^{T^p^ds)P{Xa,du^){dp), 

and with du'^{ds) = 1 if a = i 5 and 0 else, this gives that maps a vector dp to 

-r^^, so 

= -r^^du^. (1.6a) 

Vector fields V on M pull back to vector fields f*V by {f*V)\p = F|/(p). The 
connection then induces a connection on f*TM by 

xf;™rv = rv%^v. (i.eb) 
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A. Preliminaries 


Let us abbreviate da := as before, and additionally di := f*-^, and /*„ := §^- 
For example, the usual coordinate representation of df is df{da) = / q ^- Therefore, 


d, = = rif\ 

“ J ,a ’ 


k _d 
ij dv 


rf.^) = f:ar^A. (i.6c) 


The connections on T*N and f*TM induce a connection on the product bundle, cf. 


JOSTI l|20iil eqn. 4.1.23): 


V^{u}<g>V) = (V^*^w)(g)l/+a;®(V^*™y) for to G n\N), V G f*X{M). (i.6d) 


1.7 Lemma. Let f : N ^ M be a smooth mapping between Riemannian manifolds, and 
let y, VF G TpN. Then consider a variation of curves 7(5, t) in N with dt'y = W, 
dsl = V o-nd Dgdt"/ = 0 (everything is evaluated at s = t = 0). Let c := / o 7 
be the corresponding variation of curves in M. Then dtc = dfV, dsC = dfW and 
{V^df){V,W) = Dgdtc. If dfV 7^ 0 this is 

{y^df){v,w) = vl^dfv, 

where V and W are extended such that VwV = 0. 

Proof. Inserting df = Padu°‘ 0 di in li.bdl we have 


Bv li.bal 


^0,df = Vll^ifladuP ®d, + Padu^ ® 




and together with 1 1.68 this gives 

^a,df = {f^pdu^^ - f:aT^p^duP ®d. + f:^du^ ® ppT'(^d, 

Icf. l.Tosilbnii . eqn. 8.1.19). We conclude that Vdf, taken as bilinear map TpNxTpN 
Tf(^p')M, acts on vectors d^ and dg as 


Vdfid0,ds) = iPapdu^^ids) - f:aT'^g^duPds)]d.+f:adu<^ids)fpT':^dk 

= if:sp-fAh + fAAA) 9 ^- 


This is, as it should be, symmetric in (I and 5 by the symmetry of P^g and the 
Christoffel symbols. 

On the other hand, let us compute Dgdtc. The derivatives of 7 are given by 9(7 = 
7')9 q, and 9^7 = 7)^9^. By the chain rule, dtC = c^ = 'y'ffadi and dsC = "fdj^^dj. By 

[iidal 

DAc = {AJ:c.),s+ttfA^sfA)k)d.. 
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1. Riemannian Geometry 


Now (7“/a),s = again by the chain rule. As we have assumed 

£>^9(7 = 0, we get 7“^ = -7,^7Sr^5 for every a, so 


Dsdtc= 


ps 


+f:o.ptti^s+iv^,Mp^]k)d^, 


= v^w^{-f%r^0s + flsp + 


q. e. d. 


Corollary 

/ : A^ —>■ R, applied twice to the tangent of a geodesic 7, is the second derivative of 
/ o 7, and it holds 

Vdf{V,W) = (Vvgrad/,H^) = {Vw f,V) = V{Wf) - df {VyW). (1.8a) 

Scalings 

In most situations, we will try to prove scale-aware estimates, i. e. estimates for co¬ 
ordinate expressions or absolute terms where both sides of the inequality scale similar 
when the coordinates or the diameter of the manifold is scaled (if both sides of the 
inequality even remain unchanged under rescaling, we call the estimate scale-invari¬ 
ant). Therefore, we will need to know the scaling behaviour of vectors and tensors. 


■Tos-df 


2011 


eqns. 4.3.48, 4.3.30). If M = R, then the Hessian of a function 1.8 


Coordinate change, fixed absolute manifold. First, consider the case where the 1.9 
abstract (absolute) geometry of Mg is fixed and only coordinates are changed. A useful 
application is when coordinates (U, x) are given and the eigenvalues of the matrix g“- 
lie between and but one would like to have eigenvalues in the order of 1 (i. e. 
between and 1 ). This is achieved by coordinates 


y = 


_d_ 

dy 



Components of vectors always scale like the coordinates: If VF = 
then = pw'^’^. This scaling indeed fulfills our requirements: 




= 9 {^ 


dy' 


dyo ) ^2 


— —) = —qA 

dx" dxi ! ,.2^i'3 

O' 


The inverse matrix obviously scales with (g®-’)^ = /r^(g®'l)“.—The Christoffel symbols 
and the components of the curvature tensor scale with 

(F^)" = i(^^)^ = ^( 4 ,)-. 

fl II 
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A. Preliminaries 


1.10 Fixed coordinates, manifold scaling. Consider a new Riemannian manifold Mg 
with g = g?g. Then diam(Mg) = /rdiam(Mg) and d(p, g) = /rd(p, g). The norm of 
a tensor that is covariant of rank k and contravariant of rank i scales with . For 
example, a vector W, a linear form w and the curvature tensor R scale with 

\W\g=g\W\,, |u;|, = ^|a;|„ ||R||, = ^||R||^. 

If coordinates {U,x) remain the same, then gij = g'^gij and g^^ = and the 

Christoffel symbols and tensor components remain fixed: 

Now suppose two manifolds Mg and Nj with a mapping f ■. N ^ M. Consider a 
scaling g for M and v for N. As df can be regarded as a linear form on TN, resulting 
in a vector in TM, it is natural that the norm of df and Vc?/ scale as 

\\dfhs = ^\\dfU,g, ||Vd/||^.^ = ^||Vd/||^.,. (i.iob) 

The scaling behaviour of ||i?|| is the reason why we never suppress curvature bounds 
as “hidden constants”. In fact, most of our results could be simply worked out in balls 
of radius 1, and their scaling behaviour could be recovered from the curvature bounds 
and the scaling behaviour of left- and right-hand side operator norms. 

1.11 Coordinate change with manifold scaling. It might also be useful to use coordi¬ 

nates for (M, gfg) where the components gij remain unchanged, for example because 
they had previously been normalised to have eigenvalues in the order of 1 . If a chart 
([/, x) is known, such coordinates are given by y’’ = /rx®, because vector components 
also scale as and then 

\W\l = w°'''’w^’'’gafi = 

as it should. The Christoffel symbols and curvature tensor components scale as 

If two manifold and Mgij are scaled with factors g, and v in this way, resulting 

in coordinate expressions = vu°‘ for N and y* = /xx* for M, then the coordinate form 
of/, which was a mapping f/u —^ C/x, becomes a mapping / : vUu —t gUx^v i—^ gf{v/i'), 
so by chain rule 

ri _ If ^ f* 

J ,OL ^ 

for the components in 

df = f%du^ 0 ^, df = f^dv^^ 0 ^ 

and Vdf = - f^^Xps + f^sf!'p^jk)du^ ® du^ ® ^ as in the proof of[T?7l 
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1. Riemannian Geometry 


The Exponential Map and Special Coordinates 

On a point p £ M (interior, if M has boundary), there is, at least for some small 
intervall [— e;e], a unique geodesic for each initial velocity X G TpM. As the geodesic 
equation I r.qbl is homogenous and the unit ball in TpM is compact, this is equivalent 
to the fact that for some small ball around 0 € TpM, the geodesic with initial 
velocity X £ exists on [— 1 ;!]. As the unit sphere in TpM is compact, there is 
some e that works for any direction X. The exponential map is defined to map 
Bg M, expp(A) := c^(l). From this mapping, normal and Fermi coordinates can 
be constructed . The form er construction can be found in every Riemannian geometry 
textbook (e. g. iLE^bonql p. 78). 

Normal coordinates around p £ M are coordinates {U,x) with x( 0 ) = p in which 
straight lines t ^ tv are geodesics (arclength-parametrised for \v\12 = 1), which implies 
5ij(0) = Sij, dkgijiO) = 0 and r'(j (0) = 0 for all i,j, k. 

Lemma. Any orthonormal basis Ei of TpM induces normal coordinates (B^^x) around 
p via X : (it^,... 1—>■ expp(u*i 5 i), where e must be so small that geodesics through 

p are unigue. 

Proof. By homogenity of the geodesic equation |i.4b[ the geodesic starting with initial 
velocity c{p) = V with V = Eiv'^ has coordinates 

c^(t) = tv^, so c^{t) = and = 0 

for all t in the definition interval of c. At the same time, c\^ = —c^cM^j. As both 
equalities must hold for every V £ TpM, this already implies F^^- = 0 . The correspon¬ 
dence between F^^ and dkgij is linear and of full rank, so the latter have to vanish, 
too, q. e. d. 

Corollary, dexpp = id at 0 € TpM, that means (io(expp)R = V. 

Proof. Consider a geodesic c starting from p with velocity V. As the differential 
operator d(expp) applied to V can be computed as tangent of this integral curve, 
d{expp)V = c( 0 ) = V, q. e. d. 

Observation. Then the metric pij in the parameter domain is 

gij\u = {dxei,dxei) = {du{expp)Ei, du{expp)Ej), (1.14a) 

where U = u^Ei. Likewise, the Christoffel operator F : {v,w) !->■ T^jV^w^dk (which 
is bilinear, but does not behave tensorial under coordinate changes) is computable 
as pull-back of the connection to the parameter domain: The coordinate expression 
XyW = -I- r(u,iy) given in ll.sai can be understood as pull-back ® of the 

connection onto R'" (not to be confused with the connection x*V^ on x*TM from 
li.hbll . and such a pull-back is defined by dx{Vf, ^w) = Xdxvdxw. The right-hand side 
was identified to be Vdx{v,w) in [T?7l and so we have 

d£/(expp)(r(ei,ej)) = Xdu{expp){E„ Ej). (1.14b) 


1.12 


113 


1.14 
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A. Preliminaries 


1.15 Jacobi Fields. Let c{s,t) be a smooth variation of geodesics t 1—>■ c(s,t). Denote 
T := dtC = c and J := dgC. As T and J are coordinate vector fields, \J,T] = 0 , so 
VjT = VtJ or, in other words, Dgdtc = DtdsC. Differentiating the geodesic equation 
DtC = VtT = 0 gives, by[T^ 

0 = VjVtT = Vt^jT + i?( J, T)T = Vt^tJ + R{J, T)T, 

which is the defining equation for Jacobi fields: 

J = R{T,J)T (i-isa) 


Conversely, every vector field J along c fulfilling |i.i5a| gives rise to a variation of 
geodesics by 

c{s,t) := expe,,p,j(o) t(Pc( 0 ) + sP jjO)) (i-isb) 

where P is the parallel transport from c( 0 ) to expsJ(O) ( JostI Is 
P roposition Icf. IKarcherI [TqSqL eqn. 1.2.5). ^6^ c : I 
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thm. 5.2.1). 


1.16 Proposition Icf. IKarcherIIiq eqn. 1.2.5). c \ I ^ M be a smooth curve 
and Z be a vector field along c. Then the map Lpt s ^ exp^^^^ tZ(s) has derivative 

= J{t) for a Jacobi field J with initial values J( 0 ) = c(s), J( 0 ) = Z{s). In 
particular, dy(expp)hF is the value J(l) of a Jacobi field along t 1—>■ expptV with 
initial values J( 0 ) = 0 and J( 0 ) = W. 

Proof. c{s,t) := is a variation of geodesics t c{s,t) for every fixed s, so 

dgC = (/jj is a Jacobi field, and the values for t = 0 are J( 0 ) = 5 sc(s, 0 ) = c(s), and 
j( 0 ) = Dtdsc{s, 0 ) = Dsdtc{s, 0 ) = DgZfs) again by fTTT^ q. e. d. 


1.17 Fermi Coordinates. Let c : ]a;&[ — )• M be an arclength-parametrised geodesic. 

Then Fermi or geodesic normal coordinates along c are an open neighbourhood U 
of 0 € and coordinates x :]a',h\x U —>■ M, in which x{t, 0 ) = c{t) and straight 

lines s 1—c(t) + sv with first component = 0 are geodesics (arclength-parametrised 
for \v\p = 1 ) perpendicular to c. This implies 

gij(t, 0 ) = Sij, rfj (t, 0 ) = 0 for alH G ]a; 6[. (1.17a) 

If c is not a geodesic, then one can still find coordinates with gij(t, 0 ) = Sij, but the 
Christoffel symbols cannot be controlled. In classical surface geometry, those are called 
“parallel coordinates” along c (we will not use them). 

1.18 Lemma. Let c : ]a; 6[ M be a geodesic in M. Any orthonormal basis E2 ,..., Em of 
c( 0 )^ induces Fermi coordinates along c by x : {t,u^,... ,u^) 1—>■ exp^.^^^ igPP^’^Ei). 

Proof, x is injective because the orthogonal projection onto c is well-defined in a small 
tube around c, and if a point q € M projects to c(t), then the connecting geodesic 
c{t) q determines the components u'^, ..., by use of normal coordinates c{t)'^ —>■ 
M. 

By definition of normal coordinates, the claim gtj = 6 ij and = 0 along (t, 0 ,..., 0 ) 
is clear for i,j, k > 2 . Because c is arclength-parametrised, 511 = 1 , the orthogonality of 
c and Ei at every c{t) gives gu = 0 for all i. Because is parallel, = VcPi = 0 

proves the vanishing of the remaining Christoffel symbols, q. e. d. 
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1. Riemannian Geometry 


Corollary. If P is the parallel transport along a geodesic in Mg running through 
p € M, then for any vector V € TpM and a vector field W around p, we have 
PS/yW = V pvPW, and for a vector field V along a geodesic t i—>■ c{f), the fun¬ 
damental theorem of calculus holds: 


t 

V{t) = P*’°V{ 0 ) + J P*’^V{r)dr. (1.19a) 

0 


The Distance and the Squared Distance Function 

The fo llowinK properties already occur in Karche^ I iq8cJ) and Jost and KarcherI 
( iq 82I1 . but so metimes only hidden inside thei r proo fs. For the same calculations in 
coordinates, see lAMBROSiO and Mantegazza ( iqq8 ). 

The geodesic distance d( • ,p) is a smooth convex function in some small neighbour¬ 
hood B of p, excluded in p itself. It therefore has a gradient Yp, and its length is the 
Lipschitz constant of d( • ,p), namely 1 everywhere. Additionally, 


0 = V{Yp, Yp) = 2 (VyTp, Yp) = 2 {Vy,Yp, V) for all F G T,M, q G B, 

by symmetry 1 1. Sal of the Hessian Vdd, so is autoparallel everywhere. The integral 
curves of Yp are hence geodesics emanating from p with dd(7) = 1, so d{j{t),p) = t 
for each such curve. On the other hand, d( •, p) is constant on the distance spheres of p, 
so Yp is perpendicular to them (Gauss Lemma). In normal coordinates (it^,..., it™) 
around p, we have d( • ,p) = \u\p and hence 

d( • ,p)Yp = u^di. 


Observation. Base and evaluation point can be reversed, and the vector field only 
changes sign: Yp\g = —P’^-^Yglp, because both are velocities of the arclength-parametri- 
sed geodesic p q or q p respectively. 

Lemma. In a small neighbourhood of p, 

Xp ■= grad i d^(p, •) = d(p, • )Yp 

is an everywhere smooth vector field, its integral lines are (quadratically parametrised) 
geodesics emanating from p, and expg{—Xp\g) = p, equivalently 

-Xg\p = PP'‘^Xp\g = (eXPp)"l(? 

for all q in a convex neighbourhood of p. Loosely speaking, one also writes this as 
PXp = exp“^p. 

Proof. Let c be the arclength-parametrised geodesic with c( 0 ) = p and c(t) = q. By 
definition of exp, we have exppC(O) = q, as well as c{t) = P*’°c(0) and c{t) = Yp\c(t) 
for all t by the Gauss lemma. The switch of base and evaluation point is justified by 

11.21L q. e. d. 


1.19 


1.20 


1.21 


1.22 
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A. Preliminaries 


1.23 Lemma. For V £ TqM, where q is in a convex neighbourhood of p, let J be the Jacobi 
field along p q with J( 0 ) = 0 and J(t) = V. Then 

VvXp = tJ (t), VyyXp = tDsJ(t). 

In particular, if V is parallel to Xp, then XyXp = V and VyyXp = 0 . 

Proof. Let s 1—>■ ( 5 (s) be a geodesic with ^( 0 ) = q and ^( 0 ) = V. Define a variation of 
geodesics by 

c{s,t) := expp {t{expp)~^6{s)). 

Then dtc is an autoparallel vector field and J := dgC a Jacobi field along t 1—>■ c{s,t) 
for every s with boundary values J(s, 0 ) = 0 and J(s, 1 ) = 5 {s). The t-derivative is 

dtc{s,t) = P*’°(expp)"^J(s) = P*’^Xp\s(^s) 

and hence J(t) = Dtdsc{ 0 ,t) = Dsdtc{ 0 ,t) = DsXp\c{o,t) = Xj(^t)Xp. Differentiating 
this once more gives the claim for the second derivative. If V is parallel to Xp, then 
use VyT = 0, q. e. d. 

Remark, (a) Variations of Xp with respect to the base point p will be considered in 

EH31 

(b) Analogously to (exp^)”^ = PXp, the derivatives of Xp and exp^ correspond: 
VyAp is the derivative of some Jacobi field with prescribed start and end value, 
whereas (iy(expp)tV = J(l) for a Jacobi field with J( 0 ) = 0 and J( 0 ) = W. 

(c) Although Yp is not differentiable at p, we have XXp = id at p, similar to do expp = 
id. 

(d) In the notation of iGROHS et al\ (l2niqh . our vector field Xp and its derivative are Xp\a = 
log(a,p) and Wp|a = V2 log(a,p). 


Submanifolds 

1.24 Extrinsic Curvature. For a smooth fc-dimensional submanifold S C M, we treat 
TpS as a linear subspace of TpM, denote the orthogonal projection TpM — >■ TpS as t 
and the projection onto the normal space TpS'^ as n. The bundle over S with fibres 
TpM is denoted as TM\s = TS(BTS^ (meaning a fibre-wise sum of vector spaces). The 
Weingarten map or shape operator with respect to a normal field z/ is Wi, := Vv, 
that means U 1—V[/y. The second fundamental form with respect to v is 

MU, V) := -{WM, V) = {VuV, v) (1.24a) 

because (y, V) = 0 and hence U{v, V) = 0 . In particular, ^u{U, V) is in fact tensorial 
in V, U and V. Sometimes 1 I(C/, V) := nVuV is also called the second fundamental 
form in the literature, although it is a bilinear map, not a form. If the orthonormal 
parallel normal fields Vk+i,... ,Vjn locally span TS-^, it holds 1 I(C/, V) = Vi ^)- 

The covariant derivative induced by is V'® = tV, hence II = V — V'®. 
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2. Functional Analysis and Exterior Calculus 


Generalised Fermi or Graph Coordinates. The “tubular neighbourhood theo¬ 
rem” states that a small neighbourhood IB^ ( 5 *) oi S is diffeomorphic to S x with 
an open er-ball C IR™”*^ around 0 ( Bredon iqq4 thm. II. 11.4). By explicitely 
constructing this diffeomorphism, upper bounds on e can be derived: For t G [ 0 ; 1 ], let 

: TS^ ^ M, (p, Z) 1-^ expptZ. (1.25a) 

If p moves with velocity p, we know bv li.ibi that d^t{p) = J{t) for a Jacobi field with 
J(0)=:P, j(0) = VpZ. 

The case where Z = v is parallel along p is particularly interesting. Then d>t 
parametrises the level sets of the distance function from 5 ”; we have J( 0 ) = W^p, 
and the parallel transport of i/ along t is normal to the image of $4, so the whole curve 
fulfills J = VFi/J. Therefore, the pull-back m etric <h*q(iP , jp) ch anges with respect to t 


as -^^^g{p,p) = 2 g{J,j) = 2 g{J,W„J), see KARCHEra ( iqSqL eqn. 1.2.7). Hence the 


maximal eigenvalue of Wi, over v G C TS-^ and over t bounds e. We will pursue 

this more explicitely in lii .81 By J = -^{WJ) = WJ + Wj = ^J + W^ J, one then 
obtains a Riccati-type equation for the Weingarten map ( Karcher iq8qL eqn. 1.3.1) 


W = R^-W‘^ for = R{v, ■ )v. 


(1.25b) 


Generally, a tangent vector U G T(^p^z)TS-^ is induced by a curve s 1—>■ exp^j-^,) tZ{s), 
where p is tangential to S and Z = tZ + nZ. The above-mentioned Jacobi field J can 
be split into two Jacobi fields Jp{s) J- Ju{t) with initial values 


Jp( 0 ) =p 
jp(0) = tZ 


T40) =0 

i,(0) = nZ. 


(i-25c) 


The part tZ is in fact tVpZ (if we assume Z to be extended parallel along t), so it is 
uniquely determined by p, and thus U has the representation (p, nZ) in the chart <I>t. 
Let -0 be the orthogonal projection Be (S') —>■ S. As Jacobi fields with orthogonal initial 
values and velocities stay orthogonal, we have an orthogonal splitting V = Jp J- W for 
V G TpM, p G Be(S), with Vp = Jp(l), W = Jy{l). This gives a simple representation 
of namely d^{Vi,) = 0 and dij^^Vp) = p. The geometric interpretation of the splitting 
is 

Vp = ppMp) 14 = ppMp) nP'l’<^p)’PV, (i-25d) 

that means Vp and I4 are the orthogonal projections onto PTS and PTS-^ respectively. 
This is proven by -^{Jp,Z) = 0 (if Z is extended parallel along t) and the initial 
conditions {Jp{ 0 ),Z) = 0 and {jp{ 0 ),Z) = 0 . 


2. Functional Analysis and Exterior Calculus 

We will quickly review the Dirichlet proble m and the Hodge decomposition in this 
section. All proofs are reformulations from Schwar 3 (jiQQtJ) . but we tried to take 
special care that not the vector bundle structure of , but only its functional analytical 
nature has been use (the only exception is l2.i6ll . 


1.25 
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A. Preliminaries 


Notation, (a) We have defined X and as the spaces of smooth vector fields and 
fc-forms on M. The pointwise scalar product g on all tensor products of TM and T*M 
naturally induces an product on them: 

{v,w) := J g{v,w) 

M 

The completion X with respect to the norm will be called L^X, analogously for 
the differential forms. The notation (•, •) will only be used for the scalar product, 
so all indices like (•, • )l 2, (■, • )L2(Mg) ^ ‘ )l2x (’ > ’ only added for 

ease of reading. 

(b) Let M have a boundary dM. The projections t and n from TM\gM onto TdM 
and TdM-^ pull back fc-forms as ..., 14 ) = v{tVi ,..., tl4) and similarly n*uj. 

The spaces of fc-forms with vanishing tangential part on dM are called 17^. 

Together with the usual exterior derivative d, the form the smooth de Rham 
cochain complex 

^ ^0 

The exterior coderivative 6 is, for forms with appropriate boundary conditions, 
adjoint to d with respect to the scalar product: 

{v,dw)^_2Qk+i = {Sv,w)^_2Qk for all w e w e 

and all v € 


2-3 


The image and the kernel of d in 17^ are called the spaces of boundaries and cycles, 
:= imd|f2fc~i and := kerdjQfc. The space of harmonic forms is := fl 
d{Q'^~^)-^. For (5, we have and defined analogously, so fl bv [2TTI 

Denote 


Lap(u, ic) := {dv,dw) + {Sv,Sw), Dir(u) : 
Definition. Define the following six norms on each 


2 ._ 

HUO ■ — 

1?; 



\dv\l2 


2 ._ 
HO,i ■ — 

1?; 



\6v\l2 


2 ._ 
HUi ■ — 

1^ 



\dv\l2 + 

\dv\l2 = 

Hi+i : = 

1^ 


+ 

|d(5u|L2 + 

|(5(ir)|L2 

- 

1^ 


-h 

|Vu|2, 


:= 

1?; 



|V2uP,2 



Lap(r!, v). 


+ Dir(r;) 


(2.2) 


Let etc. be the completion of with respect to these norms. The L’’, W^’’’ 

and norms are the usual modification of the L^, and norms for exponents 
r ^ 2. 

2.4 Observation, (a) (H^’°D,(i) is a cochain and is a chain Hilbert complex, 

that means that d 01 5 are bounded linear operators with d^ = 0 or = 0 respec¬ 
tively (to be notationally precise, a Hilbert complex requires d or 6 only to be closed 
operators). 
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2. Functional Analysis and Exterior Calculus 


(b) The norm dominates the and the norm; the norm dominates all 
the other norms. For functions, and = H^. 

Remark. ISchwar 3 iIiqqfJ . sec. 1.3) uses a different definition of which depends on local 
choices of orthonormal brises. He then uses the and norms to control the exterior 
(co-)derivatives. We consider the use of a nd sim i lar no rms a sharper tool for this, as only 
the actually needed derivatives have to exist. Ijos 20 111 , eqn. 3.4.4) writes for what we 
call H^+h 


Fact ( Bruning and LeschIIiqq 21 corr. 2.6). The spaces are closed in 

if and only if are closed in If this is the case, and if is finite-dimen¬ 
sional, then is called a Fredholm complex. 


2.5 


Laplace Operator and Dirichlet Problem 
Observation. Directly from Green’s formulaone gets for v G w G 2.6 

Lap(u, w) = {{dS -I- Sd)v, w) 


in either of these four cases: 

t*w = 0 , n*w = 0 t*w = 0 , t* 5 v = 0 

n*w = 0, n* dv = 0 n*dv = 0, t* 6 u = 0 . 

Definition. The strong Laplacian is A := d 6 + Sd : L^O. The weak 2.7 

Laplacian is L : (H^’^D^)*,u 1— Lap( •,?;). The strong Dirichlet prob¬ 
lem is to find u G with 


Am = /, t*u = 0 , t* 6 u = 0 . 


(2.7a) 


The weak Dirichlet problem is to find u G with Lu = / in that 

means 

{du,dv) + {Su, 6 v) = {f,v) for all M G (2.7b) 

Such a M is called a Dirichlet potential for /. 

Fact (Dirichlet principle). A form u G is a solution of the weak Dirichlet 

problem if and only if it minimises Dir(M) — (/, v) over all m G H 


Remark. For sections of smooth vector bundles over M, the trace of the second 
covariant derivative gives a “metric” Laplace operator tr V^, co nnected to o ur Laplacian 
or “Laplace-Beltrami” operator by the Weizenbock formula ( Jost|[20ii . thm 4.3.3.), 


which we do not use, and only mention to avoid confusion. They agree if and only if 
Mg is flat. 


Proposition. If u G is a solution of the weak Dirichlet problem and is in 2.8 

addition contained in then it solves the strong Dirichlet problem. 
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Proof. Suppose Lap(w, u) = {f,v) for all v G Then a fortiori this holds for 

V € „ and so, bv l 2.61 {Au,v) ~ (/, u) for all v G which shows Au = f 

by the fundamental lemma. As the vanishing boundary values for v were only needed 
in the use of Green’s formula, not in the testing, we can infer {Au, v) = (/, v) for 
all u G by continuity. But as contains functions whose normal trace 

does not vanish, (Am, v) = Lap(M, v) can only hold if tSu = 0 , q. e. d. 

Remark. If we had used, for the weak Dirichlet problem, for the space of 

test functions instead of , then the space of Dirichlet potentials for / = 0 

would agree with But by our definition, this requires the additional assumption 

u G 

2.0 Observation. By dehnition of the spaces involved and Green’s formula [^TTl one di¬ 
rectly obtains: 

D ( 5 (hi+^d^+^)-l d (2.9a) 

_L _L 

As and are completions of spaces with = 0 and < 5 ^ = 0 , this property 

carries over: 

d(Hi’°D'=) c c (2.9b) 

2.10 Proposition (Poincare ineqnality). Let d) be a Fredholm complex where the 

inclusion map is compact. Then: 

(a) Dir is V\^- coercive on that means there is Cq > 0 with 

I'I'Ihi ^ C'q Dir(M) for all v G v _L 

(b) If the trace operator v 1—>■ t*v is a continuous mapping ^ \J‘Vl^{dM), 

then Dir is V\^-coercive on that means there is > 0 with 

klhi ^ C'q Dir(M) for all v _L t*v = 0 . 

(c) //2t is a closed affine subspace in that does not include constant forms ^ 0, 

then there is Cg > 0 with 


< Cq |Vm|l2 for all u G 2t. 

Examples are 21 = {u G u\gM = 0 } if M has a boundary, or 21 = {u G 

/m u = 0 } in cases where the integral u makes sense. By linear translation, 
the latter one gives the consequence |r — m|l 2 < Cg |Vm|l2 for functions v G 

Proof. It obviously suffices to show the last claim with |r|^i instead of |m|l 2 on the 
left-hand side. Then the proof always follows the same lines: If the claim is wrong, 
there has to be a sequence (vi) C with Irilni = 1 , and the right-hand side tends 
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to 0 . Because this sequence is bounded in there has to be a weakly convergent 

subsequence, which we again denote by (vi). This will suffice to extract a contradiction 
in all three cases. 

ad primum: Because is finite-dimensional, it is closed, and so is that means 

V _L Sj^. At the same time, Dir(t;) = limDir(r;i) = 0 , so u S Therefore, r; = 0 , 

but at the same time |u|l 2 = lim \vi\\_2 > 0 because the imbedding is 

compact. 

ad sec.: By assumption, is closed, so the argument works again, as t*v = 

lim t*Vi = 0. 

ad tertium: Here the convergence of the right-hand side means Vfi —>■ 0 strongly in 
L^, hence Vu = 0 , so u has to be constant almost everywhere, so = 0 by assumption 
on 21, q. e. d. 


Remark, (a) It is common to prove the last part constructively, see |ld.7| We are 2.11 
not aware of a constructive proof for the first and second case. 

(b) By a scaling argument, one can see that Cq = diamM with a constant Cg 
that does not depend on the size of M. 

Proposition. Situation as in \2.inb[ Let f G with f T Then there is 2.12 

exaetly one u € with Lap(u,?;) = (f,v) for all v € 

Proof. By the Lax-Milgram theorem, there is exactly one u € withLap(M,u) 

= {f,v) for all u G Now observe that not only 

but that the second summand must also have zero boundary values, so = 

0 So everything that is missing is to prove this equality also for 

V G But that is not difficult: On the one hand, Lap( • ,u) = 0 for such v, on 

the other (/, w) = 0 by assumption on /, q. e. d. 


Remark, (a) For each u* G , one also has L{u + u*) = f. So the solution is 

unique up to harmonic components. Thi s non-uniquene s s for m anifolds of higher genus 
can indeed be observed in numerics, cf. Arnold et al\ l|2ninl section 2.3.3). 


(b) If / is not orthogonal to then there is an orthogonal projection p of / to 

this space and a Dirichlet potential u for f — p. 


Hodge Decompositions 


Proposition (weak Hodge decomposition, IBruning and Leschi liuoa lemma 
2.1). There is an orthogonal decomposition 

where := HbOfi:''n . 

Proof. By definition of this sum exhausts and the last summand 

is orthogonal to the two other ones. It only remains to show the orthogonality of 
and Hi’ 0 (e:'=)-L. So let u = da with a G Then byja^w G 

hence it is perpendicular to each element in q. e. d. 


2.13 
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2.14 


Proposition ( strong Hodge decomposition, Bruning and LeschIIiqq 3 . 


__ _ _ cor- 2.5; 

SCHWARa iQQt^ . thm. 2.4.2). For a Fredholm complex, there is an orthogonal decompo¬ 
sition 


In other words: each u G can be decomposed as u = da + Sb + c with ta = 0, 

Tib = 0, dc = 0, and Sc = 0. The parts a and b can be computed as minimisers of 
F[u]{a) = {da,da) — 2{da,u) over a G and G[u\(b) = {5b,Sb) — 2{Sb,u) 

over b G respectively. 


Proof. By the Fredho lm propert y, is closed and hence convex. By the projec¬ 
tion theorem (cf. e. g. IAlt 2006L thm. 2.2), has an orthogonal decomposition 

0 gQ everything we have to show is 




The inclusion C is clear by the last part of p.pai For the other direction consider 
u G such that 

{u,dv)^Q for all i; e 

{u, 5 w)=Q forallwe 


Because of Green’s formula, this means 

{ 6 u,v) = 0 for all v G 

{du, w )=0 for all w G 


These v and w suffice to test for du = 0 and (5u = 0 in the interior of M, so u G 
For the variational property, observe that the orthogonal projection daoiu G H^d 
onto fulfills {u — da, dv) = 0 for every v G which is exactly the 

optimality condition for F. The minimiser is unique up to elements of for which 
reason we only seek a in (£^)-*-. The analogous argument applies for G, q. e. d. 

2.15 Remark. By assumption, the ranges of d and <5 are closed, so a continuity argument 
also shows that 0 (5(H°dJ7j5)]-*- is the completion of which gives 

the Hodge decomposition 


L^n’^ = d(Hd0l7fe-i) 0 ^(H°dH^+i) 0 r/ 


2.16 Proposition (Hodge Friedrichs decomposition, SCHWARZ iqq .4 thm 2.4.8). For 


a Fredholm complex that is -regular, that means the Dirichlet potential of an 

right-hand side is in there is an orthogonal decomposition 




Hidi4j0Hdii4'=nHdo<B^, 
Hidi4^0 Hdiy/n HOdvB^ 
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Together with \27T4\ this gives the deeomposition 

= d{H^+^n';-^) © n 

Proof. We only prove the first decomposition, the second one is literally the same. By 
the last statement of p.gai So if the decomposition exists, it will be 

orthogonal. We are done if we can show that every u € with u © jg ^ 

coboundary, that means there is some v G with u = Sv. 

Therefore, suppose u G and u © Then by 15.1:51 there is a Dirichlet 

potential w G for u. Let v := dw. 

(I.) It holds u — 5v E due to the assumption on u and 5v G 

(II.) As i5m = 0 and S^v = 0, we have S{u — Sv) = 0. Because the Dirichlet problem is 
Hi+i-rggular, u = {dS + Sd)w and hence d{u — Sv) = d{dSw + Sdw — 6dw) = d{dSw) = 0. 
This shows that u — Sv G And because u — Sv = dSw is a boundary, it has 

vanishing tangent component: Take any {k — l)-dimensional domain U C DM. There 
is a domain U' C M with U = U' (1 i9M, and t*{u — Sv) = d{u — Sv) = 0. For 
this reason, u — Sv G 

Now w — (5r; is at the same time in some space and its orthogonal complement, which 
can only hold if rt — (5r; = 0, q. e. d. 


Mixed Form of the Dirichlet Problem in D'" 


Observation. The D irichlet problem also has a mixed form without coderivatives 2.17 
( Arnold et a/.ll2006L sec. 7.1): If one introduces the auxiliary variable cr, which replaces 
Su in a weak sense, i. e. which fulfills (cr, r) = (dr^u) for all r G then 

Lu = f — p can be written as 

{da, v) + {du, dv) = {f — p, v) for all v G 

(ct, t) — {u,dT) = 0 for all r G 

{u, q) =0 for all q G . 

Let 6 := x x By computing the Euler-Lagrange equation, 

one sees that (cr, u,p) G & solves the equations above if and only if it is a minimiser of 

I{a, u,p) := ^{a, cr) — {da, u) — ^{du, du) + {f — p, u). 


In such a critical point, every (cr, 0,0) € © is a descent direction, and every (0, v,q) G 6 
is an ascent direction, so / has a saddle point at {a,u,p). 


Proposition ( Arnold et aZ.|[ 20 o 6 L thm. 7 . 2 ). In a Fredholm complex, the weak Dirich¬ 
let problem in mixed formulation is well-posed, that means: There is a solution (cr, u,p) G 
© for every f G and there is a constant c only depending on M such that 

\a,u,p\e < c\f\i 2 . 


2.18 
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Proof. For s = {a, u,p) and t = (r, v, q), let 

b{a, u,p] r, V, q) := (cr, r) — {u, dr) + {da, v) + {du, dv) + {p, v) — {u, q), ( 2 . 18 a) 


let i? : 6 — 7 ^ 6* be the linear operator with b{s,t) = {Bs,t)Q, g and F {f,v)\_ 2 . 
We have to show that there is a solution to the operator equation B{s) = F in ©*. 

There is an inf-sup condition for b: According to |2.i3[ decompose u = da + x + y with 
a € n and x _L Then bv l2.iobl we have IuIhco < Cg |a|L2 < Cq |m|l2 

and by i2.ioaT I^Ihuo < Cq |(i2;|L2 = Cg |(iu|L2. With t = a — C^^a, v = u + da F p 
and q = p — y, one obtains (all norms are norms here) b{a,u,p;T,v,q) = \a\^ — 
C^^{a, a) + \da^ F \da^ F \du^ F |pp -f |yp, which is (for Cq > 1 ) greater than 

^□^(I'^Ihuo + I^^Ihuo + IpIl^)- a® < oi\a,u,p\e for some a G R, we have 

shown that there is a constant 7 with 

sup > 7 |s| 6 for all s G 6 . (2.18b) 

tee die 


Once this inf-sup-condition is established, the way is well-known (|Babuska and Aziz 
1 Q 72 L thm. 5 . 2 . 1 ): Consider the linear operator B : 6 ^ belonging to b. Because 
of |i?s| 6 * > 7|s| 6; it must be injective. And as we have that for any t G © there is 
some s G © with b{s, t) ^ 0, the image of B must be the whole space ©* by the closed 
range theorem. Using the inf-sup condition once again, we get that |q;| 6 * > l\B~^a\e, 


which shows that B ^ is continuous, 


e. d. 


2.10 Remark. The only thing that cannot be inferred from our high-level point of view is 
that (H^’°U,(i) indeed forms a Fredholm complex, and that the imbedding 

is compact. The finite-dimensionality of and he nce the Fredholm property 
is proven by the inequality |i:|hi < C'(Dir(?;) -I- |i’|l2) of Gaffne^ (jiQt^ill . see e. g. 


(jiQQq . cor. 2 . 1 . 6 ). The compact embedding is Rellich’s inequality (see e. g. 
Alt! 20o(: . A 6 . 1 ), which carries over from the Euclidean case without modification. 


Both need the vector bundle structure of (H^’°f2, d), but they are obviously also true if 
(H^’°f2,(i) is a complex of finite-dimensional vector spaces. Therefore all proofs above 
literally carry over to the “discrete exterior calculus” from section A very prelimi¬ 
nary version of this attempt of formulatio n exterior calculus wi thout recurrence to the 
vector-bundle structure has been given in VON DeylenI ( 2012 II . 


2.20 Definition. When we speak of variational problems in the forthcoming sections, 
we will always refer to the Hodge decomposition, the Dirichlet problem and other 
strongly elliptic problems. 


3. Geometry of a Single Simplex 

As typical domain for the parametrisation of simplices, the numerical community 
mostly uses the n-dimensional unit simplex 

D := conv(0,ei,... ,e„) = {p G R>o : P • < 1}, 

































g. Geometry of a Single Simplex 


where is the vector in IR" with all entries 1, and is the i’th Euclidean unit vector. 

We will investigate parametrisation of simplices over D for given edge lengths and 
see how the Riemannian metric over D changes when those edge lengths are slightly 
distorted. In contrast, geometers tend to employ the standard simplex 

A := conv(eo,..., e„) = {A G R>y : A • l„+i = 1 } 

for the same purpose (here and in the following, we will use the enumeration eg,..., e„ 
for the canonical Euclidean basis of R"+^). Although the parametrisation over A is 
not a parametrisation in the strict sense, as not some R" itself is used, but some linear 
subspace of it, we will see that there will be no problems with this additional direction. 

The Unit Simplex 

Metric on TD. Consider points po,... ,pn G R" that are supposed to be vertices 3.1 
of a simplex s. As we are only interested in its isometry-invariant properties, we can 
assume that pq is the origin of R". Then the matrix P := [pi| • • • \pn] represents a 
linear map D ^ s. The first fundamental form has entries 

■“ ~ Pi ' Pj’ b J “ 1 ) • ■ • ) 

and the volume of s is computable as vols = ;^(det If £ij = \pi — Pj\ are the 

edge lengths of s, we have by the cosine law 

Cy = 5 (4 + ^oj - 4)■ (3- la) 

If only a system of prescribed “edge lengths” £ij is given, then there is a simplex with 
such edge lengths if and only if the matrix with entries + ^Oj ~ i® positive 
definite. 


Metric on T*D. Let Di, i = 1, ... ,n, he the facet (subsimplex of codimension 1) 3.2 

of D opposite to the vertex e^. The vector is normal to Di, and as normal directions 
transform with the vectors v® := P~^ei, i = 1,... ,n, are normal to the facets st 
of s. In other words, P~^ has the normals r;® as rows. At the same time, these v'‘ are 
the gradients of the barycentric coordinate functions A®, defined by the representation 
p = A®pi for any point in s. The length of these gradients decreases as the simplex’ 
height hi above pi decreases, more precisely 


u* = grad A* T Si 


(3-2a) 


this in particular implies Isil/lw*! = — vols for each i. This formula does not only hold 
for i = 1 ,..., n, but also the appropriately-scaled normal opposite to the origin po 
is the gradient of A°. Define V := • • • |u”]. As the barycentric coefficients sum up 

to one, = —Vln- By definition of V, we have V^ P = 1, which means that the 


and pj form a “biorthogonal system” (IFiedler 
gy := {v*vy^ = P ■ v\ 


20111 


], thm. 1 . 1 . 2 ). The matrix 


= 1 ,... 
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is inverse to Cy and hence represents the scalar product of the cotangent space. Clearly 
y* = PQ, that means = pj{v^,v^), so = —yi„ = Pj{v^,v^). As pj also equals 
the edge vector e^o (due to po = 0) j this shows the translation-independent form 

u* = Ckj for alH = 0 ,..., n (S-sb) 

(where i = j might be included in the summation or not, which does not matter due 
to eu = 0). In the following, we will only consider the (n -f 1) x (n -I- 1) matrix Q, 
which extends Q by a O’th row and column: 

QP := v'‘ ■ , i, j = 0,..., n. 


It is made to have vanishing row-sum and co lumn-sum, i. e. Qlr,.+^ = 0 
In the special case n = 2, we kno w from PiN KALL and PoLTHlEra 


and = 0 . 


nowadays dates the formula back to Duffim iqf;qI or even MacNeai] 


( iqq 4 Pinkall 
iQud . but as the 


formula itself is easy to discover by classical trigonometry, we value 
to computational mathematics higher than the first occurence of two 
cotangents) 

3 » 1**^1 

\ijk\ V-^ • V 


l*jl 


= cot ay. 


(as is frequently used in discrete exterior calculus, see lHlRANllboo'J) . 
straight line from jjfc’s circumcentre to the midpoint of ij, and is 
opposite to vertex k. 


the application 
opposite angles’ 

( 3 - 2 c) 

where *ij is the 
the angle in ijk 


3-4 


3.3 Definition. We say that s is [d, /i)-small if all edge lengths £y are smaller than h 
and s has volume greater than 'dh^an, where (t„ := •yrr+T/( 2 "/^n!) is the volume of 
the regular n-simplex with unit edge length, i. e. the scaled standard simplex -^A. In 
terms of the first fundamental form, detC > (d/i"n! 

Remark, (a) The standard simplex has maximal volume among all n-simplices with 
the same edge length bounds, so "d < 1 . 

(b) The parameter is l/a„ times the fulln ess 0 (s) in I Whitney! (liQt=;7l . sec. IV.14). The 
fullness parameter from Ivon Deylen et all 1I2014I I is n\ a„-d in the notation employed here. 
It would be equivalent to require a lower bound on the angles between subsimplices. 

(c) Weaker requirements on the simplex quality that still ensure well-posedness of the 
i nterp olation problem, like the famous maximum angle con dition of iBABUSKA and Aziz 


i nterp olation problem, liKe the tamo us r naximum angle con d 
(I1Q76I) . are circumstantially treated bv ISHEWCHUKl l|2002h . 


3.5 Lemma. Let a„ := n! (T„n^ for all n G N. Then the eigenvalues Xi of C fulfill 

'dha^ < \fXi < hn. 

Proof. We have to estimate ||P|| and ||P~^|| from [ 37 ^— Recall that the n-dimensional 
unit simplex has interior and boundary measure 


vol„(D) = vol„_i(5L>) = 


\/n 


(n — 1)! (n — 1)! 
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(the latter one because the (n — l)-dimensional standard simplex conv(ei,..., e„) has 
volume )■ For any n-simplex s, the radius r of the insphere is connected to 

volume and surface via vol„(s) = ^ vol„_i(9s). This can be easily seen by considering 
the simplices s* := conv(si,c), where Si is a facet of s and c is the circumcentre 
of the insphere. These s* all have volume vo1ti(s*) = ^vol„_i(si), and vol„(s) = 
vol„(si) H-h vol„(s*+i). Now, solving voln{D) = ^vo\n-i{dD) for r gives 

1 1 

r = - — > 

n + y'n 2n 

This means that any vector v G TD with length ^ < 2r can be represented as p — q 
with points p,q € A. Its image in s is Pp — Pq^ which must be shorter than the 
diameter of s. So ||P|| < nh. On the other hand, 

> Ai„inAn“x > dct C > [■dlpn\anf 


q. e. d. 

Corollary. For the norm |rc|g := w'^w^Cij on TA holds 'dha^\'w\g2 < |w|g < hn\w\e2. 3.6 
In particular, all edges are longer than Hha^. The columns u® of V form a g-ortho- 
normal basis, and \v'^\p < for all i. 

Proof. We have |r(;|g = \C^/'^w\^2, and the extremal eigenvalues of are A^jj^ and 
Amax, which shows the first claim. The second is clear from V^P = 1 and the fact that 
C~^ has eigenvalues between (hn)~‘^ and (dha^)~'^, q. e. d. 

Lemma. Assume two symmetric matrices C,C G with C being positive definite 3.7 

and \{C — C)v ■ v\ < eCv ■ v. Then also \{C — C)v ■ rc| < £\Cv ■ v\^/'^\Cw ■ 

Proof. The claim is independent of scaling v and w, so let Cv ■ v = Cw ■ w = 1. We 
will first show the claim for C-orthogonal vectors and then for linear combinations. So 
assume Cv • w = 0 for the moment. Then C{v + w) ■ {v + w) = C{v — w) ■ (v — w) = 2, 
and the parallelogram identity (polarisation formula) 

4A(n, w) = A(y + w,v + w) — A{v — w,v — w) 

for any symmetric 2-tensor A gives 4|(C — C)v ■ wj < eC{v + w) ■ {v + w) + eC{v — 
w) ■ {v — w) = ie = 4e|C'u • v\^/'^\Cw ■ Now for a linear combination, we obtain 

\{C — C){v + w) ■ up < \{C — C)v ■ v\'^ + \ {C — C)w ■ up < 2e = e|Cn • n| \C{v + w) ■ 
(u-l-in)l, q. e. d. 

Remark. The polarisation argument is also feasible for higher-order symmetric ten- 3.8 
sors, as e. g. 

6 C{v, v,w) = C(v + w,v + w,v + w) — C(v — w,v — w,v — w) — 2 C(w, w, w), 

which (together with the usual paralellogram identity) shows that all evaluations of a 
symmetric 3-tensor can be reduced to linear combinations of equal argument evalua¬ 
tions. 
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3.9 Lemma. Let C,C G be symmetric matrices, where all eigenvalues of C are 

larger than Amin > 0 (in particular, C is positive definite), and \Cij — Cij\ < eXmin/n. 
Then \{C — C)v ■ ty| < £\Cv ■ \Cw ■ 

Proof. Due to [ 3 ^ the case v = w is sufficient. Then 
\iC-C)vv\ = - C.Mv, 

* 

and the square is, by Jensen’s inequality, smaller than = n|up. As C is 

positive definite, Cv ■ v > Amin|wp, q. e. d. 

3.10 Remark. Compa re the classica l eigenvalue distortion theorem oI IBauer and FikeI (liohd l in 
the formulation of IIpse 3 (IiqqSI, theorem 2.1): If C can be orthogonally diagonalised, A is an 
eigenvalue of C, then there is an eigenvalue A of (7 with < \C~^{C — C')||. 


< 




■E 


£An 


Vi\\Vj \ = 


■(En)‘ 


The Standard Simplex 


Thing s are less obvious in barycentric coordinates, for which reason we refer to FiedlerI 


(1201 iL chapters 1 and 2) for proofs of the following statements. 


3.11 Metric on TA. We drop our assumption po = 0 and let pQ,...,pn be arbitrary 
points in R"*. Their convex hull s has a parametrisation over the standard simplex A, 
represented by the matrix P.^ = [po| ’ ’ ’ \Pn\- The Riemannian metric on TA = {v G 
Rn+i. y . _ Q| jg given by 

Ezj ■= -^\Pi-Pjf, i,j = 0 ,...,n. (3-iia) 

By 13.la) we know Cij = Eij — E^i — Eqj for i,j = 1 ,... ,n. The volume of s can be 
computed as 

vols = 4 (-detM+)"/=, where M+= ( 1? g ]f^(n+2)x(n+2) 

V—^ / 

(3-iib) 

is — i times the usual Cayley—Menger or extended Menger matrix. The volume 
element hence is G := 2 (— det 


3.12 Metric on r*A. Assume detM+ ^ 0 . The cotangent space on A consists of all 
linear combinations VidX^ with v ■ In+i = 0 . As , i,j = 0 , ... ,n, already contains 
the correct scalar products between all possible linear combinations of the dA®, it is in 
particular the correct representation for the scalar product on T*A. The ambiguity in 
the choice of gij and pb ig especially visible in the fact that normally both are inverse 
matrices, whereas for the choices made here, we only have 


Af-i 


/ 4^.2 _2g*\ 

V-29 Q )’ 


(3.12a) 


where r is the circumradius and q are the barycentric coordinates of the circumcentre 
of s. 
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Remark. IFiedle^ (l2niih uses the symbols M and Mo for the objects — 2 E and — 2 M+. 

Definition. From now on, we will only consider barycentric coordinates on the stan¬ 
dard simplex and will always use 

gtj=E,j, g^^=Q^\ i,j = 0,...,n 

and say that 5 is a (i?, h)-small metric if s is {d, /i)-small. The latter one is of course 
only possible if s is non-degenerate, in particular n < m. Note that for any c € IR, the 
matrices gij + c and g^^ + c induce the same scalar product on TA and T*A as gij 
and g^^ respectively, so we can assume that both are positive definite on with 

the same eigenvector bounds as in [g7g} 

Proposition. Although gij and g^^ are not inverse matrices, the tangent-cotangent 
isomorphism is given by the usual identities: 

{v'‘dif = gijV^dX\ {ajdX^)^ = g'^^aidj. 

Proof. By | 3 .iib| and [gn' 2 a| Ejk = + 9* for all i, fc = 0,..., n. In other words, 

QE = t + [q\---\q]. 


But if a ■ In+i = 0, then [g| • • • \q\a = 0, irrespective of the entries of q. So on TA (and 
equally on T*A), both matrices are indeed inverse to each other, q. e. d. 


Lemma. Assume real numbers x,y > 0 with \x — y\ < ex for some e < 1. Then 


- y^\ < Cpcx^, 


where Cp = 


p (1 -I- e)P ^ for p > 1, 
p (1 — for p G [0; 1]. 


In particular, 


|\/x — \/p| < 2ev^ and \x‘^ — y^\ <‘iex^ for e < ^. 


Proof. Let f ■. x ^ x^. Then \f{x) — f{y)\ < max^g[ 3 ,.y] f'{f)\x — y\. So we only have 
to find the maximum of /'. Suppose p < 1. Here /' is monotonously decreasing. If 
X < y, then max/' = f'{x) = px^~^. li y < x, then max/' = f'{y) = py^~^ < 
p(l — £)'P~^xP~^, and this case dominates the first one.—The argument for p > 1 is 
litterally the same, but with inversed roles of x and y, q. e. d. 


Proposition. Let iij be the edge lengths of a simplex, defining a (d, h)-small met¬ 
ric g on A, and let Iij be a second system of desired edge lengths with \iij — | < 

^£n~^Qf,'d^£ij where e < ^. Then there is a simplex s C IR” with edge lengths iij, and 
its Riemannian metric g over A fulfills \{g — g){v,w)\ < £|u|g |w|g. 

Proof. For the existence claim, it suffices to show that gij = —^iij is positive defi¬ 
nite on TA. By 13.151 and the assumption, we have \iij — hence 

\Eij — Eij\ < £n~^af,'d‘^h'^. Now apply [ 3 :^ to get \{g — g){v, u;)| < e|u|g|iu|g. In partic¬ 
ular, g is positive definite for e < I, q. e. d. 
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3.17 Remark. Although the proof of [3^ seems to be fairly rough, the fullness parameter d 
is essential to |3.i6| and cannot be weakened to a bound e. g. on the minimal edge length 
(which would be equivalent to a minimal angle bound , which is a sharp tool for some 
interpolation estimates, see Babuska and AzizjliQ 72 h In fact, edge lengths 2, 1 + e 
and 1 + e of a thin triangle may only be relatively distorted by a factor S G ] 1 + £[ 

to guarantee the existence of a corresponding simplex. 


3.18 Corollary. Situation as Then g and g are equivalent norms, (1 — e)g{v,v) < 

g{v,v) < {l+e)g{v,v). 

3.19 Proposition. Situation from \g. 1 6‘[ Then the estimate carries over to all higher tensor 
spaces TJA with constants depending on r and s. 


Proof. Let us first consider the extension of g to bivectors, that means elements of 
T2 A. Without regarding the tensor-product structure, we could just say this it is a 
scalar product with components gijgki, where (i, k) and (j, (.) are the indices of the first 
and second argument respectively. And these components are almost equal to those of 
g (we abbreviate \gij — gij\ < Sgij): 


\ 9 ij 9 kl 9 ij 9 ki\ "S: \ 9 ij 9 ki 9 ij 9 ki\ + \ 9 ij 9 kl 9 ij 9 kl\ 

— ^ 1 9 ij 9 kt \ “1“ 1 9 ij 9 ki I 

— T d{\ -\~ 6 ^\ 9 ij 9 k^\ Si ‘^^\ 9 'ij 9 kE\• 


For the x matrices 9ij9ki and 9ij9ki, applyagain, q. e. d. 

3.20 Corollary. Situation as in \g.i6] The volume elements of 9 and g fulfill |G — G| < ceG 
with a constant c depending only on n. 

Proof. The usual proof would be to use G = 2 (—det together with the first 

order approximation o f the determinant : det(F -|- tA) = (l -|- t tr F~'^A + 0{t^)) det F 
(a classical reference is BELLMANN||m6oi pp. ghsqq., although the nicest proof that we 
know is in Eschenburg and .foSTlbnoTl lemma 8.2.1). However, we want to control 
more than the first order. So observe that if di is a coordinate-induced basis of the 
tangent space, then G = |cfi A • • • A ( 3 „|g, and g — g on the space of u-forms is controlled 
by the previous proposition, q. e. d. 


3.21 Proposition. Let g,g be two metrics with \{g — g){v,w)\ < e\v\g 
Then the metric g'’^ on the cotangent space fulfills \g{a,P) — g{a,fi 
a,l3 GT*A. 


w 


where e < 


g> 

< 2e \a 


1 

2 • 

for 


Proof. Again, we do not use the conventional approach to bound Mj_ by the differen - 


tial of the matrix inversion ^(A) = —F 


bonj . lemma 2.8). Instead, the definition of 


i:(v) N 

'a - ^“RT 

lain. On the other hand, if w realises 


For every u S TA, we have |a|| > (x^)^ > 




1 

H-e 

then 


~^AF~^ (IDeuflhard and HohmannI 

g on T*A as operator norms will help: 

^ in particular for the v realising 

2 _ / a(w) ^2 ^ 1 f a{w) \2 

a — i-elTunrl ■ 


1 h|„ > 


1 

1 -I- £ 


\a\ 


< 


< 


1 

I — £ 


|a| 


2 
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Now as e < i by assumption, we have < l + 2e and > 1 — e > 1 — 2e, therefore 
||a|g — |q;|?| < 2e|Q!|g, which suffices due to [ 3 ? 7 l q. e. d. 

4 . Simplicial Complexes and Discrete Riemannian Metrics 

In computational geometry, it is common to describe simplicial complexes as the union 
of simplices in Euclidean spaces with appropriate conditions on their intersections. We 
consider these intersection conditions as tedious and use the more abstract definition 
via barycentric c oordinates, as is usually done in topology. We follow the lines of 
Munkr^ (|iq8ziI1 . but we repeat the definitions in order to directly deal with abstract 
simplicial complexes as (almost everywhere smooth) Riemannian manifolds. 


Non-Oriented Complexes 

Definition. An n-dimensional combinatorial simplex (n-simplex) is a set of n+1 el- 4.1 
ements, its ^-dimensional subsimplices are subsets of cardinality £-1-1. An n-dimensio¬ 
nal combinatorial simplicial complex is a collection ^ = (.S°,..., .ft"), where each 
.ft^ is a collection of Adimensional simplices such that if i is a fc-dimensional subsimplex 
of s € ft^, then t € .ft^. The complex is called regular if each simplex is contained in 
an n-simplex and each (n — l)-simplex is the subsimplex of at most two n-simplices. 

When we speak of simplicial complexes, we always mean regular ones. 

An (n — l)-simplex f is called a boundary simplex if there is only one e G ft" with 
fee. The (n — l)-dimensional complex formed out of these boundary simplices and 
their subsimplices is called the boundary complex 5ft of ft. 


Notation. We use special notations for the most interesting dimensions (here k is any 
dimension between 0 and 1, kept fixed inside an argumentation): 


vertices or * G ft° t G ft*^ ^ 

edges ij G ft^ s G ft^ 


facets f G ft"-^ 
elements e G ft" 


Sometimes we will also use the convention t G ft^ and s G ft*^’*'^. In every case 1 will be 
one dimension smaller than s. 


Definition I TOM Dieck ( 2 nnnll . p. 63). Let s := {po, ■ ■ ■ ,Pk} be a combinatorial 4.2 
fc-simplex. For a function A : s —R, abbreviate X{pi) as A*. The geometric realisa¬ 
tion of s is rs := {A : s —>■ [0; 1]: A • l„+i = !}• For a complex ft = (ft°,..., ft"), the 
realisation is defined as rft := re. 


Remark, (a) This definition is equivalent to, but much more elegant than the usual 
way of “annotating” the vertices of the euclidean standard simplex A with the elements 
of ft*^, considering the disjoint union of |ft"| many such annotated simplices and glueing 
them whenever two sides have equal annotations. 

(b) By setting A* = 0 for all unused vertices pi G ft°, the elements of rft can naturally 
be considered as functions A : ft° —>■ [0; 1] with supp A = s for some s G ft^. 
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(c) We say that some property is fulfilled piecewise on rA if it is fulfilled on each 
re, e € 

4.3 Proposition. Let A. he an n-dimensional simplicial complex (regular, as always). 
Then rA is an n-dimensional manifold, which is smooth everywhere except at (n — 
2)-simplices. 

Proof. For each f G belonging to e, e' G A", define a chart Xf : re U re' —^ R" 

in the following way: Without loss of generality, assume e = {po,... ,pn} and e' = 
{pi,... ,pn+i}- Let ei,..., e„ be the usual euclidean basis vectors in R", let eg be the 
origin and e„+i := ^(ei + • • • + e„). Then the convex hulls D := conv(eo,..., e„) and 
D' := conv(ei,..., e„+i) are isometric (up to a change of orientation). Now define 


Xf(A) := X'ei = 


A°eo H-h A”e„ on re, 

A^ei H-h A”+^eji+i on re'. 


(4-3a) 


This Xf is a bijection re ^ D and re' —>■ D'. For any other chart xy that also covers re, 
the chart transition is an affine map that maps E to either e or e', hence rA is smooth 
in the interior of each n-simplex. 

Around an (n — 2)-simplex, we do not give a chart, but we just remark that a 
topological manifold is sufficiently defined by a finite cover of closed chart domains. 
Open charts are only needed for the definition of smooth functions. 

Note that another choice than eg,..., Cn+i would have led to the same differentiable 
structure on rA. (as long as the convex hulls are full-dimensional simplices in R”), 

q. e. d. 


4.4 


4-5 


Observation, (a) Consider 5 € Af. By rs is a full-dimensional subset of the 
/c-dimensional affine space {A : s —>■ R: A • !„+! = 1}, so its tangent space is = 
{u : s —>■ R: u • l„+i = 0} at every internal A G rs. As rs is an affine space, we 
deliberately drop the foot point A in most cases, just as we do with TA. 

(b) The obvious linear isomorphism A —>■ rs is 1 —^ rpi. This means that A G A and 

V £ TA are mapped to X'rpi and v'rpi. 

(c) The realisation of the boundary complex is the boundary of the realisation: dr A = 
rdA. In particular, rA. is a manifold without boundary iff each (n — l)-simplex in .h 
belongs to two n-simplices. 


Definition ( Wardetzky 2006! or, similar but shorter, Hildebrandt et al. 


2nnt 


Define a differentiable structure on rA by the requirement that some function is 
smooth (or of class C^’“) if it has this smoothness property piecewise and is continuous 
up to the boundary. Consequently, define H^' as the completion of C^' with respect to 
the norm. 


4.6 Definition ( Bobenko et a?.ll 20 ioll . Let .h be a simplicial complex. A function t •. A} ^ 
R>o with the property that Cij from [37^ is positive semidefinite for each e G .S" is 
called a discrete Riemannian metric. In particular, (. fulfills the triangle inequality 
on each t G . 
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On each T\r5, s £ the discrete Riemannian metric i induces a usual Riemannian 
metric gi{v,w) = v^w^gij by gij := cf. |3.iia| As this metric does not change 

with A, rs will be flat. When we deal with a piecewise flat metric, we always assume 
that it is defined via a discrete Riemannian metric. 


Observation, (a) Let t be a facet of s. The restriction of £ to edges in t is a discrete 
Riemannian metric for itself, and its induced Riemannian metric gi^t on rt is the restric¬ 
tion of g(. So the glueing of two supersimplices s,s' of t along t is done isometrically 
with respect to gi. 


(b) Consequently, every set [/ C rA that does not contain any (n — 2 )-simplex is 
flat. In fact, also (n — 2 )-simplices might be included if they have some flat neighbour¬ 
hood, which is eq uivale nt to requiring that their internal dihedral angles as defined by 
p. 412) sum up to 1 . In this sense, curvature of piecewise flat 


Cheeger et al 


spaces is concentrated in the (n — 2)-simplices. 


(c) If Pi are points in Euclidean space with iij = \pi — pj|^2, then gi coincides with 
the pull-back metric gs of s := conv(po, • ■ ■ ,Pk) to the standard simplex A. Hence, 
: {rs,ge) —>■ {A,gs), X !->■ A'ci and x® : {rs,gi) —>■ {s,P), A >->■ A®pi are both 
isometries. 


4-7 


(d) In the constructionone may use points qi with distances \qi — qj\i2 = iij 

instead of the points Ci (of course, j^o ~ Qn+i\p does not undergo any restriction). Up 
to Euclidean isometries, these qi are unique. This defines an atlas {xf. f £ of 

isometries. 

(e) Consider a triangle {pi,Pj,Pk} £ shortly written as ijk. By the usual trigono¬ 
metric formulas incorporating only edge lengths, one can define angles opposite to 
the edge ij and area \ijk\ on the basis of £ only, without using gi. The metric on 
r.fl obtained by the requirement 


|jfcA| = X'lijk], \ki\\ = X^\ijk\, \ijX\ = X'^\ijk\. 


is the same as the metric induced by gi. The generalisation of this approach to higher 
dimensions is of course feasible and natural, but notationally tedious. 

Proposition. Let r^ he a realised simplicial complex with a piecewise flat metric g. 4.8 
Consider two adjacent elements e = {po ,... ,p„} o-nd t' = {pi,... ,p„_|_i} £ .fl" with 
common subsimplex f. Then for any A in the interior of rf, the differential of the 
transition map Trt ^ Trc' has dual 


idT,,^,f : TVe ^ TVe', 


dA* !->■ dA* for i = 1,... ,n, 

dX°^-j^dX-+\ 


(4.8a) 


Proof. It is clear that the common differentials dA^,..., dA" remain unchanged. Under 
an isometric embedding as in |4.7d[ (dA°)® and (dA"+^)® are normal to the common 
facet (cf. pointing in opposite directions, which gives f^^dA® = — dA"~*~^, 

q. e. d. 
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Remark, (a) We have chosen to give and not just to obtain a nicer 

formula. One could as well say 

grad A° |j^nji| gradA”+^ 

(b) Formally, is only defined on re U re', where it is the identity. But the charts 
and from [477^ can be extended to some neighbourhood of the standard simplex, 
as (rz^gi) and {rz',gi) are glued isometrically. 

Oriented Complexes 

4.9 Definition („Nanu, Sie kennen .hor noch nicht?“). Let F be a set. Define an equiv¬ 
alence relation ~ on the set F”+i of (n -|- l)-tuples over F by a ~ 5 iff there is a 
permutation with positive sign that maps a into b. Let [uq ,..., On] be the equiva¬ 
lence class of a G F”. The quotient of under ^ is called the set of oriented 

fc-simplices with vertices in F and is denoted by [F”]. 

If [6] G [F^] is an oriented simplex, its facets are the oriented {k — l)-simplices 
obtained by dropping one of its vertices: [bo,... ,bi,... ,bk] < [t>o,..., 6fc]. The sub- 
simplices of [6] are obtained by dropping one or more vertices. If dimensions do not 
matter, we also abbreviate [a] < • • • < [6] as [a] < [6] if [a] is a subsimplex of [6]. 

An oriented simplicial complex .hor of dimension n with vertex set F is a collection 
..., where C [F^], such that [a] < [b] for some [6] G implies [a] G 
The complex is regular if no vertex occurs twice in any simplex, each simplex is 
contained in an n-dimensional simplex, each (n— l)-simplex is the boundary of exactly 
one n-simplices, and each two n-simplices in have different vertex sets. 

If Ao-c is a regular orientable simplicial complex, we denote the corresponding complex 
made out of non-oriented simplices by A. The realisation of an oriented complex Aa 
is defined as rAo^ ■= rA. 

4.10 Remark, (a) There are exactly two distinct oriented simplices with the same set 
of vertices oq, ..., a„, which we write [oq, ..., an] and [oq, ..., an]~■ As non-oriented 
simplices were defined as sets, each non-oriented simplex corresponds to two oriented 
simplices. So the last condition on a regular complex says that not [a] and [a]“ G [F"] 
can belong to an n-dimensional complex at the same time. 

(b) The vertices of a non-degenerate euclidean simplex s = conv(po, ■ • ■ ,Pn) C R™ 
can be ordered such that P = [pi — po \ ■ ■ ■ \pn — Po] as in [^TT] has positive determinant. 
This is what we call the canonical orientation of {po,... ,Pn}- (On the other hand, if 
Po,... ,Pn are not taken out of some oriented space, there is no canonical choice.) 

4.11 Proposition. rAor is an orientable piecewise smooth manifold for any regular oriented 
simplicial complex .Sor • 

Proof. We will show that if we use only those charts from the proof of [4:3] that respect 
the orientation of n-simplices, we obtain an oriented atlas of rAor- 

Suppose there are two simplices 5 , 6 ^ G that share n vertices, say pi ,... ,pn. Then 
because t := [pi,... ,p„] can only be contained in one of them, we can assume that 

s = [P 0 ,Pl, ■ ■ -^Pn], s' = [pn+l,Pl, . . . ,Pn]" 
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for two vertices poIPra+i £ -^or- Now let Xi be the chart as in |4.3a[ Obviously, [eo,..., e„] 
and [cn+i, ei,..., e„]“ are both canonically oriented. As there was no choice in this 
construction, every other chart that covers s must also map rs to a euclidean simplex 
with this orientation, therefore every transition map is orientation-preserving, q. e. d. 

Barycentric Subdivision 

Definition. Let be a simplicial complex, regular as usual, and := U • • • U A" 
be the set of all its simplices. An (ascending) fe-flag in is a set a := {uo,. .., afc} C 
(^*)fc+i that, if its elements are ordered by magnitude, C Uj+i. In other 

words, a fc-flag is a sequence of fc -I- 1 nested simplices. If £ .8"% we also write 
a = {{no), ..., {rik)), meaning that (j) is a “generic j-simplex”. 

Remark, (a) The notation (z) is uncommon, but not more ambigous than other 
authors’ notations such as ctL Our notation is made to save double subscripts. 

(b) Of course, flags are simplices, only in some special complex. But having a different 
name will (hopefully) prevent confusion. The term “flag” is borrowed from algebra, 
where it signifies sequences of nested linear spaces, whereas set theory mostly speaks 
of “ascending chains” for nested sets. But the term “chain” already has a canonical 
meaning in simplicial homology theory, and in sectionj^we need to use both at a time. 

(c) All elements of a fc-flag lie in a common n-simplex. An n-flag contains exactly one 
fc-simplex for each k. 

Example. Suppose A consists of one triangle ijk, its edges and its vertices. Then the 
0-flags are the elements of .S* (to be totally precise, the 0-flags are singletons containing 
elements of .ft*). The 1 -flags are of the form (( 0 ), ( 1 )), that means combinations of a 
vertex and an edge containing it, or of the form ((0), (2)), i. e. a vertex and the triangle, 
or ((1), (2)), an edge and the triangle: 

{i, ij}, {i, ik}, {i, ijk}, {ij, ijk}, {ik, ijk} and similar for the vertices j and k. 

The 1 -flags (( 0 ), ( 1 )) are interpreted as straight line segments from the point r( 0 ) to a 
point somewhere on the edge (1), and the flags ((1), (2)) connect the points 
to the “barycentre” A^2) • The 2 -flags consist of a vertex, an edge containing this vertex, 
and the triangle, they are all of the form ((0), (1), (2)): 

{i, ij, ijk}, {i,jk, ijk} and similar for other vertices. 

Definition. The (barycentric) subdivision sdft of the complex A is a complex of 
the same dimension whose /c-simplices are the fc-flags in K. 

Suppose there is some A^ £ rs given for each s £ ft*. Because of |4.13e[ the mapping 

r(sdft)° —?► rft, r{s} As 

can be uniquely extended to a continuous, piecewise affine mapping i : r(sdft) —>■ rft, 
mapping the realisation of a flag r(ao,..., a^) to the convex hull of Auq ,..., An,. If £ is a 


4.12 


413 


4.14 


415 
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discrete Riemannian metric on M., then r(sd A.) can be endowed with the induced metric 
•^{5},{j3'} = |As — Xs'lgn i becomes an isometry. Let r' := ior he the “realisation of 
sd M. in rA!\ 

If .^r is an oriented complex, one can obviously define an oriented subdivision by 
considering the n-flags as tuples instead of sets and using the orientation induced by r'. 

4.16 Observation. There are several obvious conclusions from the fact that r'{oo, ..., at} 
= conv(Aao,..., Aaj,). Most prominently, one can decompose the realisation of a fc-sim- 
plex ri into the realisations of fc-flags ending at t. The boundary 9 (rs) of a realised 
[k + l)-simplex s is covered by (the realisation of) fc-flags ending at facets of s: 

rt= (J r'((0),...,(fc)), 9 (rs) = |J r'((0),..., (fc)). 

(A:)—t {k)(Z 5 

Definition. For s G , aggregate the n-flags containing s in the neighbonrhood 
t/(s) of s and the (n — fc)-flags starting with s in the dnal cell *s: 

C/(s):= IJ r'((0),...,(n)), r(*s) := |J r'((fc), ..., (n)). (4.16a) 

{k)—5 (A:)—s 

4.17 Observation, (a) The flags occuring in t/(s) must obviously be different from the 
flags occuring in f 7 (s') for s 7^ s', so these neighbourhoods form a covering of rA. with 
disjoint interior for each fc. 

(b) The set of all n-flags “running through s” can be represented as a product of two 
flag sets: fc-flags ending at s, whose union is rs, and the (n — fc)-flags beginning at 
s, whose union is r(*s). For us, the latter is just a way to write this uni on, we will 
not de fine the combinatorial dual of .S. The interested reader is referred to IMunkresI 
()iq8ziL § 64). 

(c) The boundary of a neighbourhood consists of those flags where (fc) = t is left out: 

dU{s)= y r'(( 0 ),...,(fc),...,(n)). 

{k)=s 


This can be seen as follows: The boundary of any n-flag a consists of (n — l)-flags a! 
where any one of the elements in a is left out. The boundary of the union U (s) now 
consists of those facets r'a' where some element is left out and there is no other n-flag 
belonging to U{s) on the other side of r'a'. This second condition is satisfied only if s 
is left out, because if (i) 7^ s is left out, there is another flag ((0)',..., (n)') running 
through s with (i)' 7^ (i). 

4.18 Lemma. Let rA he a realised simplicial complex with piecewise flat metric, and let 
As be the circumcentre of rs for each 5 G A*. Then for each n-flag (( 0 ),..., (n)), the 
vectors := A^^+i) — A^^) are perpendicular to r{i) and thus pairwise orthogonal. 

Proof. Consider the two-dimensional case: If Xijk is the circumcentre of r(ijfc), then 
\'Vi,ijk\ = The “circumcentre” of the edge ij is Xij = ^{ri rj). So we have 
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two equilateral triangles (Ai,Ay , A^fe) and (Aj, Ajj, Ayfc), which must hence have the 
same angle 7r/2 at Ay . The same argument applies in higher dimensions: If At is 
the barycentre of ri, then all u^t have the same length. If t is a facet of s, then all 
triangles (i, At, A^) are equilateral and hence have the same angle at At. This can only 
be (because the vectors Vi^t span the supporting plane of rt) if ut,s is perpendicular to 
the supporting plane of t, q. e. d. 

Corollary. If the complex is well-centred, i. e. if the circumcentre A^ always lies 
inside rs, then the volume of a k-flag a € (sd.S)^' can be computed as 

k ^1 = l^ao.Qi I ■ ■ ■ l^^afc-i.Ofc I = 1^00,01 A • • • A \. 

Together with \4.ij^ we get for an n-dimensional complex 

|s| l*s| = |t7(s)| foriGM.’^, 


where we have written |s| instead o/|rs| for short, as we will always do in the following 
(no ambiguity will occur, as the magnitude fc + 1 of 5 is always indicated by saying 
5 e 

Remark. This last volume equation is, to the best of our knowledge , not yet used 


i n dis crete exterior calculus, but frequently in Regge calculus ^ see e. g. IMiller et al. 
( 201‘i ), and its use for discrete calculus seems to date back to MillerI ( iqqtII . 
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5. The Karcher Simplex: Definition 

Notation. Let d be the geodesic distance on M. For points pq, ... ,pn G M and the 
n-dimensional standard simplex A, consider the function 

i? : M X A —>■ R, (a, A) I—>■ A° d^(a,po) + • • • + A" d^(a,p„). 


Convexity. Let cvrM be the largest radius such that all geodesic balls BcvrM(p), 5-i 
p G M, are convex in the sense offTT^ It can be estimated by 


cvr M <\ min {, inj M }, 
where inj M is the injectivity radius of the manifold ()Cheeger and Ebin 


5-14), 

B 


and d(p, •) is convex in B, 
m{p), all functions d{q, ■), q G 


107.4 thin. 
Consequently, for a smaller ball B := 


cvr m{p)- 

M, are convex in B. 


Remark. One knows that inj M > min{mininial distance between conjugate points}, {|min. 
length of a closed geodesic}, and because conjugate points need to have distance greater 
than ■kI\/Co by Rauch’s comparison theorem (which is true also if the sectional curva- 
ture is somewhere negati ve, so the restriction of nowhere positive sectional curvature in 
ICheeger and EbinIIiq';^ . corr. 5.7, is only historically determined and factually unnecces- 
sary), 


cvrM < i|—imin. length of a closed geodesicj. 

The probably best known estimate for the latter term in arbitrary dimension is 


L(7) > 2?! 




volM 


vol S’" ^ sinh( yUo diam M) 


for a closed geodesic 7 iIHeintze and K archer] |ic[ 78|L where S"^ is the m-dimensional unit 
sphere in For even dimensions, KlingenbergI (Iiqc;(J 1 showed L{'y) > 2 'k\/Cq for 

orientable and 1/(7) > 'K'/Cq for non-orientable M. 


Observation. Local minimisers of E{ - ,X) for fixed A are zeroes of the section F : 5.2 

M X A ^ TM, 


F(a, A) := XXi\a with Xi = i gradd^( • ,pi) from FTT^ 


If the points pi lie in a common ball B := B1 cvrAiip) some p G M, then E{ - ,X) 
is convex, and hence there is a unique minimiser in B. But this can be sharpened: 

Proposition. Let pq, ... ,Pn be contained in a ball B = Br((?) C M with r < ^ cvr M. 5.3 

Then for each X in the standard simplex, E{- , X) has a global minimiser in B. 
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Proof. Let x be the minimiser of E{-,X) over a G B. The key o bservation i s that 
this function cannot have more than one minimiser in ]B2r(AL), 


as 


Kendall 


thm. 7.3) has shown. His approach uses convex functions on M; but IAfsarii (I2011I) 
has given a direct proof: K( •, A) points inwards on the boundary dlB2rip), and the 
Hessian of E is positive definite at all critical points. So by the Poincare-Hopf index 
theorem, E needs to have exactly one zero. 

Outside lB2r(p), there cannot be any minimisers: If a £ M \ lB2r(9), then d{a,pi) > 
r > d{x,pi) and hence E[a, A) > E{x, A), q. e. d. 


Assumption. From now on, we only consider pq, ... ,pn that lie in a common ball B 
of radius smaller than ^ cvr M, in particular Cor < 

5.4 Definition. For a given A £ A, let a;(A) be the minimizer of E{ - ,X) in B. We call 
this map a; : A — >■ M the barycentric mapping with respect to vertices pi, and its 
image s := x{A) the corresponding Karcher simplex. 


5-5 


Remark, (a) In case M is the Euclidean space, x is just the canonical parametrisa- 
tion A I— >■ A*pi, because d^(p,a) = |a — p|^2 gives Xi\a = a — pi. 

(b) For A* = 0 , the value x{X) is independent of So the subsimplices of the standard 
simplex are mapped to “Karcher subsimplices” which only depend on the vertices of 
the subsimplex. 


(c) If Ci is the i-th Euclidean basis vector of ]R"+^, then xifCj + (1 — t)ei) = 7(t), 
where 7 is the unique shortest geodesic with 7(0) = pi and 7(1) = Pj. 

(d) Because x is continuous, the Karcher subsimplices form the boundary of a Karcher 
simplex: d{x{A)) = x{dA). 


(e) Concerning the definition of x, we will not make use of the fact that all A® are 
positive beside in I6.16I It was only nee ded to have an easy access to the well-defined¬ 
ness of the minimiser. Sander (|2ni‘jl showed that negative weights also lead to a 
well-defined minimum if the pi are contained in a ball whose radius is bounded by a 
constant depending on inj M, the curvature of M and max |A® — A-^ j. 


5.6 Proposition. If all pi lie in a totally geodesic submanifold S, then x{A) C S. There¬ 
fore the usual no tion of simplices i n spaces of constant curvature as convex hull of the 
vertices (cf. e. g. TiLURSTONlliQQTl . ex. 3.3.6) is recovered. 


Proof. Let ki := ^ d^( • ,pi). If B is convex, then so is i? fl S. So E{X, ■ ) is convex on 
H n S', and hence there is a unique minimizer a ol E{X, ■)\s, so there are coefficients 
A® with A® grad(fci|5) = 0 at a. As S is totally geodesic in M, it holds grad(fci|s) = 
(gradfci)ls = A'ils. Hence A®Ai = 0 at the point a £ S, q. e. d. 


5.7 Proposition. Define bundle maps a, Ay (for v £ by 

a\x. TxA^T^(^x)M, V ^ - v"Xi\^,^x), 

Av\x ■ T^(x)M ^ T^(x)M, V ^v'"VvXi\x(^x)- 
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5- The Karcher Simplex: Definition 


If X is smooth at X G A, then its first and seeond derivative there fulfill 

Axdxv - a{v) = 0 , (5.7a) 

A\Vdx{v, w) + Ayjdxv + Aydxw + ^^'^dxw,dxv^i = 0- (5-7b) 

Proof. Similar to the proof of the implicit function theorem: The derivatives of f in a 
direction {V,v) G T^f^x)AI x T\A are 

'^{v,v)F = + v'^Xi = A\ V — ax{v). 

Now consider a curve 7 : 11—>■ X+tv with derivative j{t) = v. We have F{x{'-f{t)),'y{t)) = 
0 , and differentiating this gives, like in the usual proof of the implicit function theorem, 
0 = DtF = V(dxv,v)F■ This shows the first claim. 

The second claim is a totally analogous computation, but involves second covariant 
derivatives Vyw VyVw — VvvW- We differentiate F once more and obtain 


{w,w)'^ {v,v)F = + A'VrvVyXi 

= + X^V^yXi + A*VvwV-^i 

Again, F neither changes in direction {dxv,v) nor in direction (dxw,w), so we get 

0 = uXX dxvXi + dxwXi + ^’^^dxw.dxv^i + V dx ^dx yXi 

AxXdxwdxv. 

And because A is flat, we have Xdxwdxv = Xdx{v,w) by fTTTl q. e. d. 


= Ay,dXV + Aydxw + yXl^y,,l^yX, 


Remark. One can consider A as a point meas ure on M that ass igns the mass A* to the point 
Pi. For a general probability measure p on M. IKarcherI (IiqvtII speaks of the minimiser of 


A, 


fa) := f dfi{a,p) Ap(p) 
J M 


as Riemannian centre of mass, but the subsequent l iterature has m ostly called it the 
K archer mean w ith respect to the measure p (cf. e. g. IJermynI bood . probably initiated 
by Kenda^ monl'l . The concept seems to go back to Cartan (se e the historic overview 


^JSA^ boii 1 . but has not been used by others until the work of I Grove and Karchei^ 


5.8 


lliQvdl . 

Karcher himself used the centre of mass to retrace the standard mollification procedure 
of Gauss kernel convolution in the case of functions that map into a manifold. Considering 
the centre of mass as a function from an interesing fini te-dimensional space of measures 
into M, as we use it, has been d one by Rustamov bpinl'I. and the baryc e ntric coordinates 
we deal with have been used by ISand^I ll2nid . bnidl and IGrohs et al\ ll2niqh . All these 
emphasise the possibility to glue the Karcher simplices along corresponding facets, but do 
not investigate the distortion properties of this mapping. Recently, we were informed that 
Dyer and Wintraecken (Rijksuniversiteit Groningen) have also proven a result similar to 
|6.i7| by Topogonov’s angle comparison theorems. However, this approach seems not to deliver 
an analogue of l6.22l 

In contrast, there i s a large literature for “ barycentric coordina tes” on general convex poly¬ 
gons in the plane, cf. IWarren et al\ ll 2007 l ): [Meyer et al\ (l2002ll and references therein. 
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B. Main Constructions 


6. Approximation of the Geometry 


Estimates for Jacobi Fields 


4 ■ 


(6.1) 


It is well-known that locally, Jacobi fields grow approximately linearly: 

\J{t) - P‘’0(J(0) + ti(0))| < Cot'lil'd J(0)| + ifl j(O)l) for Cot^ < 

In fact. IJo^ ( 20iiL thm. 5.5.3) proves that the left-hand side is smaller than |J(0)| 
(coshct—1)-|--^| J| (0)(i sinhct—t) for c = y^Co. By Taylor expansion and -^\J\ < | j|, 
this estimate is weakened to our form. Clearly, if the values J(0) a nd J(t) are given , 
one can expect J to behave like — J(0)), but as Richard IDedekindI ( i8q4 

p. 11) said, “nothing that is provable ought to be believed without proof in science.” 


6.2 Situation. Suppose 7 : [0; r] —>■ M is an arclength-parametrised geodesic with 7(0) = 
p and 7 (t) = q, and V G TqM. Let s 6{s) be a geodesic with (5(0 ) = 7(t) and 
i5(0) = V. Define a variation of geodesics by 


c{s,t) := expp (y(expp) M(s)). 

Then T := i9jc is an autoparallel vector field and J := dsC a, Jacobi field along 1 1—>• c(s, t) 
for every s with boundary values J(s, 0) = 0 and J(s, r) = (5(s). 

6.3 Proposition. Situation as m l 6.2! Define V{s,t) := P^’'^S{s) and £{s) := r|T|(s), the 
distance from p to S{s). (By construction, |R(s,t)| is constant in s and t, and |r(s,t)| 
is constant in t, so we drop the unneeded arguments.) //C'o^^(s) < ^ for all s, then 

I J(s,t) - ^V{s,t)\ < 2Cof{s)\V\, 

\jis,t) - ^Vis,t)\ < ICoe{s)\T\is)\V\, 

|J(s,t)| <Co\T\\s)\V\. 


If the derivatives of R up to order k are bounded by constants, then so are the t-derivati- 
ves of J up to order fc J- 2. 


Proof. From the usual Jacobi field estimates, e. g. I.IostI (2011, thm. 5.5.1), we get that 
|J| is increasing for all t < r in case Cq£^ < By the Jacobi equation |1.15it[ this 
already shows the last claim. Now observe J{s, 0) — 0 j(s, 0) = 0 and 


Dt [j{s,t) - t j(s,t)) 


t\J{s,t)\<Cot\T\^is)\V\. 


So the vector field U : t J{s,t) — tJ{s,t) vanishes at t = 0, and we have bounded 
its derivative. The fundamental theorem of calculus ll.lpai gives 

\j{s,t) -tJ{t,s)\ < \Cot^\T\'^{s)\V\. (6.3a) 

By J{s,t) = V{s,t), we have 

|R(s,r)-rj(s,r)| < \Co£\s)\V\. 
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6. Approximation of the Geometry 


Now J{s, t) — J{s, t)\ < (t — t) max | J| by the mean value theorem, and thus 

\V{s,t) - r j(s,t)| < |P*’'^F(s,r) - j{s, r)| + t\P*''^ j{s,T) - j(s,t)| 

< \Co\V\e{s) +Co(r-t)r|h^||T|2(s) 

< lCoe{s)\V\. 

This proves the comparison between J and . For the comparison to J, consider 

|d(s,t) - ^V{s,t)\ < \J{s,t) -tj{a,t)\ +t\j{s,t) - ^V\ 

< ^Cot^\V\ \T\\s) + ICotT\V\ |r|2(s) 

< 2C^e{s)\V\. 

The statement about higher derivatives of J is justified by the fact that one can easily 
give linear ode’s for them by differentiating the Jacobi equation, e. g. J + i?( J, T)T + 
R{Jj T)T = 0 as T = 0, q. e. d. 


Remark. These estimates are scale-aware with respect to reparametrisations of 7 : If 6.4 
t is replaced by At, then also r becomes At, whereas \T\ becomes y|T|. So ^ and hence 
the whole first inequality in | 6 . 3 | is scale-independend. As J = VtT (loosely speaking), 
the second inequality scales with 1/A and the third one with 1/A^. 

Lemma. Consider some function U : [0; r] —>■ R"* satisfying the linear second-order 6.5 
differential equation U = AU -\-B with smooth time-dependent data A(t) G and 

B(t) G R™ as well as boundary conditions U{0) = U{t) = 0. Then, provided that 
1 ^( 011 '’’^ ^ 1 everywhere, it holds 

|f7(t)| < 3|i?|T, |t/(t)| <6|i?|t(T-t). 


Proof (by David Glickenstein). Denote the maxima of ||A|| and |i?| over [0;t] as a and 
b respectively. As t/ is C^, there is an upper bound K for \U\ on [0;t], attained at 
t = d. As this point is critical for |17p, we have {U, U) = 0 there. So 


+'d‘^\U{'d)\^ = \U-tUf(d) 


■d 


\j tUdt'^ 
0 


< 


I f t{aK -\- b) 


{\d'^{aK + h)f. 


This shows K < ^'d^(aK -|- b), so K < r'^b by assumtion and hence |D| < 26. (Note 
that this argument, which first roughly bounds |?7| and then re-inserts this bound 
into the differential inequality to get a sharper estimate, is the same as in iii.iyj iq.) 
Furthermore, the inequality chain also shows 'i?|f/(i?)| < 6 d^, which means |t/(d)| < 6t. 
For other values of t, we have \U{t)\ < br -\- / |t/(t)| < 36t and, by integrating once 
more, \U{t)\ < Sbrt as well as \U{t)\ < 3bT{T — t), whose minimum is dominated by 
66t(T — t), q. e. d. 
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6.6 Proposition. Sitation as in 1 6. gl Co^^(s) < Then 


\DsJ{s,t)\ < 90Co,i(s)^ |Pp \T\{s) 

<90CoAs)t\V\^\T\{s), 


\DJis,t)\ <50C'o,i(s)|Pnr|(s). 


with ( 70 , 1 ( 5 ) := Cq+ i{s)Ci. If derivatives of R up to order k are bounded by constants 
Cl,... ,Ck, then t\D^ ^Jj < c((7o,..., (7fc) |Pp. Under reparametrisations of'y as in 
the first estimate remains unchanged, the second one scales with j. 

Proof. Our approach is to derive some differential equation for DgJ = VjJ, which 
has boundary values DsJ{s,0) = 0 and DsJ{s,t) = 0 for all s because J(s, 0) = 0 is 
constant in s and J(s,t) = 6{s) is the tangent of a geodesic. 

ad primum: Because J and T are coordinate vector fields, |i. 3 a| gives DgDtU = 
DtDgU + R{J,T)U for every vector field U, so we have 


DsJ = DsDtDtdsC = DtDsDtdsC + R{J, T)j 

= DtDtD.dsC + DtR{J, T)J + R{J, T)j 
= DlDM + R{J, T)J + R{j, T)J + 2i?( J, T)j 


whereas the (negative) left-hand side is, due to the Jacobi equation, 

-DJ = Ds{R{J, T)T) = {DsR){J, T)T + i?(L», J, T)T + R{J, j)T + R{J, T)j 

(note DgT = DtJ = J). From now on, we consider J and J as being part of the given 
data (which is allowed, as we have already sufficiently described their behaviour in 
I 6 . 3 P . So we have a linear second-order ODE for (7 :=DsJ'. 

U = AU + B, 


where both sides scale with 1 /under reparametrisation, and the norm of A is 
bounded through ||A|| < (7o|Tp(s). For ease of notation, we will thus assume that 
we consider a t-line with |r| (s) = 1 and rescale our results afterwards. By assumption 
on the smallness of r, 

\B\ < 2Ci|Jp + 5Co|J| |i| 

< 2 Ci|P |2 + 5Co|P|(i + |Cor)|P| 

< 15(7o.ii|Pr 

Now consider Fermi coordinates along c(s, • ) as in | 1 . 17 | to obtain an ode in Euclidean 
space. For any smooth vector field V = V^di, the covariant derivative in direction 
T = dtc is just VtE = V\di. Hence, our ode has the coordinate expression 

U\A = 

As we only need to know the values of U on a; = (t, 0,..., 0), this gives a euclidean 
differential equation for the components (7* of the same form as above. The claim on 
U is then contained in | 6 . 5 [ 
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ad sec.: With U = DgJ as above, we have DgJ = U + R{J, T)J and thus 

\Dj\ < \U\ + Col Jp < 45Co,i|t"P + 4Co|Fp. 

ad tertium: For higher s-derivatives, one can proceed by induction: The statement 
is true for k = 0,1, as we have shown above. By analogous computations, one can 
control ,, J by a linear second-order ODE, in which all lower derivatives might enter 
as “given data”. This data is bounded by a constant, and hence the solution will be 
bounded as well, q. e. d. 


Estimates for Normal Coordinates 


Situation. Fix some p € M and consider normal coordinates around p as in n. 121 6.7 

X : (u^,..., u”*) !->■ expp u^Ei 


for some orthonormal basis Ei of TpM. Recall from [TTT^ that the deviation of metric 
and connection from its Euclidean counterparts can be described by differential and 
Hessian of the exponential map. In the following, we let r d(p, •) be the geodesic 
distance to p. 


Lemma. Situation as above, Cgr^ < Then \gij — 5ij\ < C^r^ 


6.8 


Proof. By [ 3 ^ it suffices to consider i = j. So we only have to show |Idx— l| < 

Cor^. By means of |i.i 4 a[ this amounts to control | \du{expp)Ei\ — \Ei\ |. From fiTT^ we 
know that du{'expp)Ei is the terminal value J(l) of a Jacobi field with J(0) = 0 and 
j(0) = Ei. Now ||(i)7(expp)F;i| - \Ei\\ < |dc/(expp)£’i - PEi\ < jC^r"^ bylHnl and the 
squared norms thus cannot differ by more than 2(1 -|- • jC^r^ < O.SOJCor^ due to 

g. e. d. 

Lemma. Situation as above, Cgr^ < Then ||r|| < lOCgr + 5Cir^. 6.9 


Proof. Again, the case i = j is sufficient. Additionally, we will only prove the claim 
for r = 1. The correct scaling is then automatically enforced by 11.101 So let T,V G 
TpM be unit vectors, assume Cq < ^, and consider a variation of geodesics c{s,t) = 
exppt(T + sR). As the exponential mapping has no radial distortion, we may assume 
R _L T wthout loss of generality. This delivers us a Jacobi field J(s, •) = dsc{s, ■ ) for 
each s, and [i.iqb tells us that dT(expp)(F(u, u)) = V(iT(expp)(R, R) = VjJ(0,1) = 
Dsdsc{0, 1) for R = v^Ei. As observed in 16.61 the vector field U := DgdsiO, ■) along 
c(0, •) obeys the linear second-order ODE 


U = R{T, U)T + R{j, J)T + ZR{T, J)j + R{T, j)J + R{T, J)J -k {DsR){T, J)T 


where the obvious notation T for dtC has been used. So again we have U = AU + B 
with ||A|| < Co and |i?| < 2Ci| Jp-|-5Co|T| |T|, but this time as an initial-value problem 
with t/(0) = C/(0) = 0. Denoting the supremum over |i?| by b again, the norm |17| 
will be dominated by the solution of u = c^u + b, c= \/Cq, which is 2 ^e““‘(e°* — 1 )^, 
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which itself is smaller than |6 for ct < This means |V(iT(expp)(y, F)| < |&, and 
our task is to estimate B against V = J(0). 

From [S?T1 we get | J(t)| < (1 + j(0)| < (1 + f6)^|</(0)| for all t < 1. On the 

other hand,gives t\J{t)\ < (1 + ^Cot‘^)\J{t)\, and combining both leads us to the 
rought, but sufficient estimate \ j{t)\ < (1 + ^)(1 + ^)| j(0)|. So we have, as F = j(0) 
is of unit length, 

\B\ < (1 + + (1 + (6-9^) 

So far, we have only estimated the norm of VdT(expp)(F, F) = dT{expp){r{v,v)) by 
|5, and this needs to be compared to r(ri,n). By li.iGL the former is the value Z{1) 
of a Jacobi field Z along c(0, •), and the latter is Z{0). Using I^Til for Z, we obtain 
|^(1)| < (1 — ^)“^|Z(1)|, and by inserting this into |6.pa] we finally get 


|r(ei,e 0 | < 


T < 

-j TT^ — 

^ 16 


SQ I ]I_)2 

8 V-^ ^ 16/ 

-I _ TT^ 

16 


Cl 


|(1 + x )(1 + 76 )^ 

-| _ TT^ 

^ 16 


Co < 5Ci + lOCo, 


q. e. d. 


6.10 Remark, (a) As one can easily see in the proof, our numerical constants are by no means 
optimal. A sharper result, but with much more technical effort, has been given by IKaulI 
(liQvfil l. This author also deals with the case that the sectional curvature might be asym¬ 
metrically bounded between co and Co, whereas we are only interested in the simpler case 
Co = — Co- 


(b) Considering the Christoffel symbols as objects that store “derivative information” 
for the metric, the classical procedure of numerical analysis would have been to first 
estimate the Christoffel symbols and then integrate this to obtain a bound for the 
metric tensor. It is a specific property of the gij that they can be bounded by a 
right-hand side which includes fewer terms than the bound for their derivatives. 

(c) Under scaling of Mg as in ll.lll the estimate | 6 . 9 | scales like i, and 16.81 is scale-in¬ 
dependent. The assumptions in both propositions are scale-independent. 


(d) Regarding 1 1 . 3 b| rises the question if derivative s of R are actually needed to bo und 


||r|| . In fact they are needed in normal coordinates (IDE TURCK and KazdanIIiq8iL ex. 
2 . 3 ), but not in harmonic coordi nates, which would l ead to estimates that only depend 
on (Zoc. cit., thm. 2 . 1 ). As IBemelmans et all (Iiq8zi1I remarked, the metric g can 
be infinitesimally abridged by a short time of Ricci flow, and the new metric g has 
||V®.R|| < C'i(C'o) for all i. Furthermore, a bound on Vi? will be needed in 1 6.6 1 any wav, so 
we decided to take normal coordinates, which make it easier to give explicit numerical 
constants in the estimates. 


6.11 Conclusion. In a normal coordinate ball {B,u) of radius r with Cqu^ < 1 and 2 r < 
inj M, g and the Euclidean standard metric are equivalent, and 

\\V\g^u)-\V\,^\<Co\u\^\Vy, ||F(m)|| <10Co|m|+5CiKP. (6.11a) 
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Corollary. In a Fermi eoordinate tube of radius r with C^r'^ < 1 and 2r < inj M, g 
and the Euelidean standard metric are equivalent, and 

\\V\git,u) - < Co\u\^\V\p, ||r(i,w)|| < lOColul +5Ci|ix|2. ( 6 . 12 a) 

Lemma. Let g and g^ be two Riemannian metrics with ||t;|g — |f|ge| < £\v\gB, e < 1. 
Then the curve lengths and geodesic distances with respect to g and g^ fulfill 


\Lg{c)-Lge{c)\ <eLg.{c), \dg{p,q) - dg.{p,q)\ <edge{p,q). 

Proof. The first claim is proven in the obvious way by integrating | |c|g — |c|ge | < e|c|ge 
along c. The second claim is a combination with Lg(c) < Lg{c*^) and Lge{d^) < Lge[c) 
if c and c® are the distance-realising geodesics for g and 5 ® respectively, q. e. d. 


Approximation of the Metric 

Corollary. Let q be in a convex neighbourhood of p, i := d(p,q) with and 

let U £ TqM he an arbitrary direction. Then 

|VvAp-y| < fCo^Vv^l < lCo^^\v\, 

\KyXp\ < 50(Co+£Ci)£|4Fp. 


Here iTy is the orthogonal projection onto the orthogonal complement of Yp\g in TqM. 
Proof. Direct consequence of| 6 . 3 | and 16.61 together withfTT^ q. e. d. 


Remark. With \U\ instead of [ vry[/| , but wi th a smaller cons tant, the first claim is directly 
proven in I Josr and KarcherI 


ven in I Josr and KarcherI (lioSj . also cf. IKarch^IiqytI . a. 5.4). For the im provement. 
KaulI ( iQvSi l. An exact computation of Vd exp for symmetric spaces is given bv iFLETCHERl 




Lemma. Let A : V ^ V be an endomorphism of a normed vector space V with 
I id—A|| < e < 1. Then || id—< e/(l — e). 


Proof. By the Neumann series 


( Afidbon^ ex. 3 . 7 ): 


A i=^(id-A)*, so ||id-A i||<^e* = -^— 

i=0 ^ ^ 


q. e. d. 


Lemma. Let pq, ... ,pn be distinct points inside a convex hall of radius h and x he 
their barycentric mapping. If QC^h? < 1, then for a tangent vector v £ T^A at any 
A £ A and a as in Proposition f^7f\ 

\dxv — cr(t^)| < 2Coh^ |o'('f)|. 


6.12 


6.13 


6.14 


6.15 


6.16 
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Proof. By [5?7l d\xv = (AA)~^cr| a(i^) and hence \dxv — (t{v)\ < —id|| |o-(n)|. By 

| 6 . 14 [ one has |VyXi — V\< |C'od^(- ,pi)|F| for all tangent vectors or, in terms 
of operator norms, \VXi — id || < |C'od^(- ,pi) < Thus, as A • l„+i = 1 and 

A* > 0, 

IIAa - id II = ||A*(VX, - id)|| < |A*| ||VX, - id II < §Coh2. 

Now if QCoh? < 1, then 1 — > |, and the claim follows from | 6 .i^ q. e. d. 

Notation. We write a < 6 if there is some constant c which only depends on n 
such that a < cb (saying “a < 6 up to a constant.”). Equivalently, we will also write 
a = 0{b). We in particular remark that our suppressed constants do not depend on 
the geometry parameters. 


6.17 Theorem. Let pq, ... ,pn be distinct points inside a convex ball and x be their barycen- 
tric mapping. Let be the flat Riemannian metric on A induced by geodesic distances 
d(pi,pj). Suppose g® is (d, h)-small, SnCp/i^ < with from^flT^ Then it holds for 
tangent vectors v,w €. T\A 


\{x*g- g''){v,w)\ <Coh^\v\\w\. ( 6 . 17 a) 

The norms on the right-hand side can be either x*g or g^ norms, as both are equivalent. 

Proof. Note that the assumption on h includes the requirements of I 6 . 16 I and | 6 .i 3 [ 
Due to [3771 it suffices to show the claim for v = w. Consider a point A G A with 
image a = a::(A). We first compare x*g to the Euclidean metric of the simplex Sa = 
conv(Ailo) C TaM , and compare this metric to afterwards. 

Parametrise Sa in the canonical way over the unit simplex via x : Xci 1 —>■ X'Xi\a. 
Now clearly dx = a from [^ 77 ] The metric of Sa is the induced metric of the surronding 
vector space, namely g\a. Now use l6.i6l to get 

\ix*g\a){v,vy/'^ - {x*g\a){v,vy^^\ = \ \dxv\g\^ - \dxv\g\^\ 

< \dxv — dxv\g\^ < 2Coh^\dxv\g\^ = 2Coh'^\v\x»g- 

And of course, the same is true for the squared norms by 13.151 \{x*g — x*g){v,v)\ < 
6Coh^\v\ge. Hence we have successfully compared x*g to the euclidean metric of Sa- If 
we can show that s and Sa have almost equal metrics, we are done with [ 3 . 1 ^ 

The edge lengths of Sa are \Xi — Xj\g\^, and the edge lengths of s are the geodesic 
distance between pi = exp^(Ai) and pj = exp^(Aj). By | 6 .i 3 [ we have for their edge 
lengths iij and iij 

\iij - iij\ = \d{pi,pj) - \Xi - Xj\ \ < Coh'^ d{pi,pg) = Coh^iij, 
so 5 ® and g® match [ 3 T 1 BI with = Cohf, q. e. d. 

6.18 Corollary. hd| ■\p ;$ | • |g ^ d-l • |f2. 
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6. Approximation of the Geometry 


Definition. We say that points pi,... £ M lie in (i?, /i)-close position, if there 

is some p £ a;(A) such that defines a (i?, h)-small metric in the 

notation ofg^} (Note that this can only be if n < m.) 

Corollary. Each collection of points pQ,...,Pn in {•d,h)-close position that fulfill 
SnCgh^ < a„ defines an injective harycentric map. 


Remark. As x*g and g are equivalent metrics, there is a self-a djoint automorphism J of 
TaA such that x*g{v,w) = g'^{Jv,w), as has been empoyed by HoLSiLand_STERN| ([201 j . 
thm. .*^■ 8 ). Fo r a com parison to the metric distortion tensor A of lWARDETZK^ 2nofih and the 
Ah of lDziu^ (lioSStl et ah, see jll.t^] 


Approximation of Covariant Derivatives 


Remark. The second-order approximation qualities of a parametrisation would usu¬ 
ally be measured by bounds on the Christoffel symbols. However, our definition of 
is not exactly the inverse matrix of qa, so the us ual definition I i. 2 bl would not work. 

Instead, we employ the idea from KaulI ( lovbll to consider the operator T = ® — 

V® , which would have the coordinate expression T{V,W) = dk in a usual 

n-dimensional chart. Recall from I 6 . 7 ] that ® is defined by dxV^ = Vdxvdxw 
We suppress the g subscript for norms. 


Theorem. Situation as in \ 6 . 1 y[ Then iVdxlAgs.Mg < C^o.ih. 

Proof. Due to [g771 it suffices to show the theorem for v = w. Similar to ~ id || ^ 
Cgh^ we have, as the u* sum up to zero, 


|A„(R)| = \PVvX,-Y,EV\ < \P\\yvX,-V\ < ^\vUiCoh^\V\g. 
Now we use | 5 . 7 b| and again that ||R^^|| 1 + Coh^: 

Y7p^^|Vda;(u,u)| < 2 \Ay{dxv)\ + \XVl^^^dxvX^\ 

< Coh'^\v\e^ \dxv\ + Coph\dxv\'^ < Co,ih'^\v\i2 |u|, 


where V'^Xi has been estimated by | 6 .i 4 [ q. e. d. 

Corollary. < CQ ih\v\ |w|. 

Proof. It suffices to consider v and w with constant coefficients, so V® w = 0. By defi¬ 
nition of ® and fTTTl |V^ ^w\x»g = \Xdxvd^'>^\g = |Vda;(u,w;)|, and so the preceding 
theorem applies, q. e. d. 


Corollary. For A,// £ A, it holds \d\xv — Pdfj,xv\ < Cq ih \\ — g] |u|. 
Proof. By the fundamental theorem [TTT^^ 

d\xv = PdgXv + J PVdx^dxv = PdgXv + J PVdx{fi,v) 


for a curve 7 : A ^ /r. 


q. e. d. 


6.19 

6.20 

6.21 

6.22 


6.23 

6.24 
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B. Main Constructions 


6.25 


Remark, (a) The piecewise flat metric gfj = is nothing more than the “most 

natural candidate” for a constant Riemannian metric. Any other 5® that is a sec¬ 
ond-order approximation of g|o for some a G s would give the same result. The par- 
ticulary interesting observation is that g can be approximated up to second order by 
a piecewise constant metric, whereas an arbitrary function would require a piecewise 
linear function for a similar approximation order. 


(b) The convergence result for the connection does not mean that if Mg is triangu¬ 
lated over a sequence of finer and finer simplicial complices rKh, the connections of x*f^g 
and gf^ would converge. In fact, gf^ would be always piecewise flat, so the connection 
would vanish and hence can never approximate the connection of a curved Mg. This 
global impossibility is consistent with our simplex-wise convergence result because the 
connection for (7® on two adjacent simplices cannot be compared to each other, as the 
metric is not continuous across the simplex boundary. The connection V® can hence¬ 
forth not be connected to the derivative of globally, but only in those matters which 


make sense in this situation, e. g. higher derivatives of real-valued functions as in 7.2b 

(c) The convergence of Lipschitz-Killing curvatures from Cheeger et al. (|iq 8 /iI 1 ap¬ 
plies for our situation, although their triangulation is defined in a slighly different way, 
see |8.c[b | It is a convergence in measure of rate Por submanifolds of Euclidean 
space, Cohen-Steiner and MorvanI (200^ give a convergence order h in measure, 
but we did not check if their arguments can be carried over to our setting. 


(d) The metric approximation result is similar to, and of the same order as the usual one us¬ 
ing orthogonal projection of a triangular surface onto some nearby smooth surface (cf. IDziukI 
lioSSh . We will reproduce this conventional approach for the approximation of submanifolds 
in section rm 


7. Approximation of Functions 

Goal. In this section, we want to apply the Karcher mean construction and the results 
from the previous section to functions between manifolds: First, we consider the case 
of functions (A, x*g) —?► R, where the preimage is to be approximated by (A, g®). After 
that, we will consider the case of approximation in the image, which means we will 
interpolate some function y : A —>■ M by the Karcher simplex parametrisation x with 
prescribed values x{ei) = y{ei). 

Approximation in the Preimage 

7.1 Situation. Suppose x*g and y® are two Riemannian metrics on A, that y® is flat 
and that I6.17I for the metric as well as I6.23I for the Christoffel operator holds. Vector 
and operator norms, if not explicitely qualified, are taken with respect to one of these 
equivalent norms. 

7.2 Proposition. Situation as in \Y 7 T\ For given smooth u : A — > R, 

Igrad"^ ® M — grad® u| < Co/i^lgrad^ ® u|, (7-2a) 

||V^*®du - V^‘du\\ < C'ojh\du\ (7.2b) 
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7 - Approximation of Functions 


Remark. It is easier to estimate the operator norm ||V“ ^du — V® du\\^ although we 
will actually need the induced norm |V® ^du — V® du\ for bilinear forms or bi-covec¬ 
tors. Recall that the equivalence constant for these two norms only depends on the 
dimension, which will be neglected as usual, so || • ||g < | • |^ < || • ||g on any tensor 
bundle over TM. 

Proof, ad primum: Represent du = UidX'. In the notation of |3.i3| we have grad“ ® u = 
Q^^Uidj and grad® u = {Q^Y^Uidj. So with u = (ui,..., u„), 

|grad“ ® u - grad® u|g = E{Q - Q®)u • (Q - 
< {C'^h^fEQu-Qu 
= {C'oh^YQu ■ u = {C'qIY \du\f. 

ad sec.: Bv ii.2ai the difference between two connections only depends on their Christof- 
fel symbols. Extend the vectors v,w G T\A to vector fields with constant coefficients. 
As 5® is flat, this gives = 0 and V®'w = 0 . Now bv li.8al 

(V®'du - V^*®du)(u,u;) = tc) - du{vf^w) = du((V®" - V“*®)„w) 


and together with |6.23} 

|(V“*®du-V®'du)(u,u;)| = \du{vfsw - w)\ < \du\ |r(u,u;)| 

< \du\CQ ih\v\ |w|. 


q. e. d. 

Proposition. Situation as in \Y 7 I\ The -norms, k = 0 , 1 , 2 , with respect to x*g 
and are equivalent for every p G [1; c»[.' 

l«irp(Ax*g) = l“ir.>(Ag<=)(l + 0 {C'^^h^)), (7.3a) 

M«lLP(Ax*g) = M“lLP(Ag<=)(l +C’(C'oCp/l^)), (7-3b) 

M^lwi>P(Ax*g) “ M'*^lwi’P(Ags)(^ 0(170 j^Cph)), (7'3<^) 


with Cp from \gTT^ The same holds, without power p and factor Cp, for the \N^’°° norms. 

Remark. Note that the estimates speak about | • 1 ^* instead of | • |. This means that 
the estimates become worse for p —>■ oo. Therefore, an additional argument for the case 
p = oo is needed. 


Proof. Case k = 0 ; The Lebesgue norms on Ax*g and Ap® only differ by their volume 
elements G and G®, which fulfills the claimed equivalences thanks to |3.2q[ So 


I lure- 1 1 


l\pG<^ 


<C',h^ 


l\pG. 


In the norm, there is nothing to show, as both norms agree. 


7-3 
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B. Main Constructions 


Case k = 1 : Here an approximation of the volume element and the gradient norm 
enter: 


ly g{du,duy^^G — J {du,duY^'^G' 

— \j g{du,duy^^{G — G^) + Cp^J {g — gy{du,duy^^G‘^ 


< 


G'oCph^ J g{du,duy/^G, 


because Cp > 1 . For the norm of du, it suffices to observe that if |dA*u|pe is maximal 
among all A £ A, then \dx*u\gt> ^ (1 + 0 (C'g/i^))|dA* w|g < (1 + O(CQh^)) max^ \dxu\g. 

Case k = 2 : We do not have an estimate of our usual form \x — y\ < £\x\ for the 
Hessian, but the proof of 13.151 also admits this situation: 


||V9dw|P - |V9"du|P| < CplV^duYg-^WV^dulg - \Vs‘du\g 

<Cp|V®du|P-i|du|p||r|| 
<Cp{^\y^du\P + ^\du\Pg)\\T\\ 
<cp{\v^du\p + \du\p)\\ri 


thanks to Young’s inequality ( Alt 2006I eqn. 1-11). Now one needs approximations of 


the volume form, the norm on covectors and bi-covectors from 13.191 as well as of the 
Hessian: 


I \V^du\PgG- I \y^‘dufg.G^\ 

< J\\Vadu\Pg-\Vs‘du\Pg\G + J\\VP‘du\Pg - \VP‘dufgelG + ^ |V^^|^s(G - ) 

< J \\V3du\Pg-\V3'‘du\Pg\G + C’oCph^ J du\PgG 

< {C'oCph^ + C'o^.cph) j\VPdu\PG + G^^iCp/i J \du\PgG, 


q. e. d. 

7.4 Theorem. Situation as in [Y 7 T\ For a function it : A —> R, let : A —> R &e its 
Lagrange interpolation, that means Uh is linear and Uh(ei) = u(ei). Then 

\u- U/t|L~(A) + h\d{u - M?i)|l->(A) ||V®'dM||L=o(Ape). 

The right-hand side can be replaced by h^'d~^{l -f Gg ih)||V“ . 

Proof. If we were only interested in this interpolation of real-valued functions, the 
easiest method of proof would be to use the interpolation estimates in Euclidean space. 
But when we come to mappings into a second manifold in [77^ these methods would not 
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7 - Approximation of Functions 


be applicable without further work. Therefore we decided to use a more “geometric” 
approach. 

ad primum: Let /i £ A be an arbitrary point, consider the tangent e^- = Cj — Ci to 
the geodesic 7^- : Cj Ci and 


n : A !-)■ {dxu - d\Uh)eij. 

This scalar-valued function has a zero along the geodesic , because ri o 7^ is the 
map t I—> [du — duhjifiij) = ~ '^h\{'yit)), and \u — Uh\ is zero at both endpoints 

of 'jij. Let ly € A he the position of this extremum. 

Now let 7 be the geodesic v p, and if{t) := ri{'y{t)) = {d.y(^t)U — d,y(^t')Uh)eij. Then 

iPit) =ff®(grad®'(M-u/,),V? e^)-k grad®° (u - e^). 

The first summand vanishes because is parallel with respect to (7®, and the second 
one is V® d{u — Uh){eij,j) due to [T?7| So 

\ip{t)\ < h||v®'d(u - Uh)\\g'^ Iflgo for all t, (7.4a) 

and because is linear, = 0 . Hence < /|■^/'(s)| ds < du\\_^{^^..g^y 

If Ek form an orthonormal basis, then \du\'^ = 'J 2 du{Ek)'^ Because of |3.6[ the Ek 
have an expression in the e^- with coefficients smaller than l/dh, which gives 

\du - dUh\gc\^ < h'd~^ ||V(iM||L=o(A,g'=)- 

As p was chosen arbitrarily, this holds for every point in A. 

ad sec.: Now consider a new arbitrary point /r £ A, the function 

ro : A |u(A) - m?i(A)P 

and a geodesic 7 : A for some vertex et of A. Then let (p{t) := ro(7(t)). 

As To vanishes at the interpolation points, we have 75(0) = 0 , and everywhere 
\(p{t)\ = \d{u-Uh)j\ < \d{u - Uh)\g’‘\'y\g'> < /i'd“^|7|ge ||Vs'du||L~(A.g«) and thus 
< f |(/3(s)| ds < ||V 9 “du||Loo(A,ge). 

ad tertium: The last statement is a direct application of [77^ q. e. d. 


Corollary. The same result also applies for the norms: 

\u - Uh\LP(A) +h\d{u- U;j)|lp(A) ^ ||V® du||LP(Ag<=)- 

The right-hand side can be replaced by h^'&~^{l Cq ih) ||V® dM||LP(Aa;*g)• 

Proof. Only the estimate |7.4a| has to be refin ed by the “Holder i-trick”, a common 
application of Holder’s inequality ( Alt! [20^ lemma 1.10): Suppose some function 
a £ L°°(A) is estimated pointwise by |a(A)| < where the integration path 
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B. Main Constructions 


7[A] : eo a is of size h. T hen as in the m ost “basic” proof (there are others, cf. l2.iob|l 
of the Poincare inequality (jADAMsIliQvd . sec. 6.26), 


= / |aP < 


ILp(A) 


61 


A 7[A] 




(7-5a) 


A 7[A] 


A 7[A] 


Then compute the A integral by first integrating over the subsimplex Aq opposite to 
the vertex eo and then over the ray cq ^ p G Aq. Then A = tp + (1 — t)eo for some 
t between 0 and 1 , and each function c G L°“(]R) with c > 0 fulhlls c(s) ds dt < 
r fg c(s) ds, we have 

I I bP<h I bC (7.5b) 

A 7[A] a 

Then observe ^ = p—1, so we have |a|Lp|.^j < for such a function a. As there 

does not occur any term in the hnal estimate, it remains valid for a,b G L^(A), 

q. e. d. 


Approximation in the Image 


Remark. For curves in M, there are already interp olati on estimates fo r high-order (quasi-) 
interpolation methods bv lWALLNER an d Dyn i|20ofJ') an d iGROHsI (l20iqll . 

During the finishing of this thesis, Grohs et al\ ll2nT?l have given a very elaborate estimate 
for higher-order “polynomial” interpolation using the Karcher mean construction. We decided 
to nevertheless publish our proof here, as we hope that our approach gives more geometric 
intuition, involves simpler constants, and is used in sections ITTIfT^ 


7.6 Situation. In the following, we assume that A carries a {d, 6,)-small Euclidean metric 
(which is not assumed to come from geodesic distances in M). We consider a smooth 
function y : Ag® —>■ Mg (and assume that y{A) lies in a convex ball of radius r as in 
16 .Ill with < 1 ) and dehne x to be the barycentric mapping with respect to the 
vertices y{ei). We will usually write x and y instead of x(A) and y{\). 


7.7 Lemma. Situation as in \ y. 6 \ Let P be the parallel transport TyM —^ T^M. Consider 
d(x,2/) and i 3 ?{x^y) as functions A —^ R. Then 


did^ix, y))v = 2 g\^{Xy, {dx - Pdy)v), 
d{d{x,y))v= g\x{Yy,{dx - Pdy)v), 


with Xp, Yp as in \ 1.221 

Proof. From[T 32 l we know the gradients of d and d^ if only one of the two arguments 
is varying. Then for (/? : M x M —>■ R, {p,q) d^(p,g), we have for tangent vectors 

V G TpM and W G TqM that 

difiiV, W) = g\p{V,Xq) + g\q{W, Xp). 


48 






















7 - Approximation of Functions 


Now d^(a;,2/) is the concatenation of the map A i—^ {x,y), which has derivative v i—>• 
{dx V, dy v) with tp, so 

d{d^{x,y))v = g\a;{Xy,dxv) + g\y{Xa;,dyv) 

As Xy is the starting tangent of the geodesic x y parametrised over [ 0 ; 1 ], we have 
PXx = —Xy, and P is an isometry, so g\y{Xx,dyv) = —g\x{Xy, Pdy v), q. e. d. 

Lemma. Let c{s,t) be a smooth variation of curves c(s, •), and let Pg’‘^ : Tc[s,a)AI —>■ 
Tc(s,b)Ai be the parallel transport along these curves. Then 

b 

jg^pb,a^ J pb,tj^^g^ g^^pt,a 

a 

(note that the integrand is always a linear map the integration 

is therefore defined without problems) and hence ||Z 3 sPs’“|| < Cq g |ds|. 

Proof. Consider a vector field F(s) G Tc(s,a)M, and let V{s,t) := P)’°-V{s). As 
first step, observe that Ds{Ps’°'V{s)) = Pg’°'{DsV{s)) + {DsP^’°')V{s). This formula 
seems obvious, but actually requires a little argumentation: It symbolically resembles 
V{AV) = iyA)V + AiyV) for linear bundle maps A from FTT^ but as P mediates 
between different tangent spaces for preimage and image, the V operator is not the 
same on both sides. Instead, consider the function / : c(s, a) i—> c(s, b) between the a 
and the b isoline A and B respectively. The parallel transport from a to 6 is a mapping 
TxM —>■ Tf(^x)M and hence an element of TM\a 0 f*{TM\B). As in li.bbl the induced 
connection on this bundle is given by 

Va. {oj ® rV) = (Va,w) 0 /V + a; 0 fX^nogV, 

and indeed df{dsc{s, a)) = dsc{s, b), giving the Leibniz rule for PV. On the other hand, 
the fundamental theorem of calculus gives 

b 

Ds{V{s,b))=Py{Dy{s,a)) + I P^’\DtDsV{s,t)) dt. 

a 

Because C(s, •) is parallel, the vector field in the integrand is 

DtDsV{s, t) = DsDtV{s, t) + R{dt, ds)V{s, t) = 0 + R{dt, ds)PyV{s). 

Because C(s) is independent of t, it can be pulled out of the integral.—The second 
claims results from 

b 

WDsP^’l < I l|i?ll 1^*1 WPs’^W lln‘’“ll dt 

a 

and IPs’* I = 1 everywhere because parallel transport is isometric, q. e. d. 


7.8 
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7-9 


Remark. Generally, it is well-known that the curvature tensor can be characterised as 
infi nitesimal version o f holonomy, i. e. the parallel transport alon g a clo s ed cu rve (see 
e. g. Peters^ I2006L sec 8.6). We found this specific version in IRaniI ( 20oqL lemma 
.■^.2.2). The estimate can obv iously be sharpened by replacing | 9 t| | 3 s| by | 9 t A 9 ^ 1 , see 
Buser and Karch^ (|iq8i[ 6.2.1). 


Lemma. Situation as m |7.6) Then ifd{x,y) < p everywhere in A, we have at every 
vertex Ci 


\\de,x - de^yjAg^Mg ^ hd-^{\\Wdy\\ L°°(Ag<‘,Mg) ^0,1 p\\dy\\iaa ■ 


Proof. First, consider v to be an edge vector Cj — Ci, so c : 1 Ci + tv parametrises the 
ij edge over [ 0 ; 1 ]. Then choose Fermi coordinates (t, ..., u™) along an arclength- 

parametrised version of 7 := cc o c. As 7 itself is not parametrised by arclength, it has 
coordinates 7(t) = (^, 0 ,..., 0 ) with a = d(j)i,pj). The image of c under y is another 
curve S which intersects 7 at pi and pj , so 


, 5 ( 0 ) =7(0), 5(1) = 7(1). 

By the intermediate value theorem, each component 7* — 5 ® must have a zero at some 
T® G [ 0 ; 1 ]. As Dtf = 0 and F^^ = 0 along 7, the second derivatives Yu of the compo¬ 
nents vanish, too. So 


|(f-< 5 *)( 0 )| < I | 5 :,,|dt<r®| 5 :,*|L==([ 0 ;i]). 

0 

By [T?7l DtS = Wdy{v,v), and together with Dt 5 = {Sf^. + from fT:^ we 

have 

\S,tt\g < \S/dy{v,v)\g + \dyv\l max||r|| 

Bv l6.i2al we have \dyv\ = | 5 |g ^ (1 + Co,ipY\d\p: which means that both norms are 
equivalent for small p. Similarly, max||r|| < C^pp. Together with \v\ge = a < h, we 
have 

\{dx - dy)v\g^^_ <h[\\Vdy\\ L°°([0;a]g7Mg) + C'o.l P || dy || L~ (cgGMg) ) bIs'- 

This shows the claimed estimate for edge vectors. And some general v that is not 
tangent to an edge can be represented as linear combination of edge tangents Cj, and 
all coefficients n® are estimated from above by \v\g/'d up to a constant, q. e. d. 

7.10 Remark, (a) For triangles, the fullness parameter d controls the minimum angle at 
each vertex. This is exactly the parameter that enters in the last argument, so there 
is a direct geometry meaning of the factor 1?“^. 

(b) There are also coordinate-free methods to prove [7:^ but we did not find any 
method that is “so intrinsic” that no curvature term like C'o,i/ip||c??/|| comes in. For 
example, one could transport 5 and 7 both to the vertex p = y{ei) and do all compar¬ 
isons there. Then the estimate 5 _tt — 5 is not needed anymore, but some VP and the 
holonomy _ pi{t),ppp,&d) j^ave to be estimated by I7.8I and 113.41 
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(c) The term in parentheses on the right-hand side of [7?^ is what Ic^OHS et aD ll20iql i esti¬ 
mate by their “smoothness descriptor”. Our computation shows that the nonlinear lower-order 
term \dy^ only enters with an additional distance factor p. 

Proposition. Situation as in Then ifd{x,y) < p everywhere in A, T-n 

||dx - Pdy||L=o(Ag%Mg) < h'd~^ (llVdy - PVda;||L~(cg«,Mg) + C'o,l/o|M2/|lL-(cgSMg)) ■ 

Proof. Let us prove the claim at some p = x{p). Consider some vector v € TA, and 
let V := (dx — Pdy)v. Along a geodesic ^ = x o c pi p, which comes from a curve 
c : Ci p in A, we have by the fundamental theorem p-igai 

1 

l/|p = pi’Vlp, +1 dt, 

0 

where P is the parallel transport along 7 (not to be confused with the parallel transport 
P along geodesics y x). Inside the integral, we have DtV = VdxcV = S/dxcdxv — 

^dxc{Pdy)v. As in the proof oCT define a mapping / : a;(A) 1—?► y{X). Then df{dx w) = 
dyw and hence X/dxc{Pdy)v = PX/dycdyv + iydxcP){dyv). Together, this gives 

DtV = Vdxcdxv - PVdycdyv - {VdxcP)dyv 
= Vdx[c,v) - PVdy[c,v) - {VdxcP)dyv. 

By |7.8[ we have \S/dxcP\ A Cop max 15 ^ |, where p is again the maximum distance 
between x and y, and dg is the vector field defined in the proof above and has values 
dx c and dy c at the endpoints x and y. Thus 

\DtV\g < ||Vdj/-PVda;|||c|ge \v\g^ -f Copmax \\dy\\ |z;|ge 

Now observe | 9 s| < |dyc|, which gives 

lAI^Ig < h{\\Vdy - PVdx\\ + CopIdylP) |u|ge |c|ge. 

I^IsIp — l^lglp- + max|PtP|, the claim is proven with help of [77^ for the initial 
value at pi, q. e. d. 

Proposition. Situation as in [7^ IfCopid ^h\\dy\\^ is small, then 7.12 

|d(a;,p)|L.»(A) < \\Vdy - PVda;||L~(cge,Mg)- 

Proof. Consider any point A G A and a geodesic c : A. Then, with [777I 

d(a;(A),y(A)) = y d{d{x,y))c< j\{dx - Pdy)c\g < h\dx - Pdy\, 
everywhere, and this norm is estimated by 17.111 There is a constant a such that 

and the assumption means that the fraction is greater than, say, ^, q. e. d. 


51 










B. Main Constructions 


7.13 Remark, (a) The smallness assumption on Co,IdyP is reasonable because the 
interesting situation is when the domain is decomposed into finer and finer simplicial 
complexes. In this case h —?► 0 , whereas (given that the subdivision is performed intel¬ 
ligently) d can be bounded from below independent of h, and Co,i as well as ||c??/|| are 
independent of this refinement (here it is important that ||dj/|| is taken with respect to 

on A, not 

(b) The estimates are scale-aware: When A is scaled by and M is scaled 

by 5 = g?g like in li.iobl then both sides of the estimate s |7.11| and | 7 . 1 : 5 | scale similarly, 
namely like ^ and like /i respectively: In fact, f = pr, h = vh, ||i?||g = ^||di||g, and 
derivatives of x and v scale like in fiTiobl ||da;||ge ^ = p||c?a:||ge _g, ||V(ia:||ge,g = ■;^||Vdj:||ge,g 
and similar for y. 

7.14 Conclusion. Taking \j 7 T 2 \ and \'y.ii\ tog ether, we get 

|d(a;,2/)|L~(A) + h\dx - T’d?/||L-(Ag«,Mg) < \\Vdy - CV(ia:||L~(AgyMg)- 

7.15 Theorem. Let N and M be Riemannian manifolds with curvature bounds Cq and 
Cl as usual, and y : N ^ M be a given smooth function. Suppose pQ,...,pn S N 
are given points in (d, h)-close position with h so small that their barycentric mapping 
A —> s is injective, where s G N is the Karcher simplex with respect to vertices pt, 
and furthermore suppose that the barycentric mapping A —>■ M with respect to vertices 
y{pi) is well-defined. Then */Cg 1 h||dy||^00 is small in comparison to the dimensions, 
there is a function yh '■ s ^ M interpolating y at the pi, with 

|d(2/;t,y)|L»(s) +h\dyh - Pdy\i_^(sM) ~ h'^'d~^\^dyh - CVd?/||L=o(s,M)- 

Proof. By |6.i9( there is a bijective barycentric mapping xn : A —>■ s with Cj ^ Pi. If 
xm : a —>■ M is the barycentric mapping with respect to vertices y{pi), which have 
distance less than h \\dyl set yh '■= Xm o The the estimate is a combination of 
[7^41 and q. e. d. 

7.16 Remark, (a) One could have proven the intermediate estimates [77^ and 17.111 for 
scaled versions of A and M, for example with instead of g^ or diamM < 1 . But we 
did not consider the situation above complicated enough to justify a separate scaling 
argument. But if one likes, the argument obviously could have been executed for A 
and M having both unit size. Then is the equivalent of the usual “transformation 
from the reference element”. 

(b) The step from the estimate 17.141 to an estimate works exactly as in [7^ so 

we save paper by not repeating all the integrals. 

(c) The estimate could be considered as “incomplete work”, as the right-hand side 
still contains a Vda; term. We decided not to estimate it bv I6.22I to make clear that 
the right-hand side tends to zero if y is an “almost barycentric” map. 

(d) For a “higher-order” interpolation, IGrohs et all pnioh use basis functions : 

A —^ M of higher order, not just 7)* (A) = A* as we did, that fulfill = 1 and 
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8. The Karcher-Delaunay Triangulation 


= Sj for control points € A. Then the interpolation with respect to points 
Pi = y{pi) is the minimiser of (p* •). 

If one chooses the such that there are k + 1 control points on each edge, as is 
usually done (e. g. for the quadratic basis functions A®( 2 A® — 1 ) and 4 A®A'^, i ^ j), then 
[7:^ and 17.111 can obviously be interated to give 

\\V^dx - PV^dyl < \\V’^dx - W'^Pdy\\ + curvature terms. 

The only point of difficulty is to show that ||V^da;|| is actually bounded by ||dj/|| and 
the geometry, which needs quite some computation, as it involves many derivative 
norms but is provable along the same lines. 


8. The Karcher-Delaunay Triangulation 


Notation. The term “triangulation” is used in various senses in (discrete) geometry, 
topology, and computational mathematics. We will use it only in the topological mean¬ 
ing as a map M for some simplicial complex M.. (To obtain the correct homology, 

one usually requires this map to be bijective. We will construct this map and give 
conditions for its injectivity, but we will also call it a triangulation without these con¬ 
ditions.) The partition of a space M into topological disks will in contrast be referred 
to as a “tesselation”. 


Goal. In this section, we want to explore how a simplicial structure can be imposed 
on a “ 5 -dense” point set in a manifold. Conversely, if the simplicial structure and the 
vertex set are given, the resulting triangulation weill be considered in section im 


Situation. In the following, the usual assumption from [TT^ that M is closed (i. e. 8.1 
complete and without boundary) is essential. Let V C M be a set of finitely many, but 
at least m points in M and <5 > 0 be such that each ^-ball in M contains at least one 
point from V. We say that V is d-dense in M. 


Definition ( Leibon and Letscher 2nnrJ) . Situation as in 18 .il Let p € V. The 8.2 
Voronoi cell of p is the set of points in M that are nearer to p than to any other 
point in V: 

Vp := {a G M: d{a,p) < d(a, q) for all g G V} 


These sets cover M, overlapping only on their boundaries. The cover {Vpi p G V} is 
called the Voronoi tesselation of M with vertex set V. 

The bisector Bpg of p and q G V is the set of points which have equal distance to 
p and q, but larger distance to all other points in V: 


Bpq := {a G M: d{a,p) = d(a, q) < d(a, r) for all r G V} 

Obviously, Bpg = Vp O Vq. Similarly, the bisector B^ of a set s C V is defined as 
fjpGs ^p- ^ generic if each non-empty bisector B^ is a disk-type 

submanifold (with boundary) of codimension |s| — 1 and B^ is empty for |s| > m -I- 1. 
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8.3 Proposition. Situation as in \H.i\. p £ V. Then Vp has diameter less than 26 . If 
25 < cvr M, then Vp is a topological ball. 

Proof, ad primum: No point in Vp has distance greater than 5 from p, so the claim is 
simply the triangle ineqnality. 

ad sec.: If geodesics starting from p are nnique, then Vp is star-shaped, cf. [5TTI 

q. e. d. 

8.4 Remark, (a) Let 5 C V with nonempty B^. Locally, B^ is a smooth submanifold, 
and its tangent space in a non-boundary point a G B^ is 


{W G TaM : g{W, Xp - Xg) = 0 for all p,qG s}. 


where Xp and Xq are the gradients of squared distances as in [5^ In fact, a curve 7 
with image in B^ fulfills d^(7(t),p) — d'^(p/(t),q) = 0 everywhere, which has derivative 
5(71 ~ ) ■ 

(b) Note that Vp will in general no t be convex b ecause Vp and Vq could only be both 
convex if Bpq were totally geodesic. IBeem (i^XS) showed that all bisectors are totally 
geodesic if and only if M has constant sectional curvature. 

(c) Gener ally, the p r opert ies of topological spheres in Riemannian manifolds are 
treated by KarcherI (1368): A topological sphere that does not meet its cut locus 
cuts M in two open sets, some “interior” ball and some “outside”. A set B C M is 
convex if and only if each p G dB has a “geodesic support plane”, that is a subspace 
Hp C TpM of codimension 1, such that all starting directions (expp)“^a for points 
a G B lie on the same side of Hp. 

(d) IBoissonnat et 'ai. 1 201 ih remark that genericity of a point set, which can be 
achieved in Euclidean space by an arbitrarily small perturbation of a degenerate point 
set, is not always removable by infinitesimally small changes of V. For this reason, we 
assume a generic V and disregard the question of sharp conditions that ensure this. The 
problem is currently treated in detail by Dyer and Wintraecken (Rijksuniversiteit 
Groningen). 


8.5 Proposition. Situation as in \ 8 .i \ with generic V and 26 < cvrM. Then ^ := {s C 
V: Bs is non-empty, 5 G £ = 0 ,...,m, define a regular simplicial complex 

without boundary, called the Delaunay complex for M, with vertex set = V. 


Proof. The only property for a simplicial complex, that some I C s with cardinality k 
is contained in the set of /c-simplices is clear, because Bi C B^ for i C s. 

It remains to show that A is regular and has no boundary. An (rn — l)-simplex 
t cannot be part of more than two m-simplices, because a non-constraint bisector 
{a G M: d(a,p) = d(a, q)} divides M into two distinct sets. On the other hand, there 
cannot be only one m-simplex containing t, because boundaries of the Voronoi cells can 
only occur where two cells meet if M has no boundary for itself. And each Gsimplex 
belongs to an (£-|- l)-simplex: Let t G Because Bi is a topological disk of dimension 
n — £, it must have a boundary, which in turn can only consist of bisectors B^ with 
t C s, q. e. d. 
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g. A Piecewise Constant Interpolation of dec 


Remark. Note that the following situations are ruled out by our assumptions: 8.6 


(a) an m-dimensional sphere with V = {po,... ,Pm}, because By would consist of 
two points of equidistance, which is not a 0 -ball, but a 0 -disk. However, the definition 
of Voronoi regions would be feasible, but its dual would consist of two m-simplices 
with the same vertices, and our notation does not allow to distinguish between them. 


(b) m + 2 equidistant points V = {po,... ,Pm+i} in an m-dimensional manifold, 
because this V is not generic: In fact. By would be the point of equidistance, but this 
set should be empty, as |H| > m -I- 1 . 


(c) The counterexample of lBoiSSONNAT et al. ( 201: 
B{p,u,v,w} is not empty. 


pp. 38sqq.), because the bisector 


Definition. Situation as in IS.il with generic V and 2 S < cvrM. Let A be the complex 8.7 
from |8.5| For e € .ft”, let Xc be the mapping from [^ 7^1 As ^dr-f only depends on f for 
fee, this piecewise definition gives a well-defined mapping a; : r.ft —>■ M, called the 
Karcher Delaunay triangulation of M with vertex set V. 


Proposition. The Karcher-Delaunay triangulation is indeed a triangulation in the 8.8 
usual sense: Situation as in un\ with generic V and 26 < cvrM. If 6 is so small that 
the requirements of \6.ig\ are met on each Karcher simplex, then x is hijective. 

Proof. The map x is surjective because its image is non-empty and has no boundary 
in M. By |6.i9| each x^ is injective, and as the Karcher simplices do not overlap except 
on their boundaries, so is x, q. e. d. 


Remark, (a) If M is not closed, but compact and with boundary, the construction 8.9 
is of course feasible, but will only be bijective if there are also points on the boundary 
and the boundary is aligned with their Karcher-Delaunay triangulation. 

(b) The construction of lCuEEGER et al\ (lioSJ l seems similar, but (of course) does not use 

our barycentric mapping. It starts with a triangulation x : rft —>■ M, considers finer and finer 
subdivisions s : rS! —>■ rft of the complex, and then compares the metric (s o xYg on rS! 
to the piecewise flat metric induced by edge lengths lij = d(a; o s{ri),x o s{rj)), for edges 
ij € ill i-ii® subdivided complex. 

(c) We know, however, that IBurago et al\ (l20iql l state that “it is now clear that in di¬ 
mensions beyond three polyhedral structures are too rigid to serve as discrete models of 
Riemannian spaces with curvature bounds”, but nevertheless there will certainly be rigidity 
results for spaces of piecewise constant curvature without counterparts in the smooth cate¬ 
gory, we are not convinced that the references they give support this statement in its full 
generality. 


9 . A Piecewise Constant Interpolation of dec 

Goal. As second main construction of this thesis, we will now give an interpretation 
of the discrete exterior calculus as piecewise constant differential forms, which turns 
variational problems in the simplicial cohomology (C^, d* ) into problems in a complex 
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{P~^n^,d). The main question will be the connection between d and the usual exte¬ 
rior derivative d on The introductional definitions are t he basics of simpli cial 

homology as the y can be found in any topology textbook, e. g. 

HatcherI (I2001II . 


MunkresI (|iQ8il or 


Discrete Exterior Calculus (dec) 

9.1 Definition. Let i? be a ring with neutrals 0 and 1 , and let .fior be a regular n-dimensio- 
nal oriented simplicial complex. For any simplex s £ let Xs '■ -^or i? be defined 
by Xb{s) = 1 and Xb{s') = 0_for any s' ^ s. 

Consider the i?-module Ck{^) that is spanned by all s G Let Ck{^), the 
space of /c-chains over (with coefficients in R), be the quotient of Cki^) under the 
identification of Xs- and —Xs- Its dual space the i?-module of all homomor- 

phisms Ck{^) —t R, is called the space of fe-cochains over ^ (with coefficients in R). 
Let be the generators of dual to 5 ^, that means f^{xs) = 1 and f^ixs') = 0 

for s 7^ s'. In the following, we will only use i? = R. 

The boundary operator is the linear map C'fe(.fi) —>■ Ck-i{^), defined on the 
generators by 

dX\po,...,Pk\ ■“ (~ 1 ) X[po,---,pi,---,Pk] 

(as usual, summation over i is intended), where pi means that this vertex is omitted. 
With respect to the basis Xs, we write d in coefficients: 

dx^ = dlXi forsGJ^'', 

where summation over i £ is intended. For the whole section, we will sum about 
indices occuring twice in a product, irrespective if they are superscripts oder subscripts. 
Volume terms like |s| or |C/(s)| do not count for this, as s is no sub- or superscript in 
them. The co-boundary operator d* is the dual of d, i. e. a map ^ (.fi) —>■ 
uniquely characterised by d*a{c) = a{dc) for all a £ and all c £ Ck{A)- 

9.2 Remark, (a) By a direct computation, or by common linear algebra knowledge, one 
obtains that the matrix representation of d* is the transposed of the matrix represen¬ 
tation of d. In other words, ( 9 *)* = 5 * for d*f'- = (d*)‘/V 

(b) It would be very natural to write 6^ instead of Xs, because Xs actually is the 
Kronecker delta on . But we will already have some operator 6 acting on differential 
forms, we will define some 6 for cochains and some 6 on piecewise constant forms. In 
the whole following section, we will not use the Kronecker symbol. 

(c) The use of functions as generators of Cfc(h) is only one possible definition. The 
other frequently encountered possibility is to speak of “formal linear combinations” of the s 
themselves (e. e. lHATCHERir200il and lHlRANllboo^ use this definition). Logically, there is no 
difference between both definitions, as the only strict way to define “formal linear combina¬ 
tions” is to use the characteristic functions Xe- However, the existence of both approaches 
introduces an unpleasant notational ambiguity that may disturb a quick reader: Linear maps 
from simplices to R are chains in our notation, whereas they represent cochains in the other. 
Our notation has the advantage to employ s only as sub- and superscript, but not as term, 
which allows for usual summation convention. 
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g. A Piecewise Constant Interpolation of dec 


Lemma. Let A. consist of one single n-simplex. Then the boundary map dk ■ Ck ^ 9-3 

Ck-i, which can be written as (^^ 1 ) ^ -matrix, has rank (Jf). 

Proof. The matrix size just stems from counting the elements in , which arise from 
choosing fc + 1 vertices out of n + 1. 

The rank of dk is proven by induction over k, starting with k = n. Here the statement 
is that dn ■ Cn —>■ C„_i has rank one, which is true because dn ^ 0 . For any k < n, the 
rank-nullity theorem gives that the rank of dk can be computed as dimension of 

its image space minus the dimension of its kernel, which is the rank of dk+i because 
the A:’th homology group of the simplex vanishes. And by assumption, dk+i has rank 
(fe+i) , which gives rank dk = (^+1^) - = (;!), q. e. d. 


Short i ntroduction to discrete exterior calculus (dec). The discrete exterior 9.4 


calculus ( Desbrun et al. 2005, HlRANll[2nn‘J) attempts to build a simple and useable 


finite-dimensional version of the de Rham cohomology based on an intelligent interpre¬ 
tation of simplical cohomology. It calls d* the discrete exterior derivative d, which 
gives that d : C'^(A) —>■ acts as 


d/‘ = dir 


with dl = dl- 


(9.4a) 


If points As € rs for all simplices s and numbers Ofc € R are given, leading to dual 
cells r(*s) as in |4.i6a[ it defines the scalar product of two discrete fc-forms as 

:= OfeOs/Ss^. (9.4b) 

Remark. The numbers Ok do usually not occur in the definition of the scalar product, 
but we will see that they must be chosen as Ofc = (^) to obtain a correspondence to 
piecewise constant forms. The points A® are classically chosen to be the circumcentres 
of the rs. 


The coderivative S is supposed to be dual to d with respect to this scalar product, 
that means {a,dl 3 )fjk = { 5 a, jd) fjk-i for all a G C'^(A) and all /? G Spelling 

out both sides for a = f^ and /3 = /* gives 


l*6l 






«fc l*s| |t| nt 


(9-40) 


Other definitions are obvious: A form a is called harmonic if {Sd dS)a = 0 etc. 


Piecewise Constant Differential Forms 

Situatiou. Let be an oriented regular n-dimensional simplicial complex with a 9.5 
discrete Riemannian metric g, let A. be the corresponding non-oriented complex, and 
let As for any simplex s define subdivision neighbourhoods C/(s). 

Defiuitiou. Situation as in [97^ Let be the space of forms that are constant 9.6 

in U (s) for each s G . 
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Any simplex s G has a volume fc-form dvolr® which can be extended to a constant 
fc-form in whole C/(s). Denoting the extension also as dvolrs, let 


dvolrs 

0 


in [/(s) 
elsewhere, 


and := span{w^: s S 1 ?'=} C 


9.7 Example. Situation as in f^Tg} dimension n = 2 . Consider two triangles rijk and rjil, 
which together contain the subdivision neighbourhood U(ij). Then is the flattened 
unit vector vector field in direction ri — rj in U{ij) and zero elsewhere, w® is the 
characteristic function of U(i), and similarly is the volume form of in rijk and 
the zero 2-form elsewhere. 


9.8 Observation. All basis elements have pointwise unit length with respect to the metric 
induced on the tensor bundles by g, and have distinct support up to null sets, so the 
V? scalar product has diagonal form in the basis w®: 


(asW®,/?s'W®')L2ofc = |C/(s)| as/ 3 s (9.8a) 

Definition. Situation as in lgT^ Let d : — >■ P~^D^ be dehned by 

du}'- = dlu}^, dl := (9.8b) 

|S| 


Let 5 : P —>■ P ^ be defined by 

Soj^ = SI 


\im,i 

|C/(t)|-® 


(9.8c) 


9-9 


Proposition. Situation as in\^ 7 ^ The maps d and 6 fuljill the 


“discrete Stokes’ formula” 



( 9 - 9 a) 


“discrete Green’s formula” {da,l3)^_2Qk = (a, ^/ 3 ) |_2Qfc-i (9-9b) 

for all s G a G , and (3 G In particular, df = 0 . 

Proof, ad primum: For any a G we have f a = Jo:- Now let a = uj^ for 

some t G Then we have 

10.* = a*' 10.* = a* |t| 

drs rt' 


(without summation over t). On the other hand, 





= d Is| 


ad sec.: If one spells out both scalar products with help of |9.8a| for a = cj* G 
and /3 = w® € one gets 5 ® |C/(t)| = |f 7 (s)| for ig.pbi to hold. 


q. e. d. 
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Remark. The discrete Green’s formula [9^^ holds without assumption on the bound¬ 
ary values because the weights | 17 (s)| and |t/(t)| already incorporate the smaller extent 
of S. For a correct treatment of boundary conditions in variational problems, one would 
have to modify [9.8c[ We decided to investigate the original dec setup here. 

Proposition. Situation as in [g 7 ^ The map ik : —>■ i—7> is a cochain 

map, i. e. each square in the following diagram commutes: 









d 

d d 

*0 

ii 




, d d ' 








If flfe = (fc); it is an isometry for each k and a chain map, i. e. each square in the 
following diagram commutes: 


C' 


0 




*0 

pn° 


*1 


Pf 7 ^ 


C" 


pf]" 


Proof. The isometry property is clear by the expressions |9.4b| and |9.8a| for the scalar 
product of and P“^r 2 ^" respectively. The properties dik-i = ikd and Sik = ik-i 5 
only need to be checked for basis elements, so it suffices to show 


1 — 1 

“s]Tf - 


xs 1 — xs 
li[ - Ot 




for &\\ s € Af ,i € A! 


k-l 


The first one is obvious from definitions |9.4a| and 9.8b The second one comes from 
l 9 - 4 bl as 


= 




fe) ak-1 il)^ 


ik^) 


q. e. d. 


Remark. It might seem a little bit queer to use piecewise cons tant forms for this construction 
and not the elementary forms introduced bv lWHlTNEvI (IiqrtI . sec. IV.27) 

^IPo-Pk] ^ PI A ■ ■ • dV ■ • ■ A dA^ 

which would also make i a cochain map. The reason is that we did not succeed to find any 
relation between the scalar prodnct of Whitney’s elementary forms and the dec scalar 
product Ig.qbl This means that although there is a worked-out interpolation estimate for the 
space spanned by a;®, s G bv lDoDZluid (lioTfil l. it gives no possibility to compare solntions 
of variational problems that were computed using the dec scalar product. 

Proposition. Suppose that rAg is a piecewise flat, ['d,h)-small and absolutely well- 
centred realised simplicial complex, that means all circumcentres have barycentric 


9.10 


9.11 


9.12 
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coordinates X\ > a, and that the circumradii are bounded by j 3 h. Then if g is a second 
piecewise flat metric with \{g — g){v,w)\ < it holds for c' := 




- w || l 2 <c'h'^\ujg\L2 
(asa;|,/3s/w|')g| < c'h^g 


Proof. The difference between dg and dg is easiest, because it only involves simplex 
volumes like |s|g and |s|g. These are close to each other by 13.201 The approximation of 
the scalar product involves comparison between the neighbourhood volumes \Ug{s)\g 
and | 17 g(s)|g. These can be estimated if we know how the circumcentres are distorted. 
By |3.12a[ these are controlled by the distortion of the Cayley-Menger matrix inverse 
and inverses of symmetric matrices are treated by 13.211 (which we apply to M_|_ 
instead of g): 

\q^-q^\<ch\\gradgX^\ 


(where r is the circumradius with respect to g) because 4 r^ and = |gradgA*p 
are the corresponding diagonal entries of By assumption, this is smaller than 

q.e.d. 


Connection to the bv Derivative 

Goal. Recall that piecewise constant functions possess distributional derivatives, which 
are (n — l)-dimensional measures concentrated on the jump sets. Their analogue for 
differential forms are the currents from geometric measure theory. (In order to avoid 
“currential derivative” or similar terms, we will speak of BV derivatives.) If our defini¬ 
tion of discrete exterior derivatives is meaningful, it should be connected to this sort 
of derivative. In fact, the BV derivative of a piecewise constant /c-form a also fulfills 
Stokes’ theorem if the jump set is transversal to the integration domain. But as their 
support is (n — l)-dimensional, we will see that its scaling behaviour does match the 
one of full-dimensional {k + l)-forms such as da. 

9.13 Definition. The comass of a fc-covector a is the absolute value of its largest compo¬ 
nent, equivalently: the maximum over all applications of a to simple unit /c-vectors: 

||q;||, = maxQ:(eio A • • • A e^^). 

For completeness, we also define that the mass of a /c-vector is the norm dual to 
the comass: ||u||» = max||Q,||^^2 A differential form a € has locally 

bounded variation (is locally of BV) if 

sup (a, ( 5 / 3 ) is finite for all U CC M, (9.13a) 

/3GC),n'“+bc/) 

Il/3||.<1 

where of course CQil^{U) denotes the space of continuously differentiable fc-forms on M 
with compact support inside U. The space of fc-forms with locally bounded variation 
is called BViocAl^. The globalisation to the space BVfl^ is as usual. 
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Fact (cf. IEvans and Gariepy iqq2 . thm. 5.1). For each a G BVioc^^^, there is a 
Borel-regular measure g, on M and a /r-integrable (fc + l)-form a such that 


(a, ( 5 / 3 ) = J / 3 ) dg for all /3 

M 


G 


(9.14a) 


We will mostly write {dr^a,fi) as abbreviation of the right-hand side. 

Remark. This formulation of the BV structure theorem is the one normally used 
for functions of bounded variation. For differential forms, one calls the supremum in 
ip.i^ai the mass of the current (linear form on 17 ^') /3 1—>■ (a, ( 5 / 3 ), and then observes 
that ever y current of finit e mass is “re presentable by in tegration” in the meaning of the 
theorem ( Federe^Hq sec. 4.1.7, or lMoRGANlbood sec. 4.3B). To obtain uniqueness 
of dF ^one usually requires it to have unit-mass everywhere, and the pointwise scaling 
then comes from g. As we are only interested in dF^a for a G it will be more 

adequate to use a non-unit-length [k + l)-form and the volume form of dU{s), s G 
for g. 

For the proof of 19.141 we refer to IEvans and GariepyI {loc. cit.), because it only 
consists of the observation that /3 1—>■ (a, ( 5 / 3 ) has a norm-preserving continuation to 
Cgf 7 ^, and the application of Riesz’ representation theorem. 

Proposition. For the basis elements tu* o/P“^r 2 ^', the BV derivative is given by g = 
dvola(7(t) d^'^UJ^ = V l\ w*, where v is the outer normal on U{i). 

Proof. If /3 G the product (uj'’FF) is supported only in U(i), where we can 

apply the classical Green’s formula because the integrand is smooth. So 

{d^'^u}\ / 3 ) SP) = j A *P + {dLo\ / 3 ) = y (1/ A w‘, / 3 ) dvolaj7(t), ( 9 -i 5 a) 

dU{i) dU{i) 

the last equality by usual multilinear algebra and diu^ = 0 almost everywhere, q. e. d. 

Proposition. There is a variant of Stokes’ theorem for the BV derivative of P“^ 
forms: If we define 

(9.16a) 


jBV, 


/ := I 

r5ndU{i) 


then 


rs drs 


for a/^aGP^^O^sG#. 


(9.16b) 


Proof. The homotopy formula is easy for constant forms: If da = 0 , then 0 = fjj da = 
Jgjj a, hence f^a = iL a if the integration domains A and P bound a common 
(k + I)-dimensional domain U, the sign depending on the orientation of B. This is the 
case for 

d(rsnU(t)) = drsnU(i') U rsn9C/(t). 

So the formula is clear for a = oj^ hy definition of the “domain integral” over rs, and 
by linearity, it hence holds for all a G g. e. d. 
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B. Main Constructions 


Remark, (a) The notation | 9 .i 6 a| might seem to obscure the actual integration pro¬ 
cess over a subdomain, but we would like to emphasise the analogue to da, which 
fulfills the same Stokes formula. 


(b) For a smoothly bounded (fc-f l)-dimensional integration domain U and differential 
forms a £ /3 € it is always true that 



(a, dvolc/)dvolc/, 


u u 



{v A I3,dvolu)dvolou, 


dU dU 


where the sign is the same as in dvolc/ = A dvolpc/. Therefore, the notation | 9 .r 6 a| 
can also be interpreted as 


J d^'^a = J (d^'^a, dvol„)dvol„nac/(t) 

rs rsndUit) 


9-17 


with d®^w‘ = z/Aw* as in | 9 . 15 [ The notational problem is mainly that the BV derivative 
is supported on a codimension-l-set, which makes the integral in Green’s formula 
(n — l)-dimensional instead of n-dimensional, and the left-hand side integral in Stokes’ 
formula fc-dimensional instead of (fc -I- l)-dimensional. Unfortunately, we do not know 
a common notation covering both. 


(c) The form ula stays correct (with an appropri ate notational adaptio n) for SBV forms 


!EE 

(intro duced bv loE Giorgi and AMBROSiolroSSl as overview we refer to lAMBROSlO et al 


I2000II which are BV forms whose derivative measure p consists of parts fiac and ps which 
are absolutely continuous with respect to the n-dimensional and to the (n — l)-dimen- 
sional Hausdorff measure in rA, if fj,s is supported on a set that is transversal to the 
integration domain U. For the proof, one can use the approximation of a £ SBVU^ 
by convolution with smooth Gaussian kernels. If the jump set of a, i. e. the support of 
fis, is transveral to U, then the convergence is uniform almost everywhere on 9t/, and 
so the integrals J^jj ai of the mollified forms a^ tend to Jgjj a and give a well-defined 
interpretation of dP‘'^a. 

(d) This means that for /? £ which is smooth inside each t7(s), s £ we 

have 


(a, (5/3) = 


S 


a A */3 -I- (d“''a,/3) for all a £ P 


dU{s) 


Proposition. {cP'^a,^) = j^{da,l3) for all a G P and all (3 £ 
Proof. Due tojgTT^ it suffices to consider, for each s £ and each s £ 

{d^^oj\P}ui,) = I A 
c/(s)n9c/(t) 


which can be spelled out by using the (n — l)-flags a in U{s) fl dU{t): 

Juj'-A*P = Y^ j a;‘A*/3= ^^^(w‘A*/3)(6a), 

c/(s)na(7(t) “ A^-i “ 
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g. A Piecewise Constant Interpolation of dec 


where ba is the pull-back of an orthonormal basis of r'a. The flags occuring in this sum 
are of the form ((0),.. ., t,s, ..., (n)), cf. 14.1701 Using the vectors from | 4 .i 8 l 

we have inside each r'a 

~ '^(0),(1) ‘ ‘ ‘ ^ '^5,(fc+2) A • * * A 

= '*^(0).<1> A • • • A A l's,(fe+ 2 ) f\ ■ ■ ■ f\ Vl^n-l),{n) i 

the factors in the latter product are all mutually perpendicular. Now observe 
/ tN# _ ^(0),(1) A • • • A V(^u-l)a , B.u _ Vs,{k+2) A ■ ■ • A Vl^n-l),{n) 

/I I ? J , I ■ 

I'*^( 0 ),< 1 > A • • • A t| l'<^5,<fc+2> A • • • A 

By orthogonality of all vectors in ba, the application {oj^ A *fi){ba), usually comprising 
all permutations of the factors, splits as 

(w* A *P){ba) = u;‘(w(o>,<i> A • • • A Vi^k-i),i) {*P){vs,{k+2) A • • • A 
= 1^(0),( 1 ) A • • • A Vi^k-i),i\ {P,^^)\v5,{k+2) A • • • A 

Summation over all flags ((0),..., t,s,..., (n)) then gives, by | 4 . 19 [ 

/.-A.,3= i.i |..i(,.=.a= 

u{5)nau{i) 

q. e. d. 


Approximation Estimate for P~^r2'" 

Lemma. Let us denote the set of multiindices I = (ii,... ,ik) with ij < ij+i and 
1 < ij < n for all j, by (which in fact is its cardinality). Suppose A. is a simplicial 
complex with only one n-simplex e with a non-degenerate flat metric. Then the x 
(fc+i) raatrix 

rc 

has full rank , where := A arbitrary basis Vj ofTrc. 

Proof. The choice of Vj does not matter, because a change of this basis only results 
in a multiplication with a non-singular x ("^-matrix from the left. Furthermore, 
it suffices to show that the matrix M(^k) •= (/ dvolrt(i' 7))/4 has full rank, because it 
only differs from by factors depending on the volumes 1 17(t) |, which must be 

non-vanishing for at least Q) of the t (which happens if the circumcentre lies on a 
facet of t). Now if we choose Vi = ri — rO, we can transform the situation onto the 
unit simplex D with vertices 0, ei,..., e„, where then Vi = Ci, and the volume forms 
of simplices t containing the vertex 0 have a particularly easy expression: 

dvolrt = dx'^ A • • • A dx'’’ for t = {0, *i,..., ik}- 

So dvolr({o}u/) (^^/') = 1 if / = U and 0 else, hence these rows of M(k) are already 
linearly independent, q. e. d. 


9.18 
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B. Main Constructions 


g.ig Proposition. Suppose rA is a simplicial complex with a piecewise flat and (fl, h)-small 
metric, a € Then there are and, if X\ >0 for all components 

of the As defining subdivision neighbourhoods U{s), there is such that the 

norms of , da^ and Sa^ are estimated by the corresponding norms of a up to a 
constant, and (with the Poincare constant Cg from fgii 1 61) 

{a-a°,P) <Cah\a\\Vl3\ for all fl £ 

{da-da^,P) < Cah\da\\VI3\ for all fl £ 

{Sa-6a'^,P) <Cah\Sa\\Vfl\ for all fl £ 

Proof. Let us introduce a space of “elementwise” P“^ forms, spanned by 


LO 


s,c 


dvolrs in C/(s) n re 
0 elsewhere 


for s C e, e £ 


Then we can find cfl, and such that, for each e £ .S", 

f {5P,P)= f (a,/3) for all constant ,5 £ 

J re J re 

f {dcfl,P)= f {da,fl) for all constant ,5 £ 
J re J re 


(g-iga) 


and similarly for Sa^. In fact, writing and inserting [3 = v\ from above 

with / £ , the equations for determining coefficients d® ^ have as system 

matrix, which has full rank by | 9 .i 8 [ The integral J{da^,fl) reduces to a boundary 
integral by 19.171 which does not include d anymore, this boundary integral is invariant 
under affine transformations, and the problem is solvable on the unit simplex D. If all 
A* are positive, the are just a row- and column-rescaling of df by non-zero factors. 

As a second step, let £ P~^n*^ be defined by the condition that, for each 

t £ A'" and all s £ 


[ {cP,P) = 

f {5P,P) 

for all constant fl' £VP, 

Ju{i) 

luit) 


[ {da^J) = 

[ {ddd,fl) 

for all constant fl' £ 

lu{s) 

fuC) 



(and similarly for Sa^), which means averaging over all parts U{t) fl re with t C e. 
These forms have the desired properties: Norm-preservation is clear by construction. 
For the approximation, observe that fl £ Xlf can be replaced by some fl that is con¬ 
stant in each element re, and the error is estimated by the Poincare inequality 12 .locl 
1/3 ~ i5h^|V/3|l 2. And similarly it can be replaced by some fl' that is constant 

in each neighbourhood U (i). This gives 

{a — a°,fl) = {a — cP ,fl) + {oP — a°,fl} 

<{a- aPfl) + (d° - aPfl') + Cah\a\ |V/3|, 
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g. A Piecewise Constant Interpolation of dec 


and the scalar products vanish by choice of and a^. Exactly the same computation 
is feasible for {da — da^, f3) and {6a — 6a ^, /3), q. e. d. 

Proposition. For a complex consisting of only one n-simplex e, let d(fc)(e) be the 
(fe+ 2 ) ^ (fc+i) representation of d : —>■ 

Suppose A. is a simplicial complex with a piecewise fiat and {D, h)-small metric. 
Assume that the x -matrix 


VM(fe+i)(e)d(fc)(e)y 


has full rank for each e € ^ . 


( 9 . 20 a) 


Let Cg he the Poincare constant from [ 2 .iih\ Then for each a G , there is a G P 
with |d|Lp ^ |q!|lp; l^dlLP ^ IdaiLP and 


{a-a,P) <C^^h\a\\VP\ for all/I G 
{da-da, 13) < C^h\da\\VI3\ for all /3 e 


( 9 . 20 b) 


Proof. The assumption |().:aOa] guarantees that we can find one single a G 

filling both equations in iq.rpalat the same time, q. e. d. 


Remark, (a) We did not succeed to verify [ 9 T 2 oa| in the general case, but there are at 
least no structural obstructions for it to hold: M(fc) has full rank (^), and d(fc), which 
is the transposed of 9fc+i from [ 9 ^ with rows scaled by |t| and columns scaled by |s|, 
has rank (j."]^), which add up to 

(b) Using an Poincare inequality (|Evans and Gariepy iqq2 . thm. 4 . 5 . 2 ) instead 
of i 2 .ioci leads to {a — a^,f3) < C^^ph\a\\_p\V f3\\j> and similar for do —do ^ andda — 

(c) The estimates | 9 . 2 ob| are formulated as norm estimates. By inserting mollified 

characteristic forms or forms with small support and |//3| = 1 (“Dirac forms”), one 
can localise the convergence. 


(d) We do not call |Q. If)] and [^720] “interpolation estimates”, as a°, a^, o? or a may 
have nothing to do with a pointwise, but only in integral mean. In contrast to interpo¬ 
lation, the integral of a and over smaller spaces like boundaries of the U (t) will in 
general not converge, as 19.171 shows. For an example, see also 110 . 30 ] The section title 
“interpolation of dec” does not refer to interpolation of smooth functions, but to the 
process of extending the simplicial definitions in to forms in | 9 . 11 [ 


9.20 


9.21 
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C. Applications 


10. Real-Valued Variational Problems 


Situation. Using the results of the preceding sections, we do not speak of a manifold lo.i 
and its triangulation, but directly suppose that M = rR is a realised n-dimensional 
regular simplicial complex (compact, as usual), endowed with a piecewise flat and 
(d,/i)-small Riemannian metric g®, as well as with a smooth metric g fulfilling I6.17I 
and |6.23| with C'^^h < 1 . Except for lro.iST jga.. we assume that if M has a boundary, it 
follows the boundary of the Karcher simplices. Therefore, the homeomorphism property 
of X remains unchanged. 

Remark, (a) Obviously, the spaces C^’“ of strongly differentiable and Lipschitz func- 10.2 
tions for g and g® (defined in the classical meaning for g and byfor g®) are different, 
but as the Sobolev norms for differentiation orders fc = 0 , 1 , 2 are equivalent, the spaces 
W^’’'(Mg) and W^’’'(Mg®) coincide. 

(b) T he convergence of curve length and geodesic distance, treated ui IHildebrandt et al. 

1 2006L sec. 4.1) for the case of embedded surfaces, is already covered by I6.13I in our 
setting. 


The Dirichlet Problem for Functions 


Goal. In this section, we will deal with approximations of the Dirichlet problem, that 
is solving a weak version of Ait = /, where A is the Laplace-Beltrami operator of M. 
The Laplacian of ^-differential forms will be dealt in the subsequent section. 

As first step, we will give a short review of the usual proof for converg ence of Galerki n 
approximations to the Dirichet problem as can be found for instance in iBRAESsI (I2002I). 
In the second step, we will add the usual error terms resulting from the “variational 
crime” to u se q® in s tead o f g. This is standard in the EE theory for geometric PDe’s 
initiated bv IPziui^ ( loSS h but often not separated from the error of the first step. 


Definition. Situation as in 110.11 Denote by Hg the space of weakly differentiable 10.3 
functions with vanishing trace on dM and by the space of globally continuous, 
piecewise linear functions (here, “linear” of course means usual linearity in the parame¬ 
ter domain re, e € by Pg the same but with vanishing boundary values. The norm 
of an operator on function spaces will be denoted by || - ||. 


Definition. For v,w G P^, recall the definition Lap(M,u) {du,dv) from 10.4 

and that the (homogeneous) weak g-Dirichlet problem is the task to find u G Hg(M) 
such that Lap(u,?;) = {f,v)j^g for all v G Hj, we shortly write Lu = / with an 
operator L : Hg — (Hj)*. The 5-Galerkin solution to the Dirichlet problem with 
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C. Applications 


10.5 


10.6 


10.7 


10.8 


10.9 


10.10 


respect to the trial space Pj is the solution Uh to La.p{uh,v) = (/, v) for all € Pj 
Naturally, there is also the notion of a g^-Galerkin solution. 


Remark. By 1 : 2 . 1:51 we know that the Dirichlet problem has no solution for general 
/ G L^, but only for / _L ii, and the solution is unique up to harmonic components, 
in other words: there is a unique solution in But the space of harmonic functions 
is one-dimensional, consisting only of the constant functions—and these are ruled out 
by the boundary value requirements. 


Fact ( Schwa^Iiqq.4 also cf. |2.19| ). The de Rham complex (H^’°n,d) of a smooth 
compact Riemannian manifold is a Fredholm complex, so the Dirichlet problem is 
uniquely solvable, and |dM|L2 < GqI/Ils with the Poincare constant Cg from l^obl 
This means that L~^ is a bounded linear operator. 

If dM is piecewise smooth or convex (that means, convex where it is not smooth), 
then M is H^-regular, i. e. there is a constant Cq depending on M, but not on /, 
with |rt|H2 < Cq I/Il^j that means that ||L“^||l2 h2 < Cq in this case. 


Lemma (Cea). Situation as in \in.n Let u be the Dirichlet potential and Uh be the 
g-Galerkin solution to f € L^. Then Uh is the orthogonal projection of u onto Pq with 
respect to Lap( •, • )• 

Proof. As Pg C Hj, also u fulfills Lap(u,u) = {f,v)^_2 for all u € Pj by which Uh was 
defined. So we have the so-called “Galerkin orthogonality” Lap(u — Uh,v) = 0 for all 
such V G Pj, which is the characterising property of the projection error, q. e. d. 


Corollary. Situation as in \io.ii Let 11 be the orthogonal projection Hj —^ Pj with 
respect to Lap( •, •) and B'*- := id—11 he the projection error. Then for any k for 
which both sides are defined, 

\u-Uh\H>^ < \ f\i 2 ||n*^L"^||L2,H'=- 

Proposition. Situation as in \iu.i \ with dimension n < 3 . Then ||n^||n2_Hi ^ 'd~^h, 
and if additionally M is -regular, then 

\u-Uh\H^ < C'0W^|/|l2. 


Proof. It suffices to show that there is one Uh € P^ with \u — Uh\\p ^ 'd~^h\u\^2, then 
the projection of u will produce a smaller error than this Uh. As usual, we take ut to be 
the Lagrange interpolation of u (which is well-defined, as C C° in dimension < 3, cf. 
Adams! lovtj . Theorem 5.4.C). And this interpolation estimate is exactly [77^ q. e. d. 


Proposition (Aubin—Nitsche). Situation as in \in.ii Then ||n'*-L < 

Under the same conditions as in \io.g\ 
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10. Real- Valued Variational Problems 


Proof. First, note that for a right-hand side g, the solution L is characterised by 
{g,v)Mg = Lap(L“^5,n) for all v € Hj. Now for a right-hand side / S L^, consider 

= sup ^ 1^9)1^ 

geL2 |g|L2 

Lap(n-^L“V:^“^ff) (*) Lap(n^L~V, 

= sup-1-|- = sup-1-|- 

\ 9 \l^ \ 9 \l^ 


where we have used in (+) that 11 and hence If-^ is a Lap-orthogonal projection, q. e. d. 


Remark. It would of course be possible to consider other interpolation procedures 
th an just nodal Lagrange interp olation, for example averaged Taylor polynomials as 


Brenner and Scott! ( 2002L section 4.1), which would circumvent the dimension 


restrictions. However, the emphasis of this thesis lies more on the different possible ap¬ 
plications of the Karcher simplex construction than on optimal results for the Dirichlet 
problem. 


Lemma. Situation as in \io.ii Let F{v) := {v,f)Mg’ Lap*^ and he defined 

similar to Lap and F, but with g^ instead ofg everywhere. Then |(Lap — Lap®) (n, w)| 
CQh'^\dv\\_2\dw\\_2 and\{F — F^)v\ C'^hf \v\\_2. 

Proof. Exactly as in [7^ q. e. d. 

Remark. In the understanding of lHiLDEBRANDT et al ] (l2006l l. the “weak Laplacian” Lg is a 
mapping —>■ (H^)*, LgU : v 1 —> Lap(u, v). In this setting. rid. 121 can be seen as a convergence 
result for the weak Laplacians: \Lg — Lgs |||.|i < Cgh^. 

Proposition. Situation as in \iu.i\ with -regular M. Letuh,uf^ € Pj be the Galerkin 

solutions to LgU = F and Lgeu*^ = F^. Then 

\uh — ul\\_2 -\- Cq \duh — du^^\\_2 < C'QC^h^\f\\_2. 


Proof. During this proof, | • | always means | • h^^Mg)- Let us first consider the deriva¬ 
tive term on the left-hand side: For some v with |i;| = 1, we have 


\duh - dull = Lap{uh - ul,v) 

< \La.-p{uh,v) - Lap®(u^,n)| -b |Lap®(u^,n) - Lap(u^,'i;)| 

< \{F- F^)v\ + |(Lap®-Lap)«,'(;)| 
<C'oh^\f\\v\+C'oh^\dul\\dv\ 


Then use \dul\ < Cq |/| from I in. HI For the estimate of \uh — ul\, use the Poincare 
inequality again, q. e. d. 

Remark. As in the euclidean setting, the proofs carry over to an arbitrary continuous, 
strongly Hg-elliptic bilinear form on instead of Lap. 


10.11 


10.12 


10.13 
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C. Applications 


Variational Problems in 

10.14 Assumption. Situation as in 110.11 For k = 0 ,...,n, let there be finite-dimensional 

subspaces of (or if needed) with V? and approximation order 

h analogous to pfT^ 

min lu — Vh\\_^ + \dv — dvh\\_2 + |( 5 u — Svh\\_^ < ah|u|H2 

DhePfi'" '- V -^ 

only for |io.i5^ 

and similar for t*v = t*Vh = 0 or nu = nvu = 0 . Furthermore, assume that the 
Dirichlet problem is H^-regular and the Hodge decomposition u = da + 6b + c is 
H^-regular, which means |c?a|Hi ^ |u|hi etc. We abbreviate (•, • )\ji(Mg^) as (•, • )e. 

10.15 Proposition. Assume \iu.i4\ Let u = da + 6b -\- c be the Hodge decomposition of 

u G which can be computed as a = argminP[u] over a S and b = 

argminG)?/] over b £ as in \2.i4\ If ah = argminP[u] over ah £ and 

bh = argmin G[ti] over bh £ Pf 2 ^+^, then 

|da — dah\\_^ + \6b — 6bh\\_2 < ahlulHi. (10.15a) 


If u = dae + 6be -P Ce is the Hodge decomposition with respect to g^, and if ah,e o-nd 
bh^e «re defined similiarly, then 

\da — dOelL^ + — 6be\\_'^ + |c — Ce|L2 < C'^hf |m|l 2, (10.15b) 

\dah - da?i.e|L2 + \ 6 bh - 6 bh,e\L'^ < Gg/i^ |m|l 2. (10.15c) 


Proof, ad primum: By the Euler-Lagrange equation (da, dv) = {u, dv) for all v £ 
and {dah,dv) = {u,dv) for all v £ PH^“'"^, we know that dah is the L^-best 
approximation of da in (i(PHj), which is smaller than Q!h|Vda| by assumption. 

ad sec.: If {dah,e,dv)e — {u,dv)e = 0 , then {dah,e,dv) — {u,dv) < Gq/i^(| da^^e| + 
|m|)|u| and hence {dah — dah^e,dv) ^ C'o^^( Ma/i,e I -P |m|) |i;| for all u G PHj. The same 
calculation is valid for doh — dah,e instead of da — doe- The c — Ce estimate comes out 
as the remainder, q. e. d. 


10.16 Remark, (a) |io.i5b| is our analogue of thm. 3.4.6 in I WardetzkyI (l2on6h . 

(b) In general, there will be no exact finite-dimensional Hodge decomposition in Pfl^', 

as we have not required any connection between d(PH^) and There i s a Hodge 

decom position in the space of Whitney forms with convergence proven by IPODZIUK 
(|iq 76L thm. 4.9). Variational problems in a specific space P”^H^ of piecewise constant 
forms will be treated in |io.24^ qq. 

(c) The FEEC setting of Arnold et al. only has a weak Hodge decomposition u = 
dah + bh + Ch oi u £ PH^ as in 12.131 but as its parts are also orthogonal projections, 
there is an estimate 


|da/i dah.eli.^ T \bh T \^h Ch,e\\_'^ ^ |'a|[_2 

corresponding to |io.i5b[ 
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10. Real- Valued Variational Problems 


Mixed form of Dirichlet problem. Arnold et all l|2on(iL bninl) have shown how 


to construct finite-dimensi onal subcomplexes (Pfl, d) of d) and solve the mixed 

3 


Dirichlet problem therein. IHolst and SternI ^|2012^ have extended this to the situa¬ 
tion where the domain of the Sobolev space and the finite-dimensional approximation 
are endowed with different, but close Riemannian metrics g and g®, which leads to the 
situation that the inclusion map \P’^Vt^{Mg) is not norm-preserving 

anymore, but only an ahnost-isometric map. Their setting directly applies to Finite 
Element computations on the Karcher-Delaunay triangulation: 

Proposition ( Holst and Stern 2ni2L thm. 3.10). Assume \10.14[ and use the no¬ 
tation from For f G let {a,u,p) G P© and {ae,Ue,Pe) G P©e be the 

solution of the mixed formulation \2.1M of the Dirichlet problem in Mg and Mg^ re- 
spectively, whe r e PG = PH^ x PH^x PJOj' is a stable choice of trial spaces from 
IArnold et al. 1 20od . eqn 7.14), and P&e differs from PG only by the last factor 
(Pi3g)t, the harmonic trial functions with respect to 5®. Then 


C 

Icr-CTelHi + Im-MeIhi + ^ 


where 7 is the inf-sup constant as in \2.i8\ (but over ). 

Proof. The solution s = {a,u,p) with respect to the “correct” scalar product g fulfills 
b{s,t) = F{t) for every test triple t = {T,v,q) G PG. On the other hand, the distorted 
solution Se fulfills be{se,te) = Fe{te) for all tg G PGe with the obvious definition of be 
and Fg. As the trial spaces only differ in the last term g, we have 


bei^ey ^e) — bgi^Sg^ -f- (tie, qe ^)e — bg^Sg^ -\- (Ug, q^ g 

(because Ug T P^g)- Now observe {ug,q) = {TTUg,q), where tt is the orthogonal pro¬ 
jection onto Pij^, and by lio.iG^ the projection of a Pijg element onto Pij^' is small. 
Hence 

{ug,q)g = {ug,q) -P {{ug,q)g - {ug,q)) < \u\\q\. 

Weakening the right-hand side, we obtain \bg{sg, t) — Fg{t)\ < C(h^ |s| Kl- By the scalar 
product comparison [TnHH also \b{sg,t) — F{t)\ < C(h‘^ |s| |t|, and taking this together 
with b{s,t) = F{t), we have 

b{s- Se,t) < C(h'^\s\\t\. 

Now, by the inf-sup-condition b.iSbl yjs — Se| < sup^ b{s — Sg,t) / \t\, q. e. d. 

Dirichlet Problems with Curved Boundary 

The case that the analytical and the computational domain actually coincide is not the 
only interesting problem. When for example a Dirichlet problem on the unit disk in 
hyperbolic space is considered, a Karcher triangulation with respect to the whole hyper¬ 
bolic space will not exactly cover the unit disk. But the treatment of such a boundary 
approximation is standard in Finite Element theory, and the main task is to carefully 


10.17 
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C. Applications 


inspect which arguments have to be modified because th ey rely on the Euclidean struc - 
ture of the domain. We give a presentation according to IDorfler and RumpfI (|iQQ8li 
and do not treat the difference between g and (7®, as this comparison can be done 
separately by using |10.13| after llo.Sll 

The usual setup for boundary approximation is that a domain 17 is replaced by 
a simplicial domain 17 ^ whose boundary vertices lie on 917 . By (n — l)-dimensional 
interpolation estimates, one then gets that 917 and dflh are only < far apart, 
where k bounds the curvature of dfl and h the mesh size of flh- We translate this, for 
17 C M, into the following 


10.18 Situation. Let M = rA. he a, piecewise flat and (d,h)-small realised simplicial com¬ 
plex. Let 17 C M be a full-dimensional domain and Q.h = rA a realised full-dimensional 
subcomplex, connected by a “normal graph map” $ : 917 —>■ P exp^ d v, where 
d : 917 —^ ]R is Lipschitz-continuous and is the outer normal on 917 , with the follow¬ 
ing properties: First, the retraction inverse {p,t) 1—>■ exp^idz/ is injective (to ensure 
that no topology change may happen). Seond, it is “short” in the send that |d| < ah^, 
Lip d < ah < 1, and |Vdd| < a (where d is smooth) for some a € R. Let all principal 
curvatures of 917 in M be bounded by n. This implies that for small h the norms of 
and Vdd) are bounded, see section rm 


10.19 


Lemma. Situation as in \io.i 8 i If v € H^( 17 /j), then |'y|L2(a;,\a) ^ cth^\dv\\_2(^Q^\Q'f 
and |n|L2(annnfe) ^ for small h. 


Proof, ad primum: It suffices to show the claim for smooth v. Consider A € Hh \ 17 . As 
d(A, dfth) ^ ah'^, there is a curve 7[A] : ^ ^ A for some p € dflh with length < ah"^. If 
h is small, this curve can be supposed to be a straight line lying entirely in one simplex 
of A. As v{p) = 0 , 



(10.19a) 


7 [A] 


Suppose 7 is arclength-parametrised. Now we can again apply the arguments from the 
proof of [ 7 ^ (keeping in mind that 7[A] has length h there, but ah^ here): 


r JlO.lfial r / r \2 17.5q1 r r 

j \v\^ < J (y Mnl) < ah'^ J y Idnp 

n/iPiO 7 [a] 7[a] 




ad sec.: Because 917 ?, is a graph over 917 , the inverse is also true: 917 is a graph (usually 
not normal) over 917 /j, so we can introduce coordinates in which a simplex f of dllh 
lies in the Sm-plane and 917 is parametrised by (a:i,... {xi,... ,Xra+iTp). 

Then 


lL2(annrf) 


r . _ ni.lObi 

y kiVi + MH" ^ 


illO.lQal 

< 


rf 


f / f \2 u .oai 

J [ J \M) ^ o:h^\dv\l2^n^\Q), 


rf 7 [A] 


q. e. d. 
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10. Real- Valued Variational Problems 


Lemma. Situation as in \in.itA For u : M —^ R, which is continuous in fl and 
M \ n, let [ti] be the jump of v across dQ. If h is small, there is a continuous ex¬ 
tension u of u € H^(r2) onto flUflh such that u\q = u, |u|h2(o^\q) |u|H2(n) o.nd 

IM“HlL2(annnh) ^ |M|Hi(n)- 

Proof. By assumption, all points in Hh \ are covered by the homotopy 


expp tv, 

where at each p G dll fl Hh, the parameter t is chosen within ]0;d(p)] (in particular, 
points with negative d(p) are excluded, as they would parametrise Q \ Qh instead of 
flh \ il). For an image point of $t, set M(expp tv) := u(expp —tv), the reflection along 
dH. This u is continuous, and [du v] = ± 2 c?m v. The norm-preservation follows from 
the assumptions on $ (but note that u is not in U fl due to the jump on dfl, 
even though is smooth), q. e. d. 


ProDosition. Situation as in \io.i 8 i Let u G H2(f]) be the solution of Lu = / with 
respect to fl, and let Uh G Pg(n;j) be the Galerkin solution over flh for an extension 
of the right-hand side f by zero onto \ fl. Then \du — duh\i2(^ci) ^ y/ah\u\y\2(^Q;^ for 
small h, where Uh has been extended by zero in Il\ Llh- 

Proof. Let u be the extension of u from nUT^ Assume we can show 


Iw-mIhho;.) ^ |m-' y|Hi(nh)+ah2|u|H2(n,,\o)-fv/ah|[du(j2)]|L2(aQnn^) (10.21a) 

for every v G P^(fl/j). Then the claim is proven by [775] and 110.201 Supposed v G P^, 
observe that in \dv — duh\ = sup {dv — duh, dw )/ |dw|, it suffices to take w G P^. So 
we have 


\du — du/j|L2(Q^) < |c?it — dv\ \dv — duh\ = \du — dv\ -f sup dw) 


= \du — dv\ -\- sup 
< 2 1 du — du I -I- sup 


{dv — du, dw) -\- {du — duh, dw) 
|du;| 

(du — duh, dw) 


|du;| 


And now, if / is the extension of / by / = 0 in fl/j \ fl, 
(du — duh, du’)L2(Q^) = J{du, dw) — fw -P J{du, dw) 


- fw 


J {—Au — f)w-\- J wduv + J—wAu-\- J wdu{—v). 


n^un 


=0 




r2h,\r2 


0(n;,\r2) 


as —V is the outer normal of fl/j \ fl. So this gives 

(du - du;i,du;)L2(n^) < l^w|L2(nfe\n) |w|L2(n;.\n) + |[dM(i^)]|L2(annnfe) |w^lL2(annn^): 

which shows, together with 110.191 for the w norms, the claimed estimate iio.2iaL 

q. e. d. 


10.20 


10.21 
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C. Applications 


Heat Flow 

Goal. As a short outlook on Galerkin methods for parabolic problems, we consider 
the approximation of heat flow under pe rturbations of m etric. We decided to exclude 
the general convergence theory (see e. g. Thome j 20061 chap, i) and concentrate on 
the difference between Galerkin approximations with respect to g and g®. 

Proposition. Situation as in \iu.n For a time interval [ 0 ;a], let Uh,Uh,e be the time- 
eontinuous Galerkin approximation to the heat flow with initial value Uq € Pj and 
right-hand side f € L°°([ 0 ; a], L^(Afg)) for metrics g and g^ respectively, that means 


{uh,v) + {duh,dv) = (/, v) 
{Uh,e,v)^ + {dUh,e,dv)^ = {f,v)^ 


for all u £ P^, 
for all u £ P^, 


Uh\t=o = Uq, 

^h,e\t—0 — ^Oj 


where (•, •) ^ is the abbreviation for (■, ■) ■ Then their differenee ean he estimated 

by 

|m/i — M/i,e|L“(L2) ^ |mo|h 1 + |/(t)|L2^ 

Pro of. The proof foll ows the line of the usual convergence proof for parabolic problems 
as in lTHOM^ l 20o6l thm. 1.2): Gonsider e := Uh — Uh,e- By the defining equations for 
Uh and Uh.e, we have 


{i,v) + {de,dv) = {f,v) - {uh,e,v) - {duh,e,dv). 


By [ 77 g] and 110.121 we have 

\{uh,e,v) - {dUh,e,dv) - {f,v)\ < C'^ih'^ {\uh,e\\v\ + \dUh,e\\dv\ + \f\ |u|), 

where all norms are norms. So we have for v = e, together with the Poincare 
constant Cq from [2TTocl 


1 "S 1^1^ + - C'oC' 0 ^^(|uh,e| + \duh,e\ + I/I) \de\. 

Then Young’s inequality gives 2 cab < c^of + 6^, hence we obtain a separated summand 
\deY on the right-hand side, which can be cancelled (the suppressed constant belongs 
to c): 

Integration over [ 0 ;a] gives, as £:|t=o = 0 , 


|ep < (GqGq/i^)^ ^ |u?i,eP + \duh^, 


I/P 


200^ eqn. 1.20, 


From the usual regularity theory for parabolic problems (|ThomeeI 
case m = 0), we know that /(|up + \u PhO< I uoIhi + / 1 ,/P, which shows the desired 
estimate for e, q. e. d. 
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10. Real- Valued Variational Problems 


Proposition. Situation as above. Let he the Galerkin approximation to the 

heat flow with implicit Euler time discretisation with respect to g and 5® respectively, 
that means 

{dul,v) + {dul,dv) = {f,v) for all v&Pl, = uq, 

+{dul^^,dv)^ = {f,v)^ for alive Pi, = uq 

for the backward difference quotient Bv^ := -. Then their difference at time 
t = nr can be estimated by \duh — dufl ^\i2 < Khft, where K depends on the geometry, 

I/Il~(l2 ) and |mo|hi- 

Proof. As before, let e" '■= Then 

{de'^,v) + {de,dv) = {f,v) - {dul^^,v) - {duh,e,dv), 
and the right-hand side is bounded by 

\{dK^e^v) - -h \(duh,e,dv) - (duh,e,dv)J -P\(f,v) - (f,v)J 

< Clh^{\dulJ\v\ + \duh,e\\dv\ + |/||u|) 

<ClC^hmdul,\+\dUH,e\+\f\)\dv\. 

Denote the whole term in parentheses as A. As before, it is bounded in terms of the 
given data. Then again the choice t; = e" gives 

|e"P - (e^-\e^) + |de’"p < C^CQ/rVAlde^p 

and so 

|e”p-h Me”P < C'^CQ/iVA|(ie"p-HC'2 |de”-^||de"|. 

And of course Ide^p is smaller than the last left-hand side, which gives |c?e”| < 
Cq j^Cq/i^tA -|- Cg Then the claim follows by induction over n, q. e. d. 


Discrete Exterior Calculus 


Observation. As we have noticed in 12.191 variational problems from section [2] 

/- 7 '\ Vw» r cn+-T'T 1 /-X-*-* /-x-p D 1 


are uniquely solvable in like in by the construction of as a 

(co-)chain complex. As just a gentle way of writing the simplicial cochain 

complex {C^,d*), its (co-)homology is isomorphic to the de R ham complex’ on e (a 


short direct pro of, called “the th eorem of de Rham”, is given in IWhitneai iq.i^tL sec. 


IV.29, although Ide R.HaI 3 IiquI proved isomorphy to singular, not simplicial coho¬ 
mology). Therefore, we can hope for approximating smooth solutions of variational 
problems by ones in 


Situation. Let rA be a realised oriented regular n-dimensional simplicial complex 
without boundary with a piecewise flat, (1?,/i)-small metric g. Let Aj, s £ A*, be the 
simplices’ circumcentres, and suppose A* > 0 for all their components (i. e. r^g is 
well-centred). Assume |9.2oa| and that the Hodge decomposition u = da 6 b c is 
H^-regular, meaning |(ia|Hi |u|hi etc. We use the Poincare inequality in the form 
b.iibl 


10.23 


10.24 


10.25 
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C. Applications 


10.26 Proposition. Situation as in \io.2i^\ For a function f € let u € be the solution 
of the Poisson problem {du,dv) = (/, w) for all v G and let ut G be the 

solution of {duh,dvh) = {f^Vh) for all Vh G P~^r 2 °. Then 

(du - duh.dvh) < Cq/i(|V/| \vh\ + \ydu\\dvh\) for all Vh G 

Proof. Let v and Vh be connected by |9.2ob[ Then 

{du — duh,dvh) = {du, dv) — {duh,dvh) + {du,dvh — dv) — {f,v — Vh) + {du, dvh — dv), 


and both terms can be estimated as claimed, q. e. d. 

10.27 Proposition. Situation as m p0.251 u = da + 5 h-\-c be the Hodge decomposition of 
u G and let u = dah+Sbh + Ch be the Hodge decomposition of its -orthogonal 

projection onto Then 


{da - dah,dvh) < C^h\u\y^i \dvh\[? 

{Sb - 6bh, 6vh) < Cq/iIuIhi 


for all Vh G 


Proof. We know that da is characterised by {da, dv) = {u, dv) for all r; € H 
Naturally, ah is characterised by {dah,dvh) = {u,dvh) for all Vh G P~^n^, but the 
right-hand side is {u,dvh) if du is the orthogonal projection onto P~^. So we can 
proceed exactly like before, but using [ qTT^ to connect only dv and dvh instead of v and 
Vh- 


{da — dah,dvh) = {da,dv) — {dah,dvh) + {da,dvh—dv) = {u,dv — dvh) + {da,dvh — dv), 

the |Vda| produced by the latter term can be estimated by |u|hi by assumption. The 
same procedure is feasible for 5 b and 6 bh (where another test form v can be employed 
such that dv is close to Svh), q. e. d. 

10.28 Proposition. Define 6^ := x x and P“^© x x 

Suppose s = {a,u,p) G is a solution of the Poisson problem in mixed 
form as in \2.ij\ and Sh = {cfh,Uh,Ph) G P^^© is the solution of the corresponding 
finite-dimensional problem. Then for all th = {Th,Vh,qh) G P^^©, 

b{s — Sh,th) |V/Il2 + |Vs|l2) 

where the left-hand side is of course not to be taken literally as in \2.iH(A but with P“^ 
exterior derivatives for Sh and th, i. e. consisting of terms like {du — duh,dvh) etc. 

Proof. In the spirit of lio.261 we start with 

b{s - Sh, th) = b{s, f) - b{sh, th) + b{s, t-th). 

As before, the first two terms are {f,v — Vh), which is well-controlled by the right-hand 
side of the claim. In b{s, t — th), there are many easy terms, which we do not explicitely 
discuss once more. Only {u, q — qh) = {u, qh) is iteresting. Estimating it actually 
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11. Approximation of Submanifolds 


means bounding the difference between and Let u = da + Sb he the 

Hodge decomposition of u (which does not contain a harmonic term), and choose doh, 

5 bh close to them in the sense of |9.1()[ Then, as qt G 

{u,qh) = {da + 5 b,qh) = {da - dah,qh) + { 5 b - 5 bh,qh) < Ci^h{\Vda\ + \ 57 db\)\qh\, 

q. e. d. 

Remark, (a) In |10.a7[ we cannot say anything about c — Ch (yet), as we would need 10.29 
to control the terms in (c — Ch, Vh) = {dau — da, Vh) + { 5 bh — 5 b, Vh), but we only have 
control over the scalar product with dvh or 5 vh respectively. 

(b) As remarked in |9.io[ the correct treatment of variational boundary value problems 
in P~^H^ would require a modification of 5 at boundary simplices. 

(c) Employing P~^ forms also as test functions is unsatisfactory, as they are no clas¬ 
sical objects to test with, so the results are not easily comparable to usual estimates. 

However, we are not sure which forms would be the right ones to test with: Perhaps 
forms that are “almost constant” in some small, but non-shrinking region would be good 
to obtain an average value for such a region. H^d forms are not the right candidates: 

If for example in the Poisson problem {du — duh, dv) converged for all “well-behaving” 
u € H^, then also \du\^ — would converge, which is not the case: 

Example. Consider an equilateral triangle mesh in the xy plane with unit edge length, 10.30 
rotated such that one of the edges is parallel to the a;-axis. Now consider the constant 
vector field v = (1,0), and let us compute its P^^H^ approximation according to |9.:ao[ 

In a triangle ijk, where ij points in x direction, we have Jj^ji.{vh, ( 1 , 0 )) = ^ \ijk\ (a^ -1- 
f Oife + f Ofcj)- So aij = 2 and aik = akj = 1 gives the correct integral mean (the scalar 
product with ( 0 , 1 ) is obvious). There are other combinations matching (for 

example the “obvious” choice aij = 3 and aik = ctkj = 0 ), but this one also captures 
its exterior derivative: As v is constant, dv = 0 , and dvh = 0 for Vh = -|- +uj^h 

This vector field v is also the gradient of the function / : {x,y) 1—>■ x. If fh € P^^H^ 
has the same values as / on all vertices, then = Vh- The homogeneous Dirichlet 
problem in is equivalent to the Dirichlet problem in P^ by |3.2C[ for which reason 

fh will be the P^^ harmonic function with prescribed boundary values fhldn^ no matter 
where the boundary is drawn. Thus, / and fh with derivatives v and Vh are the solutions 
compared in lio.26l As |u| = 1 everywhere, we have = \ijk\ in each triangle, 

but \dfh\l2(^ijk) = I \ijk\ ( 4 -tl-Pl) = 2 \ijk\. The discrepancy \dfh\l2r,jk) = 2|rf/lL2(yfe) 
does not shrink with smaller and smaller edge lengths, so the Dirichlet energy of the 
P~^ approximations will always be twice as large as of the analytical solution. 

11. Approximation of Submanifolds 

Extrinsic and Intrinsic Karcher Triangulation 

Definition. Let the piecewise smooth n-dimensional submanifold S' C M be given by 11.1 
a bijective triangulation y : rA —> S with vertices pi = y{ri). If j/j^c for each e G ^ 


77 














C. Applications 


is a barycentric mapping with respect to the induced metric g\s, then y is called an 
intrinsic Karcher triangulation. If moreover each j/j^c is also a barycentric mapping 
with respect to the metric of M, then y is called an extrinsic Karcher triangulation. 


Goal. The possibility for the existence of an intrinsic Karcher triangulation has been 
dealt with in section |51 The question of this section will now be how well an extrin¬ 
sic Karcher triangulation, induced by the same complex A. and the same vertex set 
{pi}, approximates S. Note that such an extrinsic Karcher triangulation is always an 
interpolation of the given triangulation of S in the sense of |7.14[ 

11.2 Proposition. Let y be a piecewise (tt, h)-small intrinsic Karcher triangulation y of 
S with \Wi,\ < K for all Weingarten maps Wi,. Suppose that all vertices pi = y{ri), 
i £ t, lie in a common convex ball with respect to g for each e G .ft". Then for small 
edge lengths £ij := ds{pi,Pj) and iij := d{pi,pj) with respect to g\s and g, it holds 


Proof. There exists an extrinsic Karcher triangulation x of some set S' C M with the 
same combinatorics and vertices as y (that means: interpolating y) because the vertices 
of each simplex are contained in a convex ball. We do not know if S' is a manifold, 
because the fullness of the extrinsic simplices is not clear a priori, but will be a result 
of the length estimate. 

Let us show the claim for the edge 7 : Cj in re with tangent v := 7. The 

estimate [7:^ can be extended to the whole edge: 


\(dx 


Pdy)v\ '^h'd ^ 


\{\ 7 dx 


- P\ 7 dy){v,v)\ 
\dyv\ 


As edges are mapped to geodesics, Vdx{v,v) = 0 and t\ 7 dy{v,v) = 0 . And as 
C^dyvdy V, v) = —{f^dyvV, dy v) for any normal v to S, we have \S/dy(v,v)\ < K\dy v\'^. 


So 




Pdy)v\ < Khd ^ 



q. e. d. 


11.3 Corollary. Situation as before, additionally \\Wi,\\ -1-/i||VllA|| < k for all Weingarten 
maps. Let iij also induce a (pD, h)-full metric on rA. Then for small h, 

\{y*9 - 9''){v,u})\ < {Coh'^ + Kh'd~'^)\v\ |u;|, (n-Sa) 

|Vf Sw - Vf w| < kl, (ii-3b) 


where Cop = (Go.i -f ^. The second estimate also holds for any other piecewise 

flat metric g® on rA. 


Proof. The metric estimate comes from the edge length comparison above, and the 
connection estimate from [6^2^ does n ot depend o n the chosen metric, as long as it is 
flat. Due to the Gau£ equation (e. g. Jost|[20iiI thm. 4.7.2), the intrinsic curvature 
tensor of S is bounded by Cq + Ikkiyp and its derivative by Ci -I- \\W,,\\ ||VIT;,||, q. e. d. 
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11. Approximation of Submanifolds 


Remark. The observation that |ii.3b| also holds for any other piecewise flat metric on 11.4 
A. means that if g is approximated up to second order by a better-suited approximation 
of edge lengths £ij than just £ij, then the approximation of the connection remains 
unchanged. 

Nevertheless, taken as it is, |11.3| says that a simple interpolation of a given triangula¬ 
tion, just as in Euclidean space, is not the best candidate for geometry approximation. 
Henceforth, the rest of this section is devoted to the normal graph mapping, which 
reveals better approximation properties. 


General Properties of Normal Graphs 

Definition. Let S' C M be an u-dimensional compact boundaryless smooth subman- 11.5 
ifold. A second submanifold S' C M is said to be a normal graph over S if there is 
a normal vector field Z on S such that 


: a !-)■ exp„ Z\o 


(11.5a) 


is a bijective mapping S —>■ S'. Where we need it, we will also consider the “smooth 
transition” S S' via the homotopy 


: a exp^tZ\a. 


(11.5b) 


Parallel transport along 11—?> ^t{p) from d>a(p) to d>b(p) will be denoted by 

Remark, (a) The term “normal graph” or “normal height map” is mostly used in the context 11.6 
of triangular approximation of surfaces in R®, e. g. in IHildebrandt et al\ ll20o6h . In the 
context of manifold-valued pde’s, it is more common to consider the geodesic homotopy 
"ht, see |i3.6dl which also IGrohs et al\ ([201^ use. In particular, their control of distortion 
along <l?t is equivalent to our control of the distortion by <E>. 

(b) Here, as usual, we do not want to treat global properties of M, so we always tacitly 
assume \Z\ < inj M. 

(c) Any other n-dimensional submanifold S' G M that is near enough to have a 
bijective orthogonal projection S' ^ S can be represented as a normal graph over S. 

Here “normal projection” means mapping some p € S' onto the point q G S minimising 
d(p, (7). The largest e such that the ort hogonal project ion Be( 5 ') —>■ 5 ” is well-defined 


IQK 


^ def. 4.1, for a recent overview 


IS c alled the rea ch of S (introduced bv iFEDERERl 
see lTHAlObnnSll . Another formulation for the same thing is that 

$ : TS-^ M, (p, Z) expp Z 

is a diffeomorphism from := {z/ G TS^: |z^| < e} onto its image. 

(d) It is well-known (see e. g. HildebrandtI I201 j . eqn 1.11) that for M = IR™, the 
map $ is locally a diffeomorphism if |Z| ||IW|| < 1 for all Weingarten maps W^- (Note 
that this bound can only capture the local geometry of S but cannot see if some part of 
S that is intrinsically far from a point p G S comes close to p in the surrounding space 
M.) Different to the usual argument involing the curvature radius of S and osculating 
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C. Applications 


spheres, one can use Jacobi fields as in 11.251 (we use the notation from there) to see 
this: 

We already know d^tP = J{t) for a Jacobi field with J( 0 ) = p and J( 0 ) = Vpv, so 
-^^*g{p,p) = g{J,J). Now let v he a, unit normal field. Then J = Wi,J. If p is the 
eigenvector in direction of the largest eigenvalue n, we have 

-^^UiP^P) = 2K(i)l-^(0l: 


where K{t) is the eigenvalue of the Weingarten map IW in direction p at $t(p). In our 
case R = 0 , the Riccati equation |i.25b| gives W, = —W^, so the eigenvalues Ki also 
evolve by ki = —nf. This differential equation has solution Ki(t) = (t — )~^ = 


Ki(0)t-1 > 


hence 


This is solved by ^lg{p,p) = (1 — Kt)^g{p,p), so g is positive definite for nt < 1 . 

(e) If we represent Z = /V^ with parallel unit normal vector fields Vi and scalar 
functions /*, then VZ = d/* 0 which splits into tangential and normal 

parts 


t\ 7 Z = fVui, nVZ = df 0 i/i. 


Hence, although ||VZ|| usually shrinks slower than \Z\ for |/| —>■ 0 , the tangential 
part is ||tV 0 || ^ K.\f\, where n is an upper bound for the Weingarten maps IW of S. 
Similary, ||tV^ 0 || < K\df \ + K'\f \ if k! bounds all ||VWy||. 


11.7 Situation. Let S' be given as a normal graph over S' by a vector field Z with d := 
\Z\ < and ||dd || < e everywhere, and let the Weingarten maps of S be bounded by 
||Wi/|| + e||VWj/|| < K. This means ||V 0 || < ke. For simplicity, assume k < k^. 


11.8 Proposition. Situation as m[TTT7| The map $ : a 1—>■ exp^Z is locally a diffeomor- 
phism if \d\{K + y/Uo) < 1 everywhere. 


Proof. Let us suppose Z has unit length at some point and see for which t the map 
is locally a diffeomorphism. Note that |i.2,t^ gives \ki\ < Cq + |k^| for any eigenvalue 
of a Weingarten map. The equation u = —Co — leads to a subsolution 


u{t) = \/^tan (^\/^t — arctan . 


Regarding -^^*g{p,p) < K{t)^*g{p,p), which has positive solutions (for positive initial 
data) as long as k is bounded, it suffices to know where the first pole of Ki(t) can occur. 
The first pole of u (which must also bound the position of the first pole of Ki) is where 
i/Cof — arctan • Now a simple function inspection shows 

—— < ? + arctan s, -^— > —-I + arctan s, 

1+s 2 1+s 2 

so there will be no pole as long as \'/Cot\ < (1 + 9- e- d. 
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11. Approximation of Submanifolds 


Observation. Bv li.i 61 the differential of is d^tV = J{t) for the Jacobi field with 
J( 0 ) = V, j( 0 ) = VyZ. Bv l6.ll this means iov Qt - V + WyZ 

\d<^tV - P^’°QtV\ < Co + Ket)\V\ (11.9a) 

and 

\d^sV-P^’°V\ <Ke\V\+Cod^t^{l + Ket)\V\. (n-gb) 

This estimate is scale-invariant with respect to scaling of Z: If Z' = aZ, then = 
So tZ = t'Z' is scale-invariant, and the estimate only contains tZ, never Z alone. 
But due to 111.Bel this distortion happens mostly in normal direction, the tangential 
change is of higher order: 

||t((i 4 >-P)|| < Atd-fCod^ (11.9c) 

Proposition. Situation as in 1 11. 7( Consider some point p G reach S with projection 
if{p) onto S. If p = exp^(p) dv for some unit normal vector v with \Wu\ < k, and if 
KE < ^, the orthogonal projection if satisfies 

\\diP - Q-^ tpp-’/’fp) II < Co d^ Wdif - II < K d +Co d^, (ii.ioa) 

where Qi, is the linear map TpS — > TpM , V V + d If dv is replaced by some 

other normal vector field Z with exp^j-p^ Z = p, and Q is replaced by V 1—>■ P 4- VyZ, 
it holds 


\\diP-{tQ)-HpP’^<^P^\<Cod'^, \\dif-tPP’'<^<^P^\<Kd+Cod\ (11.10b) 


Proof. Let us first show that tQ does not depend on how Z is chosen at points 7^ ip{p). 
If Z = d 1/ in a neighbourhood of if{p), where z/ is a parallel unit normal field and d is 
constant, then Z is parallel, and so VZ = tVZ = dW,^. For any other Z, tVZ stays 
the same, and only some part nVZ 7^ 0 is added. That means tQi, = Qi, on T.^(^p'^S. 
So we will only prove fii.iobi 

ad primum: Observe that the operator tQ : TpS —?► TpS fulfills || tQ — id || ^ k£ < ^ 
by assumption, hence is invertible with ||(t(5)“^ ~ id || ^ by |6.i5[ which gives 

||(tQ)“^|| < 1 - 1 - < 2 . Hence the claim is proven if we can show || t{Qdilj — 

P)\\<Cod^. 

Consider some vector V € TpM and split V = V), -I -14 as in |i. 2501 Then tQdif{V) = 
p -I- tVpZ = Jp( 0 ) -I- jp( 0 ) on the one hand, and tPV = PVp = PJp(I) by |i.25d| on 
the other. So 16 . 1 1 gives 

|Jp(I)-P(Jp(0) +jp(0))| <Cod2 IpI <Cod2|H|. 

ad sec.: We have \dif — tP\ < Idif — (tQ)“^tP|| -|- \\{tQ)~^t — t||. The first norm has 
been estimated above, and the second is < Kd because ||tQt — t|| < k d due to lii.6el 
and the boundedness of q. e. d. 


Remark. For Cp = 0 , this (exact) representation of the pr o jectio n differential is the one in 
IWardetzkyI (l2oofil . thm. 3.2.1) and iMoRVAN and Thibet ( 2nnJ . lemma 4). 


81 


11.9 


11.10 




























C. Applications 


11.11 Proposition. Situation as in \ii.'/\ Let t' he the orthogonal projection TM\s' —>■ TS'. 

Then for small e, we have \\Pt' — tP\\ ^ ks + Cod^. This means that the angles 
A{T^(^p'jS', PTpS) and PTpS^) between the corresponding tangent and 

normal spaces must be bounded by this factor, too. Therefore, normals Vi to S can he 
extended to normal fields Vi^t along with \vi^t — P*’^Vi\ ne + Cq . 

Proof. For the time of this proof, let us write the terminal value J(l) of a Jacobi field 
along t I—>■ $t(p) with initial values J(0) = q and J(0) = z> as T{q,v). Linearity of 
the Jacobi equation translates into linearity of T. In this notation, the splitting from 
II.25CI says that a vector V G TM\s> can be represented as P = T{p, tv) + T(0, nv). 
We argue that its projection t'V onto TS' is almost T{j),VpZ). 

In fact, all tangent vectors on S' have the form T{q,VqZ) for some q G TS. Now 
consider 

\V-T{q,VqZ)\^ = \T{p-q,v-VqZ)\^. 

This is minimal among all f if t'V = T{q,'VpZ), this means its norm has vanishing 
derivative in direction (r,'Vf.Z). Because T is linear, this gives 

0 = (T(p -q,v- VqZ),T{f, VrZ)) for all f G TpS. 

Now recall that T{U, W) = P{U + d W) + 0{Co d^), hence this is 

= {p-q,r + d VrZ) + d{v - VqZ, r) +d^{v- WqZ, VfZ) + 0(Co d^). 

li p = q, the first term vanishes, and (using that k d is small) the remaining ones are 
estimated from above by Kd|P| |f|. Because the minimisation is well-conditioned at 
this position, the optimal q is p + 0 {{k d +Cq d^)|P|). 

Now recall from|i.25d|that PtPV = T{p, tVpZ), which gives that the claim \Pt'V— 
tPV\ = \{t' -PtP)V\ = \T{Q,nVpZ)\+0{{Kd+Cod^)\V\) < {Ke + Cod^)\V\ is just 
the usual Jacobi field estimate O q. e. d. 

11.12 Corollary. Omitting the last paragraph of the proof, one gets || tPt'— tP\\ < k d +Co d^. 

Remark. This is analogous to the classical statement \\P{Ph — 1 )P|| < d up to constants 
depending on the geometry from Dziuk et al, where P is the projection onto TS and Ph the 
projection onto TS'. 


Geometric Distortion by the Graph Mapping 

11.13 Lemma. Situation as in \ll.f\ Then for Q : U t-G U + VjjZ, 

\{QU,QV) - (U,V)\ <{K^d+Cod'^)\U\\V\ for all U,V G TpS. 

Proof Just because {U + VuZ,V + VyZ) - {U,V) = {VuZ,V) + {VvZ,U) + 
{VuZ,'S/vZ) and pi.gci q. e. d. 

11.14 Conclusion. Situation as in \ll.f\ Pulled back to S, the S' metric g\p{U,V) = 
g\iS>{p){d^U,d^V) fulfills 

I ($*5 - g){U, V) I < d +Co d^) \U\\V\. 
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11. Approximation of Submanifolds 


Proof. Is a direct application of [TTTg^ and [TTT^ We especially remark that the differ¬ 
ence between g\p and g|$(p) does not need to be handled explicitely, as pP’'^^P'> is an 
isometry with respect to these two metrics, q. e. d. 

Remark. As ^*g and g are equivalent metrics, A := (where dtl?* denotes the g-adjoint 

of d$) is a self-adjoint autom orphism of T„ S such that ^*g{U,V) = g{AU,V), called the 
metric distortion tensor by I WardetzkyI (l20o6l . p. 53). In the numerical literature, it is 
common not to compare the Riemannian metrics, but to estimate di rectly ||^A — i d ||, which 
already includes the volume element change (cf. the proof of [7!^ , see lPziUM (fioSSh : IDemlo'VS^ 
(1200 3);1 EINeI ll20oJ ). 

For a comparison with the tensor J from 16 . 2 ^ consider 'I' := If M is the Euclidean 

space IR” and S' is a piecewise flat submanifold, its metric g" := pis/ is piecewise flat. The 
metric p|s pulls back to a metric 5 '*p on S', and there is J such that ^*g{U, V) = g’'{JU, V). 
So the transformations A and J perform inverse tasks. 

Proposition. Situation as in \ii.'y\ For a given vector U and a vector field V on S, 
define the “connection distortion” W := This vector field 

obeys the differential equation 

W = R{Z, W)Z + F forF := i?(Z, d^tU)d<PtV + 

with initial values W{Q) =0 and W{ 0 ) = F{ 0 ). 

Proof. Let us abbreviate U* := d^tU, := d^tV, and denote the parallel translation 
of Z along t I—>■ expptZ also as Z. Let K := VptI/* and J := Then we 

want to determine W = K — J. 

By |i.25[ J is a Jacobi field, i. e. J = R{Z, J)Z. An inhomogeneous Jacobi equation 
describes K: 

K = R{Z, K)Z + DtR{Z, U*)V* + DtVltytZ. 

In fact, consider a variation 7(r, s) of geodesics (in M), i. e. we assume that s H/ 7(r, s) 
is a geodesic for each fixed s, with 9^7(0, 0 ) = V and 9^7(0, 0 ) = C/. Transport this 
along t as c{r, s, f) := exp..^^^ tZ\.yi^.,. gy Then we want to determine K = Drds, so we 

consider 


K = DtDtDrds = DtDrDtds + DtR{dt, 9,)9« 
= DtDrD,dt+DtR{dudr)ds, 


the first term of which is 

VzVptVyt Z = V ytZ J- 

— V yt Z + \7y^^Yt\/zZ + R{Z,VutV*), Z 
= DtWfjyytZ + 0 +R{Z,K)Z. 

The initial value is computed in exactly the same way, q. e. d. 


11.15 


11.16 
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C. Applications 


11.17 Proposition. Situation as before. J/Co|^P + ||VZ|| < then 

tVuV] < \U\\V\{\\V^Z\\+Co\Z\) + Co\VuV\\Z\^. 


If we only consider the tangential part tW, then even 

tVc/P| < \U\\V\(\ty'^Z\\ + \yZ\\+Co\Z\^)+Co\yuV\\Z\'^. 

(11.17a) 

Proof. Preparatory step one: Let us first establish the boundedness of VL = K — J and 
show |fL| < \VuV\ +at, where a := |C/*| |F‘| (||P Co\Z\): The t-derivative of LiT 
is, as dt = Z ., 

S/z^u^V* = Vu^VzV^ + R{Z, U*)V* 

= Vu^Vv^Z + R{Z,U*)V* = Wlj,y,Z + Vv^,v^Z + R{Z,U*)V\ 

so < \K\ < O'+llVZI \K\. As \Z\ is short by assumption, we have \U*\ < |f7| and 

\V^\ < \V\. The differential inequality of the form ft < a + bu gives u < (uo + f )e^* — f. 
For bt < ^, this function is dominated by uq + 2bt{uo + f ) < 2(uo + at). 

For the bound on J, we have | J| < | J(0)| + t| j(0)| as usual, and J(0) = Vf/P, 
J(0) = V 2 J(0) shows that these terms are already contained in the K estimate. 
Preparatory step two: Now let us show 


Z 

W{t) - J P^’^F{t) dr] < Cot^\Z\^{\VuV\ + at). 


The proof idea is from ljosil ( 201 iL thm 5.5.2). Let A := PF. This is a vector held 
fulhlling A = F with the same initial values as W, namely A(0) = 0 and A(0) = F{0). 
Furthermore, let w : [0;t] —^ R be the solution of «) = Col^PllTl with initial values 
r/;(0) = i(;(0) = 0. Then, for some parallel vector held E along t, dehne v := {(W — 
A,E) — w)/t and obtain 

^ (ut^) = -^(((bF — A,E) — w)t — (W — A,E) +w) = iiW — A,E) — w)t < 0 . 

df ^ ^ 


dt 


This means that vt'^ < 0, hence u < 0. Now {W — A, E) — v has a double root at t = 0, 
so v(0) = 0 and thus ?; < 0 everywhere. And because E was arbitrary, this already 
means 

\W — A\ < w. 

Therefore we are done if we can bound u by the right-hand side of the proposition. 
But as we know that |kF| < iVf/P] -I- at, we can simply integrate w = Co|^P|VF| twice 
and obtain the desired estimate (the argument is the same as in the proof of|^- 
ad primum: The hnal estimate for |kF| comes from |/ PE\ < / |A|, together with 
l^(^)l ^ l^(0)l < a for t < 1, because the same holds for and V*, and the norm of 
is the same as the norm of V‘^Z\p because Z is parallel along 11—>■ d)t(p). 
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11. Approximation of Submanifolds 


ad sec.: For the estimate of | tW\ let t* be the orthogonal projection onto the tangent 
space of $t( 5 ') and consider 

t t 

dr = + dr. 

0 0 

Then |tPP| < | 17 | |y| || tV^Z||, and the projection difference is estimated by tllVZl + 
Cot^\Z\^ in iii.in q. e. d. 

Theorem. Let y : rA S be the triangulation of a smooth submanifold S G M 
with Weingarten maps bounded by \\Wi,\\ + /i||VlTy|| < k and x : rA. —>■ S' an 
extrinsic Karcher triangulation with the same vertices pi and y = ip{x). Suppose is a 
{id, h)-small piecewise flat metric on rA. induced by edge lengths ds'ijii,Pj) = A{pi,pj). 
Then for small h, it holds d(a:, y) < h'^id~^ || Vda; — P^dy\\L°° and 

\{y* 9 - f){v,w)\ < (KM+Cod^jlwllwl, 
iV-j — V® wl < K,h-d~‘^ \Vdx — Pydy\\_<=^ |t;| |u>| + ho., 

where “ho.” stands for higher-order terms whose coefficients depend on Cq, k, d, |z;|, 
|w| and |V® w|. The norm on the left-hand side may be induced by either x*g, y*g, or 
g®, because all three are equivalent. 

Proof. By the estimate |7.1:2| for d(x,?/), S' is a normal graph over S for small h with 
\Z\ = d{x,y) < h^'d~^\\ydx — PVdyj. Morally, it is clear that ||VZ|| must be con¬ 
trolled by \Vdx — PVdy\, too. In fact, we can precisely compute this for V = dyv. 

S/yZ = y dyv{ — Xx\y) = ^(1) 

for a Jacobi field along x y with J( 0 ) = dxv and J(l) = dyv by combining 11.231 
and 112.31 (’"^6 postpone the computation to the next section because it is more relevant 
there), so bv l6.1l |Vi/Z — {dx — Pdy)v\ < Cq d?{x, y), hencegives 

||VZ|| < hd-^ \\Vdx - PVdyW -f Co d^ 

Now because x = ^ o y and hence dx = d^dy, [TTTT^ gives 

\{x*g-y*g){v,w)\ < {n"^ d+CQd?)\v\y*g\w\y*g. 

The comparison of x*g and g" is done in I6.17I — Analogously, we only compare ® 
to ® and refer to |6.23| for the rest: For a vector v and a vector field w, | 11 . 17 | gives 
(together with In.Gel) 

\V^„^.„dxw — d^V^yj^dyw\ < Kh'^d~^ jVdx — PVdyW -f- ho. 

By definition of the pull-back connection, X^^^^^dxw = dxVi^ and thus 
|dd>d2/(Vf - Vf 9 ie)| < nh'^d-^ ||Vda; - PXdy\\ + ho. 

Together with < 1 and ||(J2/“^|| ^ 1 /hd, this shows the claim, q. e. d. 


11.18 
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C. Applications 


The Weak Shape Operator 


11.ig Lemma. Let S G M be a smooth submanifold with boundary dS in M. Let Lf be a 
smooth vector field on S, not neccessarily tangential to S. Then Lf may be extended to 
a vector field on M in such a way that div^ tU = div'^ tLf and div'^ nLf = —{U,H), 
where H = for any orthonormal basis Ci ofTpS is the mean curvature vector 

ofS. 

Proof. It suffices to find a local extension of [/ to some small neighbourhood of S. Let 
Cl, , e„, Vn+i, ..., Pm be an orthonormal basis of TpM . Then div^ Lf = (V^Lf, ef) + 
(ffJ^.Lf, Vj). If Lf is extended parallel in normal direction, the latter term vanishes. 

Regarding the tangential part, observe e^) = (tV^tt/, e^), and because 

= V^, this is div'^ tU. 

Now consider nU = a^Vj. Again, if U is constant in normal direction, div*^ nil = 
{{dia^)iyj, ei)+{a^V^ Pj, ef), the former term vanishes, the latter one is ef) = 

—(V^:^ei, Pj), q.e.d. 


11.20 


Lemma. Let S G M be a smooth submanifold with boundary dS in M. Then for 
smooth vector fields V and W on M, 

J {W, p) div"^ R + {W, p)(R, H) + p, W) + (p, Vf IT) = J (IT, p)(R, t), 
s as 


where r is the outer normal of dS in S. 

Proof. Let / := (IT, p). By the divergence theorem ( Lee! 200^ thm. 14.23), we have 


I div^(/R)= l{fV,T). 
s as 


Now by product rule and |ll.fgi div'^(/T) = f div^ V + Vf = /div^ V — f{V,H) + 
Vf, q.e.d. 


11.21 Corollary. Let S G M be a smooth submanifold without boundary. The operators 
SvtCTv '■ 'X.{TM\s) X X{TM\s) —>■ IR, defined by 

s,(TIT):= y'(V(^p,IT) + (IT,p)(T,R'), a,(T, IT) := - ^(IT, p) div'^V+(p, VfiT) 
s s 

eoincide for smooth S and each normal field p. On tangential vector fields, Si,{V, W) = 
CTiy(TIT) = — f II;y(T, IT). If S were only piecewise smooth, a^, would still be well-de¬ 
fined. It is called the weak shape operator of S. 

11.22 Proposition. Situation as in^TTT^ Then there are normal fields v for S and v' for 
S' such that approximates the weak shape operator cr),, of S', which is 

a),, (T W) = -J (IT, p') div"^ T + (p', Vf IT), 

S' 

up to first order: |(dj,(T, IT) - cr),,(PT, PIT)| < e|T|Hi(TM|s) |bb|Hi(TM|s)- 
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12. Variation of Karcher Simplex Volume 


Proof. The derivatives of V and W are not distorted at all, PVylT = ’S/pyPW, see 
|i.i7[ hence div^ PV = div^ V. According to in.in the normals n and n' can be 
chosen such that they do not differ by more than e, q. e. d. 

Rema rk. This is our analogue of the weak s hape operator conver gence result from iHiLDEBRANDT and Polthi^ 
(l20iil . thru. 8) (similar, but more extensive. iHlLDEBR.AND'Jboi j . thm. 2.4). However, Hilde- 
BRANDT defines cr(,/ with div® instead of div**^, which did not meet our own calculations. 


12. Variation of Karcher Simplex Volume 


Notation. Different from other sections, we will mostly write volS” instead of for 
the volume of some set S C M, but where we find it meaningful, we will mix both 
notations. 


Goal. Consider a full-dimensional Karcher simplex x(re) with vertices pq, ... ,pm- If 
Pi is moved with velocity pi = dxui, leading to a family xt of barycentric mappings, 
we would like to compute the derivative -^|t=o volg(a;t(rs)) of a subsimplex volume. In 
a second step, we want to show that this derivative is close to -^|t=o volge (rst), where 
the vertices of rs are moved with velocity ut. 


Fact (see e. g. I.Iost 


2011 


eqn. 4.8.3). Let <i)t : S' —>■ S'* be a normal variation of the 
smooth submanifold S C M with a velocity Z = dt^t that has compact support. If H 
is the mean curvature vector of S, then 


12.1 


d 

dt 




volS* = - 


div'^ Z 





where the last equality only holds if Z vanishes at the boundary of S, whereas the 
divergence formula is also correct for velocities with support up to the boundary. 

Situation. Let A be a regular m-dimension simplicial complex, a; : rA —>■ M be a 12.2 
Karcher triangulation with respect to vertices pi £ M such that driX^ is bijective for 
every vertex of every element and x induces metrics x* g and as usual, where is 
(i?, h)-small, fulfilling I6.17I and |6.23! Let A be an n-dimensional subcomplex of A such 
that S := x{rR) is an n-dimensional submanifold of M. 


Variation of Karcher Triangulations 

Lemma. Let p £ M be some point and X be its half squared distance gradient from 12.3 
[37^ Recall from that if the evaluation point a is moved on a eurve a(r), then 
= aJr { cr ), where Jr is a Jacobi field along the geodesie p a(r) with Jr { 0 ) = 0 
and Jr{er) = a{r). 

Now if p moves on a curve p{t), this induces a new variation veetor field Xpi^p^ := 
DtXp(^iy It satisfies = aJt{(j), where Jt is a Jacobi field along p(t) a with 

boundary values Jt( 0 ) = p{t) and Jt{cr) = 0. Combining both variations, we obtain 
^ d{r)^p{t) — (jDrJris^^ — erDiJr^ei^ 
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C. Applications 


Remark, (a) Note that J and X denote derivatives with respect to different direc¬ 
tions. In this section, the parameter r takes the role of t in section [51 

(b) In the notation oI IGrohs et al\ ll2ni j h X = Vilog(a,p). 

Proof, ad primum: Without restriction, we can assume that p(t) and a{r) describe 
geodesics, as no second derivatives of them enter. Analogous to |l.S3[ consider a varia¬ 
tion of geodesics 

c(r,s,t) := expp(t) (f(expp( 4 ))-ia(r)). 

As in|l.:53[ one can see that 


AppJ |c(r,s,t) — (j ^sC(r, S, t), Jr — dr^ Jt — Qt 


because dr and dt are Jacobi fields with the desired values at s = 0 and s = a. 

ad sec.: We have ^'^a{r)Xp(t) = DrDtdsc{r,s,a) = DtDrdsc{r,s,a) + R{dr,dt)ds, 
and because c(r, s, a) = a{r) is independent of t, we have dt = Q there. And the former 
term is DtDsdrc{r, s, a) = DtJr[cr), q. e. d. 


12.4 Lemma. Situation as before. Abbreviate I := d(a(r),p(<)), V := d(r), W := p{t) and 
T:=ds. Then if Cq^ < 


\DtJr\ < 90Co.i(r)2(2^ |I/| |W| |T|(r) 
^gOCo.iWsiyilWI |T|(r), 


iDjrl < 50Co,i(r)|y| |W| \T\{r). 


Proof. The claim is similar to 16.61 and so is the proof: We again have to find a differen¬ 
tial equation for U := Dt Jr to apply [^7^ By the usual laws of covariant differentiation 
and dr = Jr, dt = Jt, we have 


DfJr — DtDsDsdr — DgDtDgdr J- R{Jt, T)jr 

= DgDgDtdr J- DsR{Jt, T)Jr + R{Jt, T)jr 
= DlU + {DgR){Ju T)Jr + R{jt, T)Jr + 2A( Jt, T) j,. 


On the other hand, using the Jacobi equation, 

-DtJr = DtR{Jr, T)T = {DtR){Jr, T)T + R{U, T)T + i?( J„ jt)T + R{Jr, T)jt. 

So with A := —R{ ■ ,T)T, we have U = AU + B, where ||A|| < Co\T\'^. Let us assume 
that we consider some geodesic with |T| = 1 (as usual, the correct power of |T| follows 
from a scaling argument). Then, because I 6 . 3 I holds for Jr as well as for Jt, 

\B\ = |(iJ,i?)(Jt,T)J, + {DtR){Jr,T)T 

+ i?( Jt, T)Jr + 2A( Jt, T)jr + R{Jr, jt)T + R{Jr, T) Jt | 

< 2Ci\V\ |W| + l^Ct)^\V\ |W|. 


The claim on DtJr follows from DtJr = DtDgdr = DgDtdr + R{dt,dr)ds, q. e. d. 
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12. Variation of Karcher Simplex Volume 


Lemma. Notation as before, and A\ as For varying vertices Pi(f) G M, let Xt 

be the eorresponding Karcher triangulations, and V = dxv. Then 

Axitlx = -X^X,\,,(x), AxVvi = -v^X, - XV\^^X, - XVyX, - A„x. 

Proof, ad primum: Along t i—>■ Xt{X), consider the vector field U{t) := X^Xi{t)\j;,_(^x)- 
As has been stated in [57^ this vector field vanishes for all r, so U = 0 . For those who 
believe, the shorthand proof is 

0 = A|t=o(A*X,(t)UJ = AA|t=o(^*WUo) + A*A|t=o(^z(0)UJ 

. . (12.5a) 

a rather uncommon use of the chain rule. For all others, this argument is justified by 
a calculation in coordinates: Let X^Xi = U^dj and x{t) = x^{t)dj be the needed coor¬ 
dinate representation. As U^if) = 0 for all k and all t, we have U^{t) = 0 , and this is 
by usual Euclidean chain rule dtU^{t, x^(t ),..., x"‘{t)) + d(,U^{t, x^{t ),..., x^{t))x^{t). 
And if all vanish, we can add T^jU^x^ without harm, which gives 

0 = (dtU^ + diU^x^ + r’^jWx^)dk = DtU + ViU. 
ad sec.: Differentiating [T 27 p[] once more leads to 

0 = AVyViX, + Vy(AA,) = + XVvy±X, + NX, + XXyX, + NX±X„ 

q. e. d. 

Proposition. Situation as in \i2.2i and let the variation ofpi be given bypi(t) = dxwi 
for some vector Wi € Tre. Then for u := X^Wi, we have 

\x — dxu\ < Cq \Xdxvi — dxV^ m| < Cq ih\u\ \x\. 

Proof, ad primum: At the vertex Ci of the standard simplex, de, x w and the variation 
Pi agree. At another point A € A, 

dxXWi = PdeiXWi + 0{CQih‘^\wi\) = Ppi + 0{CQih?\wi\) 

by |6.24[ which means dxu = X^Ppi + OiC^ ih\u\) by definition of u, and 

X =-XXi + 0 {CQh^\x\) by [T^ 

= X Ppi + 0 {C'Qh^\x\) by [T^ and[6(^ 

ad sec.: The derivative of u is, by Euclidean calculus, just V® u = v'^Wi, so the latter 
term is dxV® u = v'^dxwi. The covariant derivative V^Xi in |i2.5| is estimated by |6.r4[ 
and the VyX term by |r2.4[ so we have 

\AxXdxvi - v^dxwil < v^\Xi + dxwi\ + X'^\Xyd:Xi\ + A*|VvAi| -|- \Ayx\ 
<Coh'^v"\pi\ +CQ ih\v\ \x\, 

q. e. d. 


12.5 


12.6 
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C. Applications 


Discrete Vector Fields 

12.7 Definition. For a piecewise barycentric mapping x : rK ^ M, let 

PXx := TM\^(^si°) = 

Pi 

be the disjoint union of all vertex tangent spaces. For U = {Ui) S PX^, define a 
piecewise interpolation: It induces a variation of x by defining xt [U] to be the piecewise 
barycentric mapping with respect to vertices exp^. tUi (where we keep t so small that 
x(A) and Xt(A) stay in a convex neighbourhood of each other). We call U U := 
a;t[i 7 ]|t=o the P^-interpolation of U and 

P^X, :={[/: C 7 g PXJ 

the set of piecewise smooth, globally continuous test vector fields. 

12.8 Observation. As a finite sum of vector spaces with scalar products g and 5®, PX^, 
carries the natural inner products 

fg{V,W)=Y,9{ViiW,), 

i i 


whereas P^X^ has the scalar products L^g and that are induced from L^X(TM): 
L^g{V,W}= [ g{V,W), LV{t",PP)=E [ 

J an Jre 


x{r^) 




where V = dxv and W = dxw. As both are isomorphic finite-dimensional vector 
spaces, all these norms are equivalent. The equivalence constants between i'^g and i^g^ 
and between \J^g and are the ones from |6.17a] and whereas the equivalence 
constant between and I'^g^ depends on the maximal and minimal simplex volume. 

12. 9 Definition. Situation as in 112.21 and U G PX^,. Inside every simplex e G the 
vector Ui, i G e, can be represented as Ui = driX^ wl (without any summation). Define 
a piecewise linear, globally discontinuous vector field u\re ■= A^w) and a piecewise 
smooth vector field u by requiring dxu = U everywhere. 

12.10 Conclusion. By definition, dxV^ = '^dxvit[U]. From, [TEUl we hence know that 

\u-u\<Ci,h\u\, \yf^u-y{u\<C'o^,h\v\\u\, 

where all norms are | • |a:*g norms. The same estimates hold for the jump [dxu]^ = 
dxcX^'wl — dx^iX^wl of u across a facet f = e fl c'. 
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12. Variation of Karcher Simplex Volume 


Area Differentials 

Observation. Situation as in I 12 .il If the vertices pi of a Karcher triangulation vary 
smoothly with velocity (Ui) € PX^, the area change of S' = x(r^) is also smooth, hence 
has a differential 

dvol^ : PXx —>■ R- (i2.iia) 

The volume is additive in s G it”, and the variations of different vertices are linearly 
independent, so it suffices to compute the differential dvol* ^ : T^irK —>■ R of |x(7’s)|g 
with respect to the variation of x{ri), i G s. Correspondingly, let dvol* be the 
analogous differential of |s|ge. 

Remark. We do not think that a notational distinction between this area differential 
and the volume form from [g.6] is neccessary. For readers who disagree, we remark that 
in ir2.iiaL the d denotes a differential and is hence written in italics, whereas it is 
upright as part of the volume form dvol. 

Proposition. Situation as in \ 12 . 2 \ Then |dvol* ^ —dvol*_ge| < Cg ^h|s|ge- 

Proof. Bv ii2.h dvol* giw) = — f div'^ Z, where Z = x is induced by the vertex varia¬ 
tion pi = dxw. livj form a p-orthogonal basis ofTrs, this is f (Vdxvjijdxvj). By the 
comparison of volume elements for g and in 13.2m 

dYollg{w) = -J g{ydxvji,dxVj)=-{l + 0{Coh'^)) j g{Vdxvji,dxVj), 

x{rs),g x{rs),g‘‘ 


and, noting that there is a p®-orthonormal basis vj of Trs with \vj — Vj\ < C^h'^ by 
|3.6[ the integrand is 

gi'^Lv^i^dxVj) = X* giV’l.u^Vj) + 0 {Co^ih\u\) by [ 12^1 

= x*g{V^;^u,Vg) + 0(C',ih|u|) + 0(C>2||Vu||) 

= + 0 {C'o,,h\u\) + 0 {C'oh^\\Vu\\), 


and the last right-hand-side term is div^*^^’® ^ u, 


q. e. d. 


Remark, (a) It is common knowledge that dvol^ = dA®|A|, proven b y insertini 
div u = dX''{w) for u = X’'w into rTXn Classically, one says for triangles ( Polthie 


2002! eqn. 4.3) that the gradient of the area functional with respect to vertex variations 


is the ^ rotation of the opposite edge vector, which is exactly (dA®)**. 

One has to take care to transfer this to the subsimplex situation, because div*-’’®’® ^ u 
= 0 if w is perpendicular to rs, so one needs a form dA® that acts like dA® on Trs and 
vanishes on Trs^. For example, if rs = conv(eo,..., e„) and vq = grad A° € TA, 


dA® (w) = dA® 


{wo,w) 


-Ug = W 


j£j,{dX\dX°). 


12.11 


12.12 


12.13 
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C. Applications 


It is easier to transfer the gradient u* of A* in rs to Trc, as its vector components stay 
the same: div*-^’® ^ u = g^{vl,w). By |3.2a[ this v\ is characterised as the vector in Trs 
that is perpendicular to the facet s \ {i} opposite i with lengths h~^, the reciprocal of 
rs’s height above s \ {i} from [gTSaj 

(b) For computational purposes. ir2.i2i is insatisfactory, as only dvol^^g<i(u) is numeri¬ 
cally accessible, not d vol^^ge (u). But this is only an easy combination with rrXriTl which 
will be spelled out for the area gradients in the following paragraphs. 


Area Gradients 


12.14 


12.15 


12.16 


12.17 


12.18 


Observation. The area differentials dvol^.g and dvol^^ge can be expressed as gradient 
with respect to different norms on PXx- Th e gradients with r espect to £^g and £^g^ 
correspond to the “mean curvature vector” of POLTHI^ l 2002h. whereas t he gradients 
with respect to or L^g and L^g^^ give the construction from lPziUKl l|iQQih . 


Corollary. If Hpg G Tp.M is the gradient of\x(r5)\g with respect to a variation ofpi, 
and H(2ge = |s|^e gradA* is the corresponding gradient of |s|ge, then \dx Hpge — Hi2g\ 
^ C'q I. 

Definition. The discrete mean cnrvatnre vector H\_2g G is the solution of 

L^g(Hi2g,V} = dvol^g(y) for all V € P^X. The approximate mean curvature vector 
H]_2ge G P^'X.x is the solution of \-^g^{H]2ge,V) = dvol^ge(F) for all V G P^X. 

Observation. The differentials of the right-hand sides can be represented as prod¬ 
ucts of linear maps: If vectors Zi = dxwi G Tp.M induce a variation x = dxz oi x, 
and if we define z = X^Wi, as well as dA : 1—>■ gradA®, then 

L g^Il\_2g., Z'^ = (dx, V ^z)\_2gpj'x^e^x*TM)^ (12.17a) 

Z) = (dA, V® z)iGg’^{TrSi®x’‘TM)- (12.17b) 

To see |l:2.17a[ we take an orthonormal basis Ei = dx pi in div'^ Z = (X/EiZ, Ei). Then 
we have an integrand of the form {a{yi), Pijji)) for two linear maps a, j3. A computation 
in coordinates easily shows that this is {a,l3). 

For |i2.i7b[ the computation is even simpler: z = A*Ui (without summation) for the 
gradient vt of A* has derivative Vz = Vi® dA\ which maps to Vk^ so divz = u) = 
(VeiZ,grad A‘). 

Proposition. Situation as before. Then \Hi2g — Hi2ge\i2 < Cg ^^(1 + h\Hi2ge\i2). 

Proof. Similar to 110.13; The functionals on the right-hand side of |i2.i7b| only differ 
by a factor of 1 -|- 0 (Cq ]^/i|S'|), and the bilinear forms on the left fulfill \L'^g{U,V) — 

l^gHU,V)\<C'oh^\U\ |V|,so 

\H\_2g Tl\_2ge\y^2g ^ P {^H \_2 g H \_2 gB ^ H y_2 g [_ 2 ^ 6 ) 

< \L'^g{Hy_2g,Ht2g - H^2gB) - L^g<^{H^2gB,H^2g - H^2gB)\ 

+ — L^d*)(^^L2ge,ddL2g — H]_2 ge)\ 

^ K'^VOlj^.g -C^VOlj|_ge)(idL2g - Hy_2ge)\ 

+ C'^h^ |ddL2g<= I \H\Jg - H\2g‘ I 

< C'Q ^h\H^2g - H^2gB\{l + h\H^2gB\), q. e. d. 
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ig. The Manifold-Valued Dirichlet Problem 


13. The Manifold-Valued Dirichlet Problem 


Goal. Let and Mg he two smooth compact Riemannian manifolds of dimension 
n and m. We have seen in 18.81 that a sufficiently dense generic point set in N gives 
us an almost-isometry N —>■ and in 17.141 that smooth functions rK ^ M can be 

interpolated by piecewise barycentric mappings rs —^ M, s G M", that are globally 
continuous. Now it is a natural attempt to consider Galerkin approximations e. g. to 
the Dirichlet problem in H^(N,M): Suppose iV has a Lipschitz boundary that can 
be resolved by the triangulation. Given a smooth harmonic function y : N ^ M, let 
X : —>■ M be the function that minimises the Dirichlet energy among all functions 
that are piecewise barycentric mappings and agree with y at the boundary vertices. 
How well is y then approximated by x? 

A popular example for such functions y is the minimal submanifold problem: If 
S' C M is a minimal submanifold, then the identity mapping S —>■ M is harmonic 
I JostI 2011 . eqn. 4.8.12). One already sees that to keep the notation consistent among 


the chapters, we denote smooth harmonic functions etc. by y and the interpolation in 
the sense of|7.6|by x. 


Rema rk. The results presented in this section are generally the same as in IGrohs et al\ 
(l20iql l. but although their interpolation procedure is the same as ours (but including high¬ 
er-order interpolation, see ly.idd] ), their functional analytic approach is slightly different: They 
do not use the distance measure pi and its corresponding Poincare inequality, but consider 
functionals that are (in a weak sense) convex along geodesic homotopies. 


The General Galerkin Approach 

Definition. The Dirichlet energy oi y G C^{N,M) is i3-i 

Dir(2/) ■=\J MyP- 

N 

Recall that the norm on TN ® y*TM induced by 7 and g has a representation in local 
coordinates for N and v® for M as \dy\^ = gijy^aU^p- 

Proposition (Dirichlet principle. IJostI I2011L eqn. 8.1.13). y G C^{N,M) is a 13.2 
critical point for Dir with given boundary values iff {dy^W) = 0 for all compactly 
supported vector fields V G X{y*TM). 

Proof. Every compactly supported vector field V along y induces a variation yt := 
expy{tV) of y that does not change the boundary values. By the usual calculus of 
variations, |t=o Dir(?/t) = {dy, Dtdy). So we only have to compute the f-derivative 
of dyt = y'‘a{i)di ® du'^ with the notation from li.fil As Dt = Vv in the usual sloppy 
notation and Dt{du°‘) = 0 , we have at the origin of normal coordinates in N 

Dtdyt = Dt{y]^di) ® du°‘ = (g) du°‘ + y]a'^vdi 0 du‘^ 

= ® du^ + y]^VM%dk ® du^ 

On the other hand, regarding I i.6d 

VI/ = {V^d,) ® du^ = 0 + V^y^T’^ydk 0 ^ 
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C. Applications 


13-3 Definition, (a) For N = rA, let P^{N,M) be the space of piecewise barycentric, 
globally continuous mappings (which obviously depends on the simplicial structure and 
not only on its manifold structure, but we do not explicitely denote this). Obviously 
the domain of Dir can be extended to include also piecewise smooth mappings, so every 
V G P^{N,M) has finite Dirichlet energy. 

(b) For a,b G M, denote a ~ 6 if there is a unique shortest geodesic a 6 in M. Say 

that x,y ■. N ^ M are close of x{p) ^ y{p) for almost p G N. 

(c) On C^{N,M), define the L’' metric po.r and the H^’’’ “distance measure” pi^r by 

poA^^y) ■= (/d’’(a;(p),y(p)) ^p) , 

N 

\dpX - Pdpyf Ax{p)^y{p)\ PA 

c„ eke ]P • 

N ' 

with the usual modification for r = oo. We abbreviate po := /0o,2 and pi := pi,2- Let 
H^(W M) be the completion of C^(A^, M) with respect to po + pi. 

13.4 Lemma. Let ^ be a closed curve in a convex region of M, and let P be the parallel 
transport along 7. If CqL'^A) < tA , then ||P.y — id || < ^CqL^A)- 

Proof. The parallel transport is continuous with respect to L°° convergence in the space 
of loops [ 0 ; 1] —^ M, so it suffices to show the claim for smooth 7. To fix notation, let 
us say 7 : [ 0 ; 1] — >■ M, 7(0) = 7(1) = p. As this curve lies entirely in a convex region, it 
can be represented as 7(t) = exp^ V{t) with a vector field V : [ 0 ; 1] —>■ TpM. Define a 
homotopy c(s,t) := expp sV(t) between 7 and the lazy loop. Denote the s-parameter 
lines by ct and the t-parameter lines by 7^. By |7.8[ we have 

1 1 

P^-id = j j dt ds 

0 0 

The coordinate vectors ct = dsC and 7^ = dtC can be explicitely computed: 9 sc(s, t) = 
Pf'^V{t) because Ct is a geodesic with initial velocity V{t)^ and dtc{s,t) = Jtis) for 
a Jacobi field J* along ct with values Jt( 0 ) = 0 , Jt(l) = j{t) and Jt( 0 ) = V{t) by 
I1.16I (any two of these conditions determine J* uniquely). So we have \ct\ = \V\ and 
1 7s I < ItI because the Jacobi field grows monotonously in s, see the proof of |6.3[ 
Because the parallel transports along 7^ are isometries, we obtain 

iL^-idl ^ J J A\ct\ ITsI < C'oL(7)max|F|. 

And |t^(t)| is the distance fromp to j{t), which cannot be larger than ^^(7), q. e. d. 

13.5 Proposition (“triangle inequality”). If x,y,z G C^(JV,M) with po^oo{x,y) + 
Po.oo(2/, z) < I, then 

Pi{x,z) < pi{x,y) + pi{y,z) + ^Col'^mA^z) 


Pi,r{x,y) ■= [J 
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ig. The Manifold-Valued Dirichlet Problem 


Proof. Pointwise, we have 


Ida; - P^'^dzW < \\dx - P^^^dyW + \\dy - - P^^ypy-^j \\dz\\. 

The difference between the parallel transports is the holonomy of the loop x y 
z X, which is smaller than by Pd -41 Now the claim is a simple application of 

Minkowski’s inequality in L^(M, R), q. e. d. 

Remark, (a) That po^r is indeed a metr i c is p roven with the same argument as the 
usual Minkowsky inequality, see e. g. Alt ( 2006L lemma 1.18). In contrast, [Ig-SI gives 
only a distorted triangle inequality for pi. Nevertheless, H^(A,M) can be defined as 
the completion of C^{N, M) with respect to po + Pi> because this term d oes not disturb 
the usual completion construction for metric spaces, see e. g. IAltI (l2nnllL no. 0.20), and 
we never need to use the triangle inequality explicitely. 

(b) Because of 

|Dir(a;) - Dir(?/)| < pl{x,y) 

(the hidden constant comes from the comparison of \dx — Pdy\ with ||da; — Pdy\), the 
definition above ensures that every u G H^(N,M) indeed has finite Dirichlet energy. 
Nevertheless, not all functions with finite Dirichlet energy are contained in our defini¬ 
tion of H^(N,M), but only those that are limits of smooth function sequences. The 
usual counterexample is the function u 1—^ from the unit ball to the unit sphere min¬ 
imises Dir in dimension m > 3 and larger iIH ildeb randt et aLlliQ77L sec. 6; general 
regularity theory is given in ScHOEN and Uhlenbeck 1082 ). So the usual definition 
of as 

{y G Wi’2(iV^]Rfe); j;(p) g M a. e.}, 

where M is embedded in R^, neccessarily has the drawback that ( A, M) i s not dense 

8^ sec. 4). In 


la 


in M) in dimension 3 and larger (ISchoen and UhlenbeckIL_ _ _ 

allusion Io Imeyers and SerrinI ( iqGtIi . ITost ( iq88L d. 266I states this as H^fA. Ml 
W^’^(A, M). By our use of H^{N, M), we restrict ourselves to functions that can be 
smoothly approximated. This space is well-suited for approximation questions, but the 
wrong one to show existence of solutions. For an overview o ver difficultie s and pitfalls 
of the harmonic mapping problem, we refer to the survey of JostI ( loSSll . 

(c) Consequently, two functions x,y G C^(A, M) are close iff pi(x, y) is finite. 

(d) li x,y G C^{N,M) are close, the geodesics x{p) y{p) give rise to a geodesic 
homotopy h : x y, i.e. & smooth mapping N x [ 0 ; I] —>■ M, {p,s) 1—?► hs{p), such 
that Hq = X and hi = y and s ^ hs (p) is a geodesic for any p. It minimises the energy 


i 

Eih) := J J \^hsip)\^ dsdp 


N 0 


over all homotopies in the same class (|JoST 
p^{x,y) if h is the geodesic homotopy x 


2011 


lemma 8.5.1). In fact, E{h) = 
y and x,y are close, because |-^/is(p)| 


is indepen dent of s in this case. This is also called the L^-width of the geodesic 
homotopy (IKoka rev 201^ the older literature mostly uses the L'^-width po,oo from 
Siegel and WilliamsIIiq8ziI1 . 


13.6 
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C. Applications 


13.7 Proposition (Poincare inequality). Suppose dN is smooth, all Weingarten maps 
of dN with respect to N are bounded by \Wv\ < k everywhere, and no point in N has 
distance larger than r to dN. Then 


Proof. Without regarding the constants, it would be very easy to reduce this case 
to the Poincare inequality for vanishing boundary values i2.ioci As a very personal 
attitude, we would like to circumvent the contradiction argument there. Let us first 
consi der a positive function q : N —» R. 

As IMantegazza and Mennucci ( 20o4 prop. 3.5) have shown, the distance field 
d := d( • ,dN) is except on an (n — 1 )-dimensional set S (in fact, they deal with 
the distance field of an arbitrary submanifold K fo r boundary less N, but the c ase of 
K = dN is also possible). By the coarea formula ( Evans and GariepyI [rooa . thm. 
3.4.2), J g can be computed by integration over the t-level sets N^ := {p G N\S: 
d = t}, where points in S can be omitted because it is a null set: 


/’ = /(/’) 

N 0 m 


dt 


(note that |gradd| = 1 , so there is no additional weighting factor). There are smooth 
homotopies h* retracting each level set to the boundary, defined on a subset Ng of 
dN, with Nq = id and hKN^) = iV*, following the gradient field of d. The intermediate 
mappings h^ cover sets N^ C N^, and the integral can be computed by the funda¬ 
mental theorem of calculus for a{s) = fj^t g as a{t) = a( 0 ) -|- d(s) ds. The derivative 
of the integrals is composed of the integrand’s change along s-lines and the changing 
of the volume element: 


_d 

ds 


J 9 = J MK) + j 

JVJ N* Nj 


gtrWs, 


where hi denotes the s-derivative of the homotopy and Ws is the Weingarten opera¬ 
tor of the distance set A^® from |i.25bl Here we ha ve used that t(s) := trWs is the 
deriva tive of the volume element (IKarcherI IiqSJ . eqn. 1.5.4 or, in a more general 
setup, IDelfour and ZOLESlol 2011I eon. Q.A.17). Now by |i.25b[ the function r(s) 
obeys t < Cq — with initial value r( 0 ) < «: by assumption. This differential inequal¬ 
ity delivers us t < K := max(«:, y/Co). (Note that not the absolute value of r can be 
bounded, only r itself—in fact r —^ — 00 where dhl becomes singular.) So we have 


■fJsilm+Kjs 

AJ AJ AJ 


OT d < b + Ka with b being the integral over \dg\. This differential inequality has the 


96 




















ig. The Manifold-Valued Dirichlet Problem 


supersolution a( 0 )e^‘ + fg b(s) ds, which hence is a bound for a(t). That means 

9 + r J\dg\. 

NO ON N 

Now for / £ H^, let g = /^. The latter term becomes |d(/^)| = 2 /|d/|, and its integral 
is estimated by 2 |/| |(i/| by Holder. Then apply Young’s inequality uv < 5 u^ + 
with 5 = ^ to obtain 

r\d{f)V.<l\f\l+ 2 r^df\l., 

q. e. d. 

Corollary. Situation as before. Suppose x,y G H^(iV, M) are close maps with d(x, y)(j>) 
< e for all boundary points p £ dN. Defining := Cn\/x, it holds po{x,y) ^ 
C'^e + rpi{x,y). 

Proof. Consider the function / := d{x,y) : TV —>■ R. It has differential df{V) = 
glfPy^ {dx — Pdy)V) by [7^ and hence \df\ < \dx — Pdy \, q. e. d. 

Remark, (a) The Poincare inequality in the form above also holds for differential 
forms, with the covariant derivative on the right-hand side. In fact, consider u £ 
and / := |u|. Then, because V is metric, \df\ = {\ 7 u,u)/\u\ < |Vu|. 

(b) By the same method of proof, the Poincare inequality of lXAPPELER et al\ ([200 j . thm. 
0.4) can be significantly shortened. They prove that if N, M are closed and compact and M 
has negative sectional curvature, then any two homotopic mappings x,y £ C^{N,M) satisfy 
Po{x, y) < 1 -f Dir(a:)^'^^ -|- Dir(i/)^/^. 

Situation. For simplicity, we assume N = rK (otherwise, concatenate the results 
below with rTTTTT^l . Suppose the metric 7 of TV is piecewise (i, /i)-small, so that we can 
omit the fullness parameter. If y : TV —>■ TW is a smooth function, we assume that its 
piecewise barycentric interpolation x is close to y, which is the case for small enough 
Coh-^. 

Proposition (Galerkin orthogonality). Situation as in \ig.g\ Let y £ H^(TV, TW) 
he a critical point of Dir with respect to compactly supported variations, and let x £ 
P^(TV, TW) be a critical point of Dir with respect to variations W £ P^Xx as in im.yi that 
vanish at boundary vertices, such that x{pi) = y{pi) on all boundary vertices. Then if 
X and y are close, 

{dx - Pdy, VW) = 0 for all W £ W\aN = 0 . 



Proof. Because x and y are close, the parallel transport induces a bundle isomorphism 
x*TM —>■ y*TM. Because piecewise smooth vector fields are in and the variation 
on the whole boundary vanishes if it vanishes on the vertices (recall that the barycen¬ 
tric mapping on a subsimplex only depends on the vertices of this subsimplex), we 
obtain that PW is an admissible variation field along y for all W £ P^X^. Therefore 
(Pdy, VW) = 0 by |i3.2[ and similarly (dx, VIT) = 0 , q. e. d. 


13.8 


13-9 


13.10 
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C. Applications 


13.11 Corollary. Situation as before. Then \dx — Pdy\\_2 < |da; — Pdy — 

because \dx — Pdy\^ = {dx — Pdy — VW,dx — Pdy) < \dx — Pdy — VW\ \dx — Pdy\. 

Approximation Properties of Karcher Triangulation Variations 

13.12 Lemma. Situation as in \ is.g\ Let V be an vector field along x. Then for any 

i € there is a variation Pi{t) of x(ri) such that the veetor field Xi from [TaT^ 
satisfies \V — < h{l + Coph)\W\\_2(^g.'^ + Cophf\V\\_2(^g.), where Si denotes 

the star of ri, i. e. the union of all simplices rs with i G 5. 

Proof. Abbreviate p := x{ri) and write X instead of Xi for the time of the proof. In 
the first step, let us consider a smooth vector fields V. Choose the variation of x{ri) 
such that p(0) = V\p. Then by 112.51 X = V at p and hence 

(V-A)l, = J PV^{V-X) = I py^v- J py^x, 

1 77 

where j : p q. The second integral should disappear in the result. By |i2.4( we have 
l^7"^l ^ |p|. So we end up with |p|, which is a point evaluation of V and 

hence undesired. Express p = V\p = PV\.y(t) — f where the integral only runs 

from 0 to t. Then 

< I ||VV||+Co,i/i I {\V\+I liVVil) <(l + Co,i/i) I ||VE||+Co,i/i J |I/|, 

7 7 

Squaring both sides and applying Holder’s inequality as in [7^5^ gives 

J\V- Ap < hil + Coph) J llVEp + Coph^ J ivp. 

Si Si Si 

So for a smooth vector field V, we have constructed an interpolation. The best ap¬ 
proximation in must of course also fulfill this inequality. And by continuity of the 
L^-orthogonal projection, this holds for every vector field of class H^, q. e. d. 

13.13 Proposition. Situation as in \ig.g\ LetQ he an section ofT*N®x*TM. Then there 
is some W G P^'X.x with jQ —VIT|l2 ^ /i|V(3|l2 + Coph?'\Q\. The hidden constant 
depends on n, m and |A|. 

Proof. It suffices to show the claim for the L^(iV, M) operator norm in the left-hand 
side instead of norm: 


||Q-VIT||l2 </i|VQ|l2+Co,i/i'|Q|. (i3-i3a) 

In fact, let v be the unit vector field on N realising IQl everywhere. Then v € 
L'^{x*TM) and hence |Qp = J \Q\^ < / IQIP = J \Qvf < < llQlPl^P- 

So let us prove [T3A3^ 


98 
















ig. The Manifold-Valued Dirichlet Problem 


In any simplex, Q can be applied to vectors ri — rj and their linear combinations. 

Choose a norm-preserving extension of Q such that it can also be applied to vectors 
ri. Then define, on each star Si, a vector field Vi := Qri. Then Qv = v^^Vi for any 
V S TN\s^. Now let Xi be the best approximation to —Vi on Si. Then 

( j\Qv+v^x,\^y^=(^ I iv^)^\v,+x,\^y^ 

< \vmv + x,\ 

< \v^\{h{l + Co.ih)\VV\+Co.ih^\V\) 

< \v\{h{l + Co.ih)\\XQ\\ + Co.ih^WQW). 

Now recall that \A\Vdxvi + Xi\ < Co^ih\dxv\ \x\ from [T2T^ in combination with |6.i4| 
and |i2.4| and || id || < (7nh^ from ifi.ifil This efives \X dxvX P v'’Xi\ ^ Coh"^\Xdxvi\ + 

\Xi \ So we have for W := x £ 

- Qv\ < IVdxvW + v^X,\ + \v^X, + Qv\ 

<Co,ih^\v\\\VW\\+h{l + Co.ih)\\XQ\\+Co,ih^ ||Q||. 

Because VIT is almost an best approximation of Q, we have ||V 1 T|| < IQ I, which 
completes the proof, q. e. d. 

Theorem. Situation as in \ig.lO\ Then 13-14 

Po{x,y) -f pi{x,y) < po^oN{x,y) -f h\Xdy\ -f Cq 


where the hidden constant depends on m and the geometry of N. 

Proof. Applying | r3.1 31 to Q = dx — Pdy, there is IT G P^Xa, with 

\dx — Pdy — XW\ < h\Xdx — XPdy\ Co,ih‘^\dx — Pdy\. 


Because C^gh? is assumed to be small, say < i, we can neglect the latter term. Due 
to I6.22I I da; I ;< Cq and VPdy — PVdy can be shown to be bounded with an 

argument like in |7.12| (spelled out in detail, this amounts to a rought estimate 
for d(a;, y), w hich is provabl e by a suitable modification of the standard first-order 
estimate as in BRAESsll2nn7l p. 89). Then the claim is proven by|i3.ii|and|i3.8] q. e. d. 
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D. Outlook 


There are several research directions that would naturally continue the course of this 
dissertation, but which could not be further investigated due to time constraints: 


(1) Whereas the weak formulation and approximation of extrinsic curvature is obvi- 
o usly bound to the embedding of a submanifold, the weak form of Ricci curvature as 
in iFRiTZ (1201*^) could be formulated intrinsically. 

(2) The measure-valued ap proximation of Lipschitz-Killing curva tures of submani¬ 
folds in Euclidean space from Cohen-Steiner and MoRVAN ( 2nn(lh could possibly be 
carried over to situations where the surrounding space itself has curvature. 

200-:!|l 


OSHER and Fedkiw 


(3) The level set approach (jOsHER and SethianiI 
that was used to approximate pde’s o n surfaces ( PznJK and Elliott! 2008 !. surface 
flows ( Deckel nick and Dziuiq 2001! or, as combination of both, pde’s on evolving 


surfaces ( PziUK and Elliott! "i 4 sec. 8 ) can directly be carried over to submani¬ 
folds. 

(4) Assumption |9.2oa| has to be verified, perhaps under additional conditions. The 
testing with forms in lio.2611251 should be sharpened or at least re-interpreted with 

the use of more classical test functions. 

{5) Our definition of the barycentric mapping x is implicit and needs to know gradi¬ 
ents of the squared distance function d^. The exact computation of geodesic distances 
is very expensive, a nd the task to find fast and accurate approx imations is a current re¬ 
search problem, cf. Crane et al. 1 20i‘jl . Campen et al. 1 20iqll and references therein. 
The use of any of these approximations to compute the barycentric mapping would 
lead to a computationally feasible approximation of x. 

(6) After this, or restricted to 3-manifolds where geodesic distances can be exactly 
computed (or s ufficiently well approximated), t he minimal surface algorithms from 
BR.AKKEj(lrQQ2h.lPlNKALL and PoLTHl^ (jroQ-jl . lRENKA and NEUBERGERldiQQtjl . and 
Dziuk and HutchinsonH iqqcJI can be applied, for example in hyperbolic three-space 


the product x IR of hyperbolic 2-space and the real line, or products with 
twisted metrics. 

(7) V ariational methods in shape space, as have been dealt bv iRuMPF and WirthI 


( 20111) . can be extended e. g. to the computation of minimal submanifolds (whose 
dimension can be freely chosen) or multi-dimensional regression. 
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ZUSAMMENFASSUNG 


(i) Sei {M,g) eine unberandete, kompakte Riemannsche Mannigfaltigkeit und A das 
n-d imensionale Stand ardsimplex. Fiir n + 1 gegebene Punkte pi G M betrachten wir 
mit IKarch^ (|iQ77h die Funktion 

FJ : M X A R, (a, A) i-A A° d^(a,po) H-+ A" d^(a,p„), 


worin d der geodatische Abstand in M sei. Liegen alle pi in einem hinreichend kleinen 
geodatischen Ball, so ist a: : A i—>■ argmin^ E[a, A) eine wohldefinierte Funktion A —^ M 
( }gTgj ). Wir nennen s := a;(A) das Karcher-Simplex mit Ecken pi. Auf A sei eine flache 
Riemannsche Metrik durch Vorgabe von Seitenlangen A{j)i,pj) definiert. Wenn alle 
Seitenlangen kleiner als h sind und vol(A,5®) > aK^ fiir ein a > 0 ist, so zeigen wir 
in I6.17I und |6.23| 

\{x*g — g^){v,w)\ < ch?\v\\w\, |(V^ ® — V® < c/i|u| |r(;| (A.ia) 


mit einer nur vom Kriimmungstensor R von (M, g) und abhangigen Konstanten c. 
Daraus folgen mit wenig Aufwand Interpolationsabschatzungen fiir Funktionen u : s ^ 
R ([7|4j) und y : s —>■ A fiir eine zweite Riemannsche Mannigfaltigkeit N ( 17.15} ). Auch 
erlaubt diese Simplexdefinition, auf Grundlage der Voronoi-Zerlegung von Leibon und 
Letscher (2000) eine Karcher-Delaunay-Triangulierung zu definieren ( [8.7^ . 

Daher konnen wir im folgenden ganz (M, g) als trianguliert annehmen. Auf jedem 
Simplex ist g durch eine Metrik y® mit lA.iai approximiert, und schwach differenzier- 
bare Funktion u G H^(M,g) lassen sic h durch stiic kweise polynomielle Uh G P^{M) 


approximieren. In der iiblichen Weise ( Dziuid iq88L Holst und Stern 2012) lassen 


sich daher Variationsprobleme wie das Poissonproblem ( jio.i^i 110.171113.14! ) oder die 
Hodge-Zerlegung ( jld.isj ) in H^(M, y) mit denjenigen in H^(M, y®) und ihren Galer- 
kin-Approximationen in P^{M) vergleichen. 

Ankniipfend an die gangige Finite-Elemente-Theorie fiir Probleme auf Untermannig- 
faltigkeiten des R™ lassen sich auch Untermannigfaltigkeiten S C M durch Karcher- 
Simplexe approximieren. Der dabei auftretende Geometriefehler ist gleich dem fiir 
Untermannigfaltigkeiten des R"* zuziiglich eines Terms ch^ (Iii.i8|) . 


(2) Sei M die geometrische Realisierung eines simplizialen Komplexes A . Die sim- 
pliziale Kohomologie (C'^(A), d*) ist von Desbrun und Hirani (|2nn4l2nr)^ als diskre- 
tes auheres Kalkiil (dec) interpretiert worden. Wir definieren Raume P~^D^ C 
und aufiere Differentiale und geben eine isometrische Kokettenabbildung P^^D^ 

an ( jp.iij ). Damit ist die Berechnung von Variationsproblemen im diskreten aufeeren 
Kalkiil auf Variationsprobleme in einem Raum von nicht-konformen Ansatz-Differ- 
entialformen zuriickgefiihrt. Wir untersuchen die Approximationseigenschaften von 
in 1 19.191 |9.2oj ) und vergleichen die Losungen von Variationsproblemen 

in ihnen (!io.26l[25|) . 

































